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Preface 



This Instructor's Solutions Manual is an ancillary for the fifth edition of Grossman's 
Elementary Linear Algebra. It contains detailed solutions to all problems in the text- 
including die MATLAB and graphing calculator problems — and in the Applications 
Supplement. Below is an overview of all the ancillaries to accompany the main text. 



Applications • one chapter each on linear programming and on Markov 

Supplement chains and game theory 

• available packaged with the text or for separate purchase 

• numerous examples and problems 

• answers to odd-numbered problems are at the back of the 
Applications Supplement 



Student Solutions • complete solutions to all the odd-numbered problems in the 

Manual text and the Applications Supplement 



MATLAB Manual: • computer laboratory exercises and applications using 

Computer Laboratory MATLAB. Each section lists objectives, prerequisites, 
Exercises and M-file and MATLAB features before the lab exercise is presented, 
disk The student is then encouraged to apply concepts 

interactively and create an edited diary session. An M-file 
disk containing programs of selected applications in the 
manual is available free upon request from The MathWorks, 
Inc., in either Mac or PC versions. 



Elementary Linear • MATLAB programs that accompany the main text in either 

Algebra Toolbox PC or Mac version are available free upon request from The 

(M-file disk) MathWorks, Inc. 



HP- 48G/GX • calculator enhancement for science and engineering 

Calculator Manual mathematics using the high-level Hewlett-Packard calculator 
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Chapter 1. Systems of Linear Equations and Matrices 

Section 1.2 

1. 4x — 12 y = 16 

—4x + 2y — 6 

— 10?/ = 22 hence, y = —11/5 and x = —13/5. 

a n a 22 — ^12021 = (1)(2) — (— 3)(— 4) = 2 — 12 = —10. 



2. 14a; — 7y = —21 

5x + 7y = 4 

19a; = -17 hence, x = -17/19 and y = 3 + 2x = 3 - 34/19 = 23/19. 

<■11-22 - 012021 = (2)(7) - (— 1)(5) = 14 - (-5) = 19. 



3. Qx — 24 y = 15 

— 6x + 24y = 16 

0 = 31 => no solution. 



-n-22 - ai2 a 2i = (2)(12) - (— 8)(— 3) = 24 - 24 = 0. 



4. 6x — 24 y = 18 

— 6x + 24?/ = — 18 

0 = 0 =£■ lines coincide. 

a u a 22 — a i2 a 2i — (2)(12) — (— 8)(— 3) = 24 — 24 = 0. 



5. 6* -f y = 3 

—4a; — y = 8 

2a; = 11 hence, x — 11/2 and y — 3 — 6a; = 3 — 33 = —30. 

011022 — 012021 = (6)( — 1) — (1)(— 4) = — 6 — (—4) = —2. 



6. 9a; + 3y = 0 

2a; — 3?/ = 0 

llx =0 hence, x = 0 and y = —3a; = 0. 

a n a 22 — «i2°2i = (3)( — 3) — (1)(2) = — 9 — 2 = —11. 



7. 4a; — 6 y =0 

—4x + 6?/ = 0 

0=0 => lines coincide; 4a; — 6y = 0 implies y — (2/3)x for arbitrary x. 
a n a 22 — 012021 — (4)(3) — (— 6)(— 2) = 12 — 12 = 0. 

8. 25x + 10 y = 15 

4x + 10?/ = 6 

21a; = 9 hence, x = 9/21 and y = (3 — 5x)/2 = 9/21. 

aiia 2 2 - 012021 = (5)(5) - (2)(2) = 25 - 4 = 21. 
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9. 8x + I2y = 16 

9a: 4- 121/ — 15 

—x =1 hence, x — — 1 and y = (4 — 2x)/3 = 2. 

aii«22 — <*12021 = (2)(4) — (3)(3) = 8 — 9 = —1. 

10. ax + by — c 

ax — by — c 

2 ax = 2c hence, x = c/ a (assuming a ^ 0) and y = (c — ax)/b = 0 

(assuming 6^0). If a = 0 and 6^0, then there are no solutions unless c = 0, in which case y — 0 
and any a: is a solution. If a ^ 0 and 6 = 0, then x = c/a and any y is a solution. Finally, if both a 
and 6 are zero, then there are no solutions unless c = 0, too, in which case any x and y gives a solu- 
tion. au a 22 — ai 202 i — o(— 6) —ba = —2ab. 

11. a 2 x + aby = ac 

6 2 z + aby = be 

(a 2 — 6 2 )x = ac — be hence, x = c(a — 6)/(a 2 — 6 2 ) = c/(a + 6) 

(assuming a 2 — b 2 ^ 0) and y = (c — ax)/b = c/(a + 6) also. If a 2 — b 2 = 0, then a = ±6; if a = 6 ^ 0, 
then the equations are the same, and y = (c/a) — x and any x gives a solution. If a = —6, then there 
are no solutions unless c = 0, in which case any x and y give a solution if 6 = 0, and if 6 ^ 0, y = x 
with any x gives a solution. 011022 — 0 i 2 a 2 i = aa — 66 = a 2 — 6 2 . 

12. a 2 x — aby = ac 

b 2 x -f aby = bd 

(a 2 -F b 2 )x = ac + bc hence, x = (ac + bd)/(a 2 + 6 2 ) 

(assuming a ^ 0 and 6^0) and y = (d — bx)/a = ad — be. 

ana 22 — 012021 — aa — (—6)6 = a 2 + 6 2 . 

13. We need —ab — ab = — 2a6 0. Therefore, we need a ^ 0 and 6^0. 

14. We need a 2 — 6 2 = 0. Therefore, a = 6 or a = —6. If a = 6, then c can be any real number; if a = —6, 
then only c = 0 gives a solution. 

15. We need a 2 + 6 2 = 0. Therefore, a — 0 and 6 = 0. We would also need either c or d to be non-zero. 

16. 3x — 3j/ = 21 

2x + 3y = 1 

bx = 22 hence, x = 22/5 and y = (1 — 2x)/2 = —13/5. 
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17. —4* + 2y = 8 

4* — 2y = 6 

0 = 14 => no point of intersection. 

18. 12* - 18y = 21 

12* - 18v = 24 

0 = — 3 => no point of intersection. 



19. 12* - 18y = 30 

12* - 18v = 30 

0 = 0 =>• lines coincide. 



20. 3* + y = 4 

—5* + y = 2 

8* =2 hence, * = 1/4 and y = 4 — 3* = 13/4. 

21. 6* + 8y = 10 

6* — 7y = 8 

15 y = 2 hence, y = 2/15 and * = (5 — 4y)/3 = 67/45. 

22. Let mi = the slope of L and m 2 = the slope of Lx- 
mi = 1; m 2 = —1. L : x — y = 6, and Lx : * -f y = 0 

Point of intersection: (3, —3). d = ^/(3 — 0) 2 + (—3 — 0) 2 = 8y/2 



23. mi = —2/3; m 2 = 3/2. L : 2* + 3j/ = —1, and Lx : 2* + 3j/ = 0 
Point of intersection: (—2/13, —3/13) 
d = 7(— 2/13 - 0) 2 + (-3/13 - 0) 2 = 0713 



24. mi = —3; m 2 = 1/3. L : 3* + y = 7, and Lx : * — 3j/ = — 5 
Point of intersection: (8/5, 11/5) 
d= 08/5- 1) 2 + (11/5- 2) 2 = 0/5 



25. mi = 5/6; m 2 = —6/5. L : 5* - 6y = 3, and Lx : 6* + 5y = 28 
Point of intersection: (3, 2) 

d= 03 - 2) 2 + (2 - 16/5) 2 = 01725 

26. mi = 5/2; m 2 = —2/5. L : —5* + 2y = —2, and Lx : 2* + 5y = —5 
Point of intersection: (0,-1) 

d= 05 — 0) 2 + (-3+ l) 2 = 09 




4 Chapter 1 Systems of Linear Equations and Matrices 



Instructor's Manual 



27. mi = —1/2; m 2 = 2. L : 3x + 6y = 3, and : 2x — t/ = 17 
Point of intersection: (7, —3) 

d= ^(8-7)2 + (-l + 3) 2 = 75 

28. 4x — 6y = 2 

3x + 6v = 12 

7x = 14 hence, x = 2 and y = (12 — 3x)/6 = 1. 

Then, we need the distance between (2, 1) and 2x — y = 6. 
mi = 2; m 2 = —1/2. L : 2x - y = 6, and Z*i : x + 2y = 4 
Point of intersection: (16/5,2/5) 
d = v/(6/5 - 2) 2 + (2/5 - l) 2 = 3v/5/5 



29. mi = —a/6; m 2 = 6/a. L : ax + by = c, and Li : ax + by = 6xj — ayi 

fac + Pxi - abyi be - a6xj + afy x \ 

Point ot intersection: ^ r= , 5 rr I 

V a z + b* a* + b* J 



Then d = |axi + by\ — c|/\/a 2 + 6 2 after a lot of algebra. 



30. Let x = number of birds and y = number of beasts. 

Then x + y = 60; 

2x + 4y = 200. Hence x = 20 and t/ = 40. 

31. If ana 22 — ai 2 a 2 i = 0, then ana 22 = ai 2 a 2 i. Assuming ai 2 a 22 ^4 0, and dividing both sides of the 
equation by ai 2 a 22 we get an/ai 2 = a 2 i/a 22 . This implies — au/ai 2 = — a 2 i/a 22 . Solving each linear 
equation of system (1) for y we get y = — an/ai 2 x + bi/an and y = — a 2 i/a 22 x + 6 2 /a 22 . These lines 
have slope — au/ai 2 and — a 2 i/a 22 respectively. Slopes are equal. Therefore the lines are parallel. If 
ai 2 = 0, then from ana 22 = 0 and an ^4 0, we get a 22 = 0. So the lines are parallel because they are 
both vertical. If a 22 = 0 similar reasoning holds. 



32. Suppose otherwise, i.e. , suppose that ana 22 — ai 2 a 2 i = 0. Then #31 shows that the lines given in sys- 
tem (1) are parallel. Thus system (1) either has an infinite number of solutions or no solution. This 
contradicts the assumption that the system has a unique solution. Result follows. 

33. If ana 22 — ai 2 a 2 i # 0 then ana 22 # ai 2 a 2 i. Dividing both sides of the equation by ai 2 a 22 we get 
«n/«i 2 7^ a 2 i/t* 22 - Thus — an/ai 2 # — a 2 i/a 22 . Hence the slopes of the lines in system (1) are not 
equal (see solution to #31). Thus the lines are not parallel. Therefore system (1) has a unique solu- 
tion. 
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34. Let x — number of cups and y = number of saucers. 

Then, 3* + 2 y — 480 Eq. 2 —10 Eq. 1: —5* = —400 

25* + 20j / = 4400; 
hence, x = 80 and y = 120. 



35. 3* + 2 y= 480 Eq. 2 -5 Eq. 1:0 = 0 

15* + lOy = 2400. 

So, these two lines coincide; hence y = (480 — 3*)/2 where 0 < * < 160, to force y > 0. 



36. 3* + 2 y- 480. Eq. 2 -5 Eq. 1: 0 = 100, 

15* + 10y = 2500 

so this system of equations has no solution. 



37. Let * = number of ice-cream sodas and y = number of milk shakes. 

Then, * + y = 160; Eq. 2 —4 Eq. 1: — y = —128. 

4* + 3y = 512 

Hence, y = 128, * = 160 — y = 32. 
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Section 1.3 

/ 1 -2 3 11 \ R 2 -+- 4 R 1 +R 2 / 1 -2 3 

1. 4 1 -1 4 R 3 - - 1 R 1 + R 3 0 9 -13 

\2 -13 10/ * \0 3 -3 

, , , r 3 -+ —\r 3 

R '-“’ +R ' ;| 3 \ *,-.*+*, 

R 3 — ♦ — 9 R\ R 3 I 0 1 -1 



0 0-4 



-40 R 2 
- 12 / . 




R 2 ^R 3 + R 2 1% 3 

Ri — *■ —R 3 4" Rl V 0 0 1 1 



/I 


-2 


3 




0 


1 


-1 


-4 


\0 


9 


-13 


-40 / 



(2, —3, 1) is the unique solution. 



(-2 1 6 18 \ R 3 -+Ri+R 3 / I 3-4 15 \ Ba - -5fli + fla / 1 3 -4 15 \ 

2. 5 0 8 -16 } Ri^Rs 15 0 8 -16 fi 3 -2fii+R3 0 -15 28 -91 

\ 3 2 -10 -3/ „ \-2 16 18/ + \0 7 -2 48/ 

R 2 — *• 2 B 3 + R 2 t\ 3 .4 15 \ B 3 -- 7 B 2 + B 3 /1 3 -4 15 \ 

r 2 R 3 [ 0 1 -24 -5 I R 3 -* j j$R 3 ( 0 1 -24 -5 1 . Use back substitution to 

, \0 7 -2 48/ „ \0 0 1 0.5/ 

find the solution (— 4, 7, |). 



3 6-6 

2-5 4 

-1 16 -14 



Ri -* f Hi 

R2-.-2R 1+ R2 f X Q l * 

R 3 -+R 1 + R 3 \o 18 -16 0 



R 3 - 2R 2 + R 3 / 1 2 -2 3 

r 2 ^-Lr 2 0 1 -8/9 0 

\ 0 0 0 0 



be arbitrary. Use back substitution to find the solutions (3 + §2:3, §x 3 , x 3 ). 



Ri - k R i 

/3 6 -6 9\ R ’ H (l 2 -2 3' 

4. 2 -5 4 6 + 0-9 8 0 

1 5 28 -26 -8/ 5R i + R 3 \ 0 18 -16 -23, 



R 3 — ► 2i?2 + R 3 



12-2 3\ 

0 -9 8 0 I . The bot- 

0 0 0 -23/ 



tom row is equivalent to the equation 0 = 23, which is impossible. So the system has no solution. 



/I 1-1 7 \ «a 

5. 4 -1 5 4 R 3 

l 2 2 -3 0/ 



-4R!+R2 / 1 1-1 7' 

- 2 Ri + R 3 [ 0 -5 9 -24 
» lo 0-1 -14, 



R 3 — ► — i ?3 

— *■ — 9H3 -|- J?2 
Ri —*■ R 3 + Ri 



r 2 ^-Ir 2 /i 0 0 -9\ 

/ti -► -« 2 + Ri 0 1 0 30 I . (-9,30, 14) is the unique solution. 
» \o 0 1 14/ 



110 21' 
0-50 -150 
0 0 1 14, 



/ 1 1-1 7\ «2 

6. 4 -1 5 4 1 R 3 

1 6 1 3 18/ 



-4fii +B 2 / 1 1-1 7' 

-6fl! + i?3 10-5 9 -24 
„ \ 0 -5 9 -24, 



R3 —R2 + R -3 

r 2 - -\r 2 

R\ — ► — /?2 "I" ^1 



1 0 4/5 12/5 
0 1 -9/5 24/5 
0 0 0 0 



. Let 2:3 



be arbitrary. Then (y- - f 2:3, ^ + §x 3 ,x 3 ) are the solutions. 
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7 . 




10 4/5 
0 1 -9/5 
0 0 0 



7 

4 J . The same row operations as in problem 6 gives the equivalent matrix 

20 , 

O' 

0 J . Since 0^1, the system has no solution. 

1 




i?2 — * — 4i?i -f R 2 
R3 — ► — 2 Hi + ^3 



1 


-2 


3 


ft? 

t 

P5 

* 

0 


/ 1 -2 


3 


0 


9 


-13 


0 r 2 ^ r 3 


0 1 


-1 


0 


3 


-3 


0/ — , 


Vo 9 


-13 



. ( 0 , 0 , 0 ) is the unique solution, by back substitution. 



Ri — ► 2R 2 -f- Ri 



(1 1 -1 


oN ! 


R 2 —*• —4 Ri 4" R 2 


(1 1 -1 


0 \ R3 — ♦ —Ri + R3 


9. 4 -1 5 


0 


R 3 ► — 6 R 1 4" R 3 


0-59 


0 r 2 - -i r 2 


\ 6 1 3 


0) 




^0 -5 9 


0 



trary. Use back substitution to find the solutions (— |x 3 , | x^ 3 , x 3 ). 




. Let X 3 be arbi- 



10 . 



R\ +Z Ri 
R 3 — ► —2R\ + R 3 





R 2 -» \R 2 



R 3 



/ 10 -2 
4 Ri + Rs 0 1 5/2 

> 0 0 -6 



34/3 \ 

-37/6 I. (34/3, —37/6, 11/3) is the unique solution. 
11/3/ 



11 . 



(1 2 -1 4\ 
V3 4 -2 7 J 



R 2 — ► —3R\ + ^2 
Ri — 1 ► — \Ri 



f 1 2 - 1 4N ) 

Vo 1 -1/2 5/2 J 



Ri — * —2R 2 4* ^1 



(i 0 
Vo 1 - 1 . 



be arbitrary. Then (—1, | + 1 x 3 , 13 ) are the solutions. 



0 

1.2 



5 / 5 / 2 )' LetXa 



12 . 



/ 1 2 -4 4\ 

V-2 -4 8 -8 J 



R 2 — ► 2R\ 4 R 2 



0 



2 -4 
0 0 



(4 — 2x2 + 4x3, X2, X3) are the solutions. 



4\ 

q J ■ Let X 2 and x 3 be arbitrary. Then 



13. 



/ 1 2 -4 4\ 

V-2 -4 8 -9/ 



4 \ R 2 -*2R 1 + R 2 l \ 2 -4 
-9/ » VO 0 0 



^ 2 ' 4 M Si 
VO 0 0 -1/ 



Since 0^—1, the system has no solution. 



14. 



12-11 
3 6-33 



Ri —* —3Ri + Ri ( 1 2-1 1 
0 0 0 0 



(7 — 2x2 + x 3 — X4, X2, X3, X4) are the solutions. 



. Let X 2 , X 3 , and X 4 be arbitrary. Then 
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15. 



2 6-42 
10-11 
\-3 2-2 0 

R 3 —* _ -ns^3 

R-2 — * 3 ^3 + R-2 

Ri — ► Rs ■+■ ^1 



Ri — * |jRi 

R 2 — ► — Ri + R 2 
R:\ — * 3 Ri + R'2 



1 3-2 1 

0-310 
0 11 -8 3 



H 2 -+ - \R2 
R \ — ► — 3 R 2 + ^1 
R 3 — ► — \\R 2 + R 3 




-1 1 
-1/3 0 
-13/3 3 



1 0 0 4/13 
0 1 0 -3/13 
0 0 1 -9/13 



20/13 \ 
-28/13 
-45/13 



• Then (ft - A X 4,-f§ + ^ 4 , -if + J3 X *’ X *) are the 



solutions, where x 4 is arbitrary. 



16. 



/ 1 -2 1 1 


2 V 


R 2 — *■ —3Ri 4- R 2 ( l i 1 1 




R 2 R 3 


/I -2 1 1 


2 \ 


3 0 2 -2 


-8 


„ „ 0 6-1-5 


-14 




0 1 -0.25 -0.25 


0.25 


0 4-1-1 


1 


Ri-tRi+Ri 1 0 4 


1 


R 2 — *• \R2 


06 -1 -5 


-14 


V-l 6-2 0 


7/ 


+ Vo 4 -1 1 


9/ 




Vo 4 -1 1 


9/ 



«3 

Ri 



- 6 R 2 4* R 3 
-4 R2 + Ri 



/I 


-2 


1 


1 


2 V 


R 3 — ► 2i?3 


/I 


-2 


1 


1 


2 \ 


0 


1 


-0.25 


-0.25 


0.25 




0 


1 


-0.25 


-0.25 


0.25 


0 


0 


0.5 


-3.5 


-15.5 


Ri — ► | Ri 


0 


0 


1 


-7 


-31 


Vo 


0 


0 


2 


8/ 




Vo 


0 


0 


1 


4/ 



. Use back sub- 



stitution to find the solution (2, 1/2, —3,4). 



17. 



z 1 

3 



2 1 1 
0 2-2 
0 4-1-1 
\5 0 3 -1 



J' 

1 

-3/ 



As in 

problem 16 



(\ -2 1 1 
0 1 -0.25 -0.25 



0 6 

Vo 10 



-1 

-2 



-5 

-6 



2 \ 

0.25 

-!4 

-13/ 



R 3 —*■ — 6J?2 + R 3 
Ri — ► — IOR 2 4 “ Ri 



/ 1 


-2 


1 


1 


2 ^ 


Ri —* —R3 4- Ri 

n , 0 n 


/I 


-2 


1 


1 


2 \ 


0 


1 


-0.25 


-0.25 


0.25 


0 


1 


-0.25 


-0.25 


0.25 


0 


0 


0.5 


-3.5 


-15.5 I 


/13 — ♦ 2/t3 


0 


0 


1 


-7 


-31 


Vo 


0 


0.5 


-3.5 


-15.5/ 




Vo 


0 


0 


0 


0/ 



. Let X 4 be arbitrary. Use 



back substitution to find the solutions (18 — 4 x 4 , —15/2 + 2 x 4 , —31 + 7x4, X4). 



18. 



( l 


-2 


1 1 


2 \ 


As in 


( l 


-2 


1 


1 


2 \ 


I 3 


0 


2 -2 


-8 


problem 17 


0 


1 


-0.25 


-0.25 


0.25 


° 


4 


-1 -1 


1 


0 


0 


0.5 


-3.5 


-15.5 


\5 


0 


3 -1 


0/ 




Vo 


0 


0 


0 


3/ 



(Just apply the row operations 



above to the changed last column). Since 0^3, the system has no solution. 



19. 





R 2 — ► “ 2R\ -f R 2 


f 1 1 


4 \ 


R 3 —* —3Ri + R 3 | 


f 0 -5 
Vo -1 


V 



R '2 R 3 
/?2 — ► —R 2 
R3 — * 5R2 + R3 




4 \ 

4 I . Since 0 ^ 19, the system 



has no solution. 



( l 1 
20. 2 -3 

\3 -2 




As in 

problem 19 



Z 1 1 

0 -5 
\0 -5 




find the solution (19/5, 1/5). 



R 3 —R 2 4- R 3 
R 2 —*■ — g-^2 



/! 1 4\ 

I 0 1 0.2 j . Use back substitution to 

\0 0 0/ 
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21. row echelon form 22. neither as the first nonzero in row 1 is not a 1. 

23. reduced row echelon form 

24. reduced row echelon form 

25. neither as the first nonzero in row 2 is too far left. 

26. reduced row echelon form 
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1 -7/2 \ / 1 — 7/2 \ 

0 31/2 I — ► [ 0 1 j row echelon form 

0 11 / \0 0 / 

/ 10 \ 

reduced row echelon form 



36. We have n = 3, 1 - a n = §, 1 - a 22 = f , and 1 - a 33 = §. Then the system is 



Using row reduction, we obtain 



' 2/3 -1/2 -1/6 
-1/4 3/4 -1/8 
^ -1/12 -1/3 5/6 

TO -4 
0 1 -3/2 
0 0 21/2 



1320/7 \ 
-300/7 
4800/7 / 




2 1 1 in 

3*! "S' 2 "6*3 =1° 

1 3 1 

-4* 1 + 4*a-8* s =l 5 

1 1 . 5 on 

-i2 IX -f 2+ 6 l3 = 30 



-360 \ 

-600 -v 
1500/ 

3560/49 \ 

2700/49 . Hence, the outputs needed for supply to equal 

3200/49 / 






demand are x\ — 3560/49, xi = 2700/49, and X3 = 3200/49. 



37. As in example 10, we have the following system: 

xi + 3x2 + 2x3 = 15, 000 

xi + 4 x 2 + X3 = 10, 000 

2xi + 5 x 2 + 5x3 = 35, 000. 

Writing an augmented matrix for the system and finding the reduced echelon form, we obtain 




15000 \ 
10000 
35000 / 




30000 \ 

-5000 I . Since xi, X2, and X3 must be 
0/ 



15000 \ /I 0 5 

-5000 0 1 -1 

5000/ \0 0 0 

greater than or equal to zero, we must have xi = —5x3 + 30,000 > 0 and X2 = X3 — 5000 > 0. 
Hence, the populations that can be supported are 5,000 < X3 < 6,000, xi = 30,000 - 5x 3 , and 
x 2 = —5, 000 -I- x 3 . The solution is not unique. 



38. Let d E d F , and d s denote the number of days spent in the respective countries. The information 
gives the following system of equations: 

30Te + 20 d F + 20^5 = 340 
20d E + 30 d F + 20ds = 320 
10d£ + lOdp + 10ds = 140 
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Upon reducing, we obtain 

/ 30 20 20 340 \ /l 1 1 14\ /ill 14\ /l 0 1 10\ /l 0 0 6\ 

20 30 20 320 U 2 3 2 32 U 0 1 0 4 I -» [ 0 1 0 4 -> 0 1 0 4 . 

V 10 10 10 140/ \3 2 2 34/ \0 -1-1-8/ \0 0 -1 -4/ \0 0 1 4/ 

Hence, dp — 6, dp = 4, and d$ = 4. 



39. Let sp , S [[ , and sm denote the respective number of shares. The information gives the following sys- 
tem of equations: 

—sp — 1.5 sp -f 0.5 sm = —350 
1.5s# — 0.5sjy + sm — 600 

Writing the system as an augmented matrix and reducing to echelon form gives 
/-I -1.5 0.5 -350 \ ?1 1.5 -0.5 350 N / 1 0 0.4545 390.9 N c . 

\ 1.5 -0.5 1 600 J ^0 -2.75 1.75 75/ 1^0 1 -0.6364 -27.27 )' Since s m can be 

chosen arbitrari y, the broker does not have enough information. If sm = 200, then sp = 300 and 

sp = 100 . 

40. Let / and b denote the number of fighter planes and bombers, respectively. The information gives the 
following equations: 

/ + b = 60 
6/ + 26 = 250 
/ - 26 = 0 

/I 1 60 \ /I 1 60 \ /I 1 60 \ /I 1 60 \ 

Then 6 2 250 ) — ► [ 0 -4 -110 -+ 0 1 20 } — » ( 0 1 20 . Since 0 ^ -30, the 

\ 1 -2 0/ \0 -3 -60/ \0 -4 -110/ \o 0-30/ 

system is inconsistent. 



/ 2 -1 3 a\ / 1 2-8 b-a\ / 1 2 -8 b-a\ „ (l 2 -8 b-a\ 

41. | 3 1-5 b ] — ► [ 2 -1 3 a ] — *■ ( 0 -5 19 -2b+3a J 3 2 ( 0 -5 19 -2b+3a |. 

\-5 -5 21 c/ \-5 -5 21 c/ \0 5 -19 5b-5a+c / + \0 0 0 3b-2a+c/ 

Hence, the system is inconsistent if 36 — 2a + c ^ 0. 



/ 2 3 -1 a\ /I -13 b\ _^ R _ m (l -1 3 b\ 

42. ( 1 -1 3 b ) -f | 0 5 -7 a-2b ) 3 3 2 ( 0 5 -7 a-2b ) . For the system to 

\3 7-5 c/ \0 10 -14 c-3b / + \0 0 0 -2a+b+c / 

be consistent, we must have —2 a -f 6 -f c = 0. 
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43 . Either an, a2i, or 031 is nonzero, otherwise, the system is either inconsistent or has an infinite num- 
ber of solutions. Without loss of generality, we may assume an ^ 0 . Elementary row operations 

( an ai2 a 13 6j \ / 1 an/an a 13 /a n *\ 

a 21 a 22 a 23 &2 I — *■ jo a?22 023 * I where 022 = a 22 ~ a 21 a 12/ a llj O23 = 

“31 “32 033 &3 / \ 0 032 O33 * / 

023 — 021013/011, 032 = 032 — 031012/011, and 033 = 033 — 031013/011- As before, either 022 

( 1 * * *\ / j * * *\ 

0 022 023 * I — ► I 0 1 * * I . Where /? = 

0 032 033 */ \0 0 /? *) 

032023/022 + 033. For the system to have a unique solution, we must have ( 3 ^ 0 . Simplify (3 to 
conclude 011022033 — 011023032 — 012021033 + 012023031 + 013021032 — 013022031 7^ 0. 
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CALCULATOR SOLUTIONS 1.3 

All solutions to CALCULATOR BOX problems will be given for the TI-85. Each will include a summary of the 
input keystrokes which were used to solve each type of problem, the output from the calculator, and the derivation of 
the answer from the output. Usually the fullest explanation will accompany the first example of each type of prob- 
lem. 

Most problems will require calculations to be made on some matrix or vector. Each such input for a problem should 
be entered into a variable before the calculations for that problem are begun. Once it is input, its’ value should be 
checked before the calculation for the problem is begun, since one of the most frequent causes of incorrect solutions 
is faulty input data. (We will not show the data input keystrokes in most cases, except in this first calculator solution 
section.) In order to allow the inputs (or outputs) for any problem to be reused or recalled in later problems, we shall 
tag the variable name(s) with numbers representing the chapter, section, and problem number; i.e. A1345 will be the 
name for the augmented matrix which is the input for chapter 1, section 3, problem 45. 

Each summary of input keystrokes follows the practice of the main text by boxing each input function keystroke 
(except for character or number keys which are displayed in the Courier font, and for which it is assumed that the 
appropriate (ALPHA] or (2nd] (alpha) keystrokes have been entered to allow alphabetic input to start or stop). 
When a keystroke sequence has selected a menu item, the named equivalent of that item will be displayed in 
Courier Bold inside angle brackets. For example the keystrokes to compute the reduced row echelon form of a 
matrix stored in the variable A1345 are displayed as: 

[2nd] [MATRX] (F4) <ops> (F5) <rref>A1345 (ENTER] . 

We will use the form MATRX ops rref to abbreviate later occurances of such a menu item entry. 

There is an alternative to all menu item function references. Since the name of any function, such as rref, is rec- 
ognized by the TI-85(even the all caps version rref), the result above can be produced by the (character) input 
rref A13 45 (ENTER) or even RREF A13 45 (ENTER) . We will ofter use this form of input in presenting solu- 
tions. 

44. To solve on the TI-85 enter the augmented matrix for the system by [[2. 6, -4. 3, 9. 6, 21 . 62] 
[-8.5,3.6,9.1,14.23] [12. 3, -8. 4, -.6, 12.61]] (STO ) A1344 (ENTER) 

Then compute the reduced row echelon form of A1344 by using the RREF command from the MATRX ops 
menu via (2nd) (MATRX) (F4) <ops> (F5) <RREF> | ALPHA) A1344 | ENTER) , 

to produce: 

[[100 86.1806588556 ] 

[010 122.285821022 ] 

[001 33.6853455595 ] ] 

which is the augmented matrix of the euqivalent (solved) system: 
xi = 86. 1806588556, x 2 = 122. 285821022, x 3 = 33. 6853455595. 

45. To solve, input the augmented matrix A1345 by [[0,2, -1,-4, 2], [1, -1,5,2, -4], [3, 3, -7, -1,4], 

[-1,-2, 3, 0,-7]] [STO» > ) A1345 (ENTER) , and verify that the matrix has been correctly entered by 
scrolling the display to examine all the entries of the input matrix. (Use the arrow keys: 0 , 0 , ar) 0 0 . 

0 if needed.) Note that the "missing" x t in the first equation is entered as a 0 in the first row. 

Now either use the MATRX ops menu, as described above, or literally enter rref A1345 (which requires 
the keystrokes | ALPHA) (ALPHA) RREF A (ALPHA) 1345 (ENTER) ) to get: 
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[[1000-3 ] 

[0 1 0 0 5 ] 

[0010 2e-14 ] 

[ 0 0 0 1 2 ] ] 

Since elementary row operations produce an augmented matrix of an equivalent system, the solutions to the 
original system can be read off from this equivalent (solved) system: 

xi=-3, x 2 = 5, x 3 = 2e- 14, x 4 = 2. 

(The calculator produced 2E-14, rather than the exact answer 0, due to rounding "errors" in its computation; it’s 
accuracy is limited (internally) to 13 significant figures, and many computations, like division by 3, may result 
in loss of accuracy due to roundoff.) 

46. Input the augmented matrix A1346: [[12.47,-2. 583, 7. 161, 8.275,-1. 205], 

[3. 472, 9. 283, 1 1. 275, 3. 606, 2. 374], [-5. 216, -12. 816, 6. 298, 1. 877, 21. 206], 

[6.812,5.223,-9.725,-2.306,-11.466]] [STOl A1346 and compute its reduced row echelon form R1346 
by entering RREF A1346 (STO ) R13 46 [ENTER | . Since we see the equations are consistent, the solu- 
tions are obtained from the last (5 th) column: (see problem 47) R1346(l,5,4,5) (ENTER) , 



[ [ 


2.22665461875 ] 


(= -tj) 


[ 


-1.93595628754 ] 


(= x 2 ) 


[ 


3.36239929557 ] 


(= x) 


[ 


-7.01511776944 ]] 


<"*4> 



47. Input the augmented matrix by [[23.42,-16.89,57. 31,82.6,2158.36], 

[-14. 77, 38. 29, 92. 36, -4. 36, -1 123. 02], [-77. 21, 71 . 26, -16. 55, 43. 09, 3248. 71], 

[91 . 82, 81 . 43, 33. 94, -57. 22, 235. 25]] I ENTER ) , and store it in A1347 by entering: 

(2nd) (ANS) I STO ) A1347 [ENTER ) . 

Then either follow the 1.44 or 1.45 solutions or produce and store the reduced echelon form in R1347 ("R" for 
reduced) via: 

RREF A1347 [STO ) R1347 

The reduced echelon form shows the equations are consistent, and the solution is the last (5th) column of 
R1347 which can be printed out by entering R1347(l,5,4,5) [ENTER ) , which yields the submatrix start- 
ing at (row=l,col=5) and ending at (row=4,col=5): 



[11.5606292935 


] 


(= *j) 


[27.8933709005 


] 


(= X,) 


[-19.8752502433 


] 


(= Xj) 


[42.3460010642 


]] 


(= x 4 ) 



48. From the input: 

[[6.1,-2.4,23.3,-16.4,-8.9,121.7] [-14. 2, -3 1.6, -5. 8, 9. 6, 23. 1.-87.7] 

[10. 5, 46. 1,-19. 6, -8. 8, -41.2, 10. 8] [37. 3, -14. 2, 62, 14.7,-9 6, 61. 3] 

[. 8, 17. 7, -47. 5, -50. 2, 29. 8, -27. 8]] [STO ) A1348 (ENTER) compute and store the reduced echelon form 
by RREF A1348 (STO ) 01348 (ENTER) . From the equivalent (consistent) system corresponding to this 
new augmented matrix, read off solution from the last (6th) column: 01348(1,6,5,6) (ENTER) 
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[ 


-4.19625216237 


] 


( = * ; ) 


[ 


3.39806581665 


] 


( = x 2 ) 


[ 


5.02950855373 


] 


( = x 3 ) 


[ 


-2.68699108976 


] 


(=V 


[ 


.651877193675 


] ] 


( = x 5 ) 



Problems 49-53 ask for the row echelon form of the augmented matrices from the equations in the previous 5 prob- 
lems. This will be computed by: 

(2nd) iMATRXl (F4) <ops> (F4) <ref> A [ENTER) 
or using only alphabetic and numeric input: 

REF A (ENTER) . 

49. The row echelon form for A1345 (saved in 1.45) is given by REF A1345 (ENTER) : 

[[11 -2.33333333333 -.333333333333 1.33333333333 ] 

[ 0 1 -.5 -2 1 ] 

[001 -.263157894737 -.526315789474 ] 

[000 1 2 ]] 



50. The requested row echelon form is given by REF A13 44 (ENTER) : 

[[ 1 -.682926829268 -.048780487805 1.02520325203 ] 

[ 0 1 -3.85314009662 -7.50853462158 ] 

[00 1 33.6853455595 ]] 

51. REF A13 47 (ENTER ] yields: 

[[ 1 .886843824875 .369636244827 -.623175778697 2.56207797865 ] 

[ 0 1 .085803613297 -.035964262264 24.6650599263 ] 

[00 1 -.125426892275 -25.1865775571 ] 

[00 0 1 42.3460010642 ]] 



52. REF A1346 (ENTER ) yields: 



[[ 1 -.20713712911 .574258219727 .663592622294 
[ 0 1 -.668756846388 -.384149034562 

[00 1 .322120802509 

[ 0 0 0 1 



-.096631916599 ] 
-1.48973311827 ] 
1.10268392998 ] 

-7.01511776944 ]] 



53. REF A1348 [ENTER ) yields: 
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[[ 1 -.380697050938 1.66219839142 .394101876676 -.257372654156 1.64343163539 ] 

[ 0 1 -.739622076066 -.258258724306 -.768456034635 -.128869813714 ] 

[00 1 1.29150456157 -1.23451620607 .754494344443 ] 

[00 0 1 -.245789833397 -2.84721587659 ] 

[00 0 0 1 .651877193675 ]] 

Problems 54-58 ask for all solutions, rounded to three decimal places, to certain systems with more unknowns 
than equations. To solve, first set the displayed precision to three decimal places by 

[2nd] [MODE] 0 0 0 0 0 [ENTER ) (EXTp . 

Then compute the reduced echelon form of the augmented matrix via rref A [ENTER) and write down the 
equivalent system for the resulting augmented matrix. If there is a new equation which says 0 = a non-zero num- 
ber, then this impossible equation shows the original system is inconsistent, i.e. has no solutions. Other wise, we 
can see that if we assign arbirary values to those variables which are not first in any of the resulting equations, 
then each equation can be solved for its first variable in terms of the arbitray variables by bringing the terms 
involving the arbitrary variables to the right side of the new equations (after changing the signs of those terms). 

This gives all possible solutions. 

54. The augmented matrix A1354 for the system is obtained by input of [[2. 1, 4. 2, -3. 5, 12. 9] 

[-5. 9, 2. 7, 9. 8, -1.6]] [STCH A13 54 [ENTER) . Then the reduced row echelon form RREF A13 54: 

[[10 -1.662 1.365 ] 

[ 0 1 -.002 2.389 ]] 

is the augmented matrix of a consistent system. In this equivalent system, we see that transposing the x 3 
terms to the right side yields the solutions (1. 365 + 1. 662x 3 , 2. 389+. 002x 3 ,x 3 ) with x 3 arbitrary. 

55. Input the augmented matrix A1355 for the system with [[-13.6, 71. 8,46. 3,-19. 5] [41. 3,-75,-82.9,46.4] 

[41. 8,65.4,-26.9,34.3]] [STO ) A1355 (ENTER) . Then the reduced row echelon form RREF A1355: 

[[10 -1.275 .961 ] 

[ 0 1 .403 -.009 ] 

[ 0 0 0 0 ] ] 

is the augmented matrix of the equivalent consistent system: 

X! -1.275x 3 = .961 

x 2 . 403x 3 = -.009' 

From this we see that transposing the x 3 terms to the right side yields the solutions (.961+1. 275x 3 , 

-.009-.403x 3 , x 3 ) with x 3 arbitrary. 

56. Since this system differs from the previous system only in two right hand side entries 

19.5 instead of -19.5 and 35.3 instead of 34.3 

the usual input can be omitted and the augmented matrix A1356 can be obtained by copying A1355 and edit- 
ing the 4’th column of the copy. (The keystrokes to copy and get to the editing stage are: A13 55 [STO 11 - ) 

A1356 [2nd] [MATRX ) (F2) <EDIT> .Then at the Name= prompt enter A13 56 [ENTER) (ENTER) 
(ENTER) . Now the 4’th column is displayed; use the arrow keys to position the cursor and make the 
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changes: delete the and change 3 4 to 3 5. Then enter [EXIT] and the new augmented matrix 

[[-13. 6, 71 . 8, 46. 3, 19. 5] [41 . 3, -75, -82. 9, 46. 4] [41 . 8, 65. 4, -26. 9, 35. 3]] is stored in A1356.) Then the 

reduced row echelon form RREF A1 3 5 6 : 

[[10 -1.275 0 ] 

[ 0 1 .403 0 ] 

[ 0 0 0 1 ] ] 

is the augmented matrix of an inconsistent system as the last row yields the equation 0=1. Thus there are no 
solutions. 



57. Input the augmented matrix A1357 for the system with: [[5,-2, 11,-16, 12, 105] [-6, 8, -14, -9, 26, -62] 
[7,-18,-12,21,-2,53]] [STO ) A1357 [ENTER] . Then the reduced row echelon form RREF A1357: 

[[100 -7.616 11.87 31.348 ] 

[010 -4.876 6.775 11.043 ] 

[001 1.121 -3.072 -2.696 ] ] 

shows x 4 and x 5 can be chosen arbitrarily. If we bring the terms involing these variables to the right side we 
get solutions ( 31.348+7.616x 4 -11.87x 5 , 11.043+4.876x 4 -6.775x 5 , -2.696-1. 121x 4 +3.072x 5 , x 4 , x 5 ). 



58. (To illustrate an alternative way to reuse data, we present a solution based on the prior entry of A1359. See 
problem 59 for a way to enter A1359.) Just copy A1359 into A1358 and change the [4,6] element by -63 
[STO] A13 58 (4,6). Then RREF A13 58 yields: 

[ [ 1 0 0 0 11.87 31.348 ] 

[0100 6.775 11.043 ] 

[0010 -3.072 -2.696 ] 

[ 0 0 0 1 0 0 ] ] 

from which we read off the solutions (31.348-11.87x s , 11.043-6.775x 5 , -2.696+3.072x 5 , 0, x 5 ) with x 5 arbi- 
trary. 

59. An alternative to the usual input of A1359 is to add a new row to A1357. This is done by A1 3 57 [STO ) 
A13 59 and then adding a new row (of zeros) by [2nd] (LIST ] <{>4,6<}> (STO) [2nd] [MATRX] 
copsxdinoAl 3 5 9 . Now use the MATRX<EDIT> function to change the fourth row elements to - 1 5 , 
42, 21, -17, 42, 63. Then the reduced row echelon form RREF A13 59 is 

[[ 1 0 0 0 11.87 50.54 ] 

[0100 6.775 23.33 ] 

[0010 -3.072 -5.52 ] 

[00010 2.52 ] ] 

from which we read off the solutions (50.54-1 1.87x 5 , 23.33-6.775x 5 , -5.52+3.072x s , 2.52, x 5 ) with x 5 arbi- 
trary. 
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MATLAB Tutorial 

The MATLAB input and output for these problems will be printed in this typewriter font. The 
symbol » is the MATLAB prompt. We will often suppress output by use of although you should al- 
ways check your input. 

1. They can be entered as 

» A = [22345; -6-1207; 12-134], b = [-1; 2; 5]; 
or as 

» A = [ 2 2345 

-6-1207 
12 - 134 ]; 

» b = [ -1 
2 

5 ]; 



2. The augmented matrix is: 



» C = [ A b] 

C = 

2 2 3 

-6 -1 2 

12-1 



4 5-1 
0 7 2 
3 4 5 



3. 



Notice that since this problem uses rand, you will get different numbers than those printed here. This 
generates a random 3x4 matrix with values between -1 and 1, and then multiplies that by 2. 



» D = 2*( 2*rand(3,4) - 1 ) 
D = 



-1.1242 0.7172 0.0777 

-1.8118 1.7388 1.3239 

0.7155 -0.4660 -1.8617 



-1.7862 
0.1188 
0 . 6846 



4. This generates a random 4x4 matrix with entries between -10 and 10, and then rounds off to the near- 
est integer. 

» B = roundf 10*( 2*rand(4,4) - 1) ) 

B = 

-10 4 1 4 

-2 2-8 8 
-9 9 3 5 

-2 7 -2 -5 




m Equations in n unknowns: Gauss-Jordan and Gaussian Elimination MATLAB Tutorial 1. 



19 



5. This first copies B into K , and then reverses the two rows in K . 



» K = B; 
K = 


K([ 1 4], 


:) = K( [4 1],:) 


-2 


7 -2 


-5 


-2 


2 -8 


8 


-9 


9 3 


5 


-10 


4 1 


4 



» C(3, : ) = 
c = 


C(3 , : ) + 


(~1/2)*C(1, 


:) 






u — 

2.0000 


2.0000 


3.0000 


4.0000 


5.0000 


-1.0000 


-6.0000 


-1.0000 


2.0000 


0 


7.0000 


2.0000 


0 


1.0000 


-2.5000 


1.0000 


1.5000 


5.5000 



» B( [2 


4] , [1 3]) 


'/, This is the 2x2 submatrix of B 


» 




V, made from the second and fourth rows 


» 




*/. and the first and third columns. 


ans = 






-2 


-8 




-2 


-2 





>. Recall that D from problem 3 was a random matrix, so your values will be different. 



» U = D( : , [3 4]) 

U = 

0.0777 -1.7862 

1.3239 0.1188 

-1.8617 0.6846 



» C(2, : ) = C(2, : ) + 3*C(1,:) 



2.0000 


2.0000 


3.0000 


4.0000 


5.0000 


- 1.0000 


0 


5.0000 


11.0000 


12.0000 


22.0000 


-1.0000 


0 


1.0000 


-2.5000 


1.0000 


1.5000 


5.5000 
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10. This will generate a random matrix: 
» T = rand(8 ,7) 



Notice that if you type 

» help : '/, Use this in raatlab version 3.5 



or 



» help colon '/, Use this in matlab version 4.0 
that 3 : 8 is the same as [34567 8], so 



» S = T( 3:8 ,:) 

will generate rows 3 through 8 of T. 

11. The reduced row echelon form of C will be: 



» rref(C) 
ans = 

1.0000 


0 


0 


-0.1915 


-1.4681 


-1 . 1489 


0 


1 . 0000 


0 


1.7447 


3.0426 


2.4681 


0 


0 


1.0000 


0.2979 


0.6170 


-1.2128 



So that an equivalent system of equations would be: 



xi -0.1915x4 -1.4681x5 = -1.1489 
x 2 +1.7447x4 +3.0426x 5 = 2.4681 

x 3 +0.2979x4 +0.6170x5 = -1.2128 
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MATLAB 1.3 

1. For each problem, first A is entered as the augmented matrix representing the system of equations. 

Next R is set to be the reduced row-echelon form of A. R represents a system whose solution, x, is 
just the last column of R. Since in each case, the system reduces to one where no variables may be 
chosen arbitrarily, there is a unique solution. 

For problem 1: 

» A = [ 1 -2 3 11; 4 1 -1 4; 2 -1 3 10]; 

» R = rref(A) 

R = 

10 0 2 

0 10-3 

0 0 11 

» x = R(:,4) '/, Equivalent system says jth variable 

x = '/. Equal to the jth entry in the last column. 

2 

-3 

1 

For problem 2: 

» A = [-2 1 6 18; 5 0 8 -16; 3 2 -10 -3]; 

» R = rref(A) 

R = 

1.0000 0 0 -4.0000 

0 1.0000 0 7.0000 

0 0 1.0000 0.5000 

» x = R(: ,4) 
x = 

-4.0000 

7.0000 
0.5000 

For problem 5: 

» A = [1 1 -1 7; 4-154; 22-30]; 

» R = rrei(A) 

R = 

10 0-9 

0 1 0 30 

0 0 1 14 

» x = R(: ,4) 
x = 

-9 

30 

14 
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For problem 8: 



» A = [1 -2 3 0; 4 1 -1 0; 2-130]; 
» R = rref(A) 

R = 

10 0 0 
0 10 0 

0 0 10 

» x = R(: ,4) 
x = 

0 

0 

0 



For problem 16: 



» A = [1-2112; 302-2-8; 



»R = rref(A) 
R = 




1.0000 0 


0 


0 1.0000 


0 


0 0 


1.0000 


0 0 


0 


» x = R(: ,5) 





04-1-11; -1 6 -2 0 7] ; 



0 2.0000 

0 0.5000 

0 -3.0000 

1.0000 4.0000 



2.0000 

0.5000 

-3.0000 

4.0000 



2. For each problem, first A is entered as the augmented matrix representing the system of equations. 
Next R is set to be the reduced row-echelon form of A. Since the bottom row of R represents an 
equation 0=1, there can be no solutions to this system. 

For problem 4: 



» A = [3 6 -6 9; 2 -5 4 6; 5 28 -26 -8] ; 
» R = rref(A) 

R = 



1.0000 0 - 0.2222 

0 1.0000 -0.8889 

0 0 0 



0 

0 

1.0000 



For problem 7: 

» A = [1 1 -1 7; 4 -1 5 4; 6 1 3 20] ; 

» R = rref(A) 

R = 

1.0000 0 0.8000 0 

0 1.0000 -1.8000 0 

000 1.0000 
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For problem 13: 

» A = [12-44; -2 -4 8 -9] ; 

» R = rref(A) 

R = 

12-40 
0 0 0 1 

For problem 18: 

» A = [1 -2 1 1 2; 3 0 2 -2 -8; 0 4 -1 -1 1; 5 0 3 -1 0] ; 

» R = rref(A) 

R = 

1 0 0 4 0 

010-20 
001-70 
0 0 0 0 1 



3. (i) (a) For the matrix we have 

» A = [ 3 5 1 0; 4 2 -8 0; 8 3 -18 0] ; 

» R = rref(A) 

R = 

10-30 
0 12 0 
0 0 0 0 

(b) The pivots have been underlined, (c) An equivalent system of equations would be 

- 3®3 = 0 
X2 +2x 3 = o 

(d) No pivot in column 3, so the solution of this system has: x 3 arbitrary, xi = 3x3, X2 = —2x3. 
(ii) (a) 

» A = [ 9 27 3 3 12; 9 27 10 1 19; 13596]; 

» R = rref(A) 

R = 

1 3 0 0 1 

0 0 10 1 
0 0 0 1 0 



(c) An equivalent system of equations would be 

Xi +3X2 = 1 

x 3 = 1 
x 4 = 1 
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(d) No pivot in column 2 so the solution of this system is: *2 arbitrary, *i = 1 — 3*2, *3 = 1 , 

*4 = 0 . 

(™) ( a ) 



» A = [101-27-4; 1 4 21 -2 2 5; 3 0 3 -6 7 2] ; 
» R = rref(A) 

R = 

10 1-203 



0 15 0 

0 0 0 0 



0 1 
1 -1 



(c) An equivalent system of equations would be 



*1 +*3 —2*4 = 3 

*2 +5*3 = 1 

*5 = -1 

(d) No pivots in columns 3 and 4 so the solution of this system is: *3 and *4 arbitrary, *1 = 
3 - *3 + 2*4, *2 = 1 — 5*3, *5 = — 1 . 

(iv) (a) 



» A = [6 


4 7 5 15 


9 








8 


59 10 10 


8 








4 


5 7 7 -1 


7 








8 


3 7 6 22 


8 








3 


2 7/2 9 -12 


-2]; 








» R = rref(A) 










R = 












1.0000 


0 


0.5000 


0 


5.0000 


1.0000 


0 


1.0000 


1.0000 


0 


0 


2.0000 


0 


0 


0 


1.0000 


-3.0000 


- 1.0000 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


(c) An equivalent system of equations would be 







*1 +. 5*3 +5*5 = 1 

*2 +1*3 = 2 

*4 —3*5 = — 1 



(d) No pivots in columns 3 and 5 so the solution of this system is: *3 and *5 arbitrary, *1 = 
1 - .5*3 - 5*5, *2 = 2 - * 3 , and *4 = -1 + 3*5. 
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4. (i) Reduce the augmented matrix representing the equations: 
» rref( [ 1 2 3 -1; 0 -3 1 4; 4 1 -2 0] ) 



1.0000 


0 


0 


0.4694 


0 


1.0000 


0 


-1.2245 


0 


0 


1 . 0000 


0.3265 



From this, the solution is x\ = 0.4694, X 2 — —1.2245 and X 3 = 0.3265. Since there is only one 
solution, these three planes intersect in exactly one point. 

(ii) Reduce the matrix as in (i): 

» rref( [2-145; 1 2 -3 6; 4 3 -2 9] ) 
ans = 

10 10 
0 1-20 
0 0 0 1 



There are no solutions, i.e. the system is inconsistent. This means the three planes do not inter- 
sect. 



(iii) 



» rref ( [ 2 


-1 4 5; 


12-36; 


43-2 17] 


ans = 

1.0000 


0 


1.0000 


3.2000 


0 


1.0000 


-2.0000 


1.4000 


0 


0 


0 


0 



The solution has X 3 arbitrary, X 2 = 1.4 + 2 x 3 , and x\ = 3.2 — X 3 . The planes intersect in a line, 

(iv) 



» rref ( [2-424; 3-636; -12-1 -2] ) 
ans = 

1-212 
0 0 0 0 

0 0 0 0 

The solution has X 3 and 12 arbitrary, and x\ = 2 + 2*2 — x z- The three planes are identical. 

5. (i) 



» A = [1 
» D=A; 

» A(2, : ) 

» A(3, : ) 

A = 

12-12 
0 0 4 4 

0 -2 -4 -6 

» A( [2 3],:) = A([3 2], :) 
A = 

12-12 
0 -2 -4 -6 

0 0 4 4 



7, Subtract 2*R1 from R2. 
'/, Subtract 3*R1 from R3. 



'/, Interchange R2 and R3. 



2 - 12 ; 2428 ; 34 - 70 ]; 

= A(2, : ) - 2*A(1 , : ) ; 

= A(3, : ) - 3*A(1 , : ) 




26 Chapter 1 Systems of Linear Equations and Matrices 



Instructor's Manual 



» A(2, : ) = A(2, : ) / (-2); */. Normalize R2. 

» A ( 1 , : ) = A ( 1 , : ) - 2*A(2,:) */. Subtract 2*R2 from Rl. 

A = 

1 0-5-4 

0 12 3 

0 0 4 4 



» A(3 , : ) = A(3, : ) / 4; '/. Normalize R3 

» A(2, : ) = A(2 , : ) - 2*A(3,:); '/, Subtract 2*R3 from R2. 

» A(l, : ) = A(1 , : ) + 5*A(3,:) '/. Subtract -5*R3 from Rl . 

A = 

10 0 1 
0 10 1 

0 0 11 

Compare this with: 

» rref(D) 
ans = 

10 0 1 
0 10 1 

0 0 11 

(ii) 

» A = Cl 2 3 2; 34-1-3; -2104]; 



» D = A; 


» A(2, : ) = A(2, : ) - 


3*A(1 , : ) ; 


*/, Subtract 3*R1 from R2. 


» A(3, : ) = A(3, : ) + 


2*A(1 , : ) 


•/. Subtract -2*R1 from R3 


A = 


12 3 


2 




0 -2 -10 


-9 




0 5 6 


8 




» A(2, : ) = A(2,:> / 


(-2); 


•/. Normalize R2. 


» A(l, : ) = A(1 , : ) - 


2*A(2, : ) ; 


*/, Subtract 2*R2 from Rl . 


» A(3 , : ) = A(3, : ) - 


5*A(2, : ) 


'/, Subtract 5*R2 from R3. 


A = 


1.0000 0 


-7.0000 


-7.0000 


0 1.0000 


5.0000 


4.5000 


0 0 


-19.0000 


-14.5000 


» A(3 , : ) = A(3, : ) / 


(-19); 


*/, Normalize R3. 


» A(2, : ) = A(2, : ) - 


5*A(3, : ) ; 


'/. Subtract 5*R3 from R2. 


» A(1 , : ) = A(1 , : ) + 


7*A(3, : ) 


'/. Subtract -7*R3 from Rl 


A = 


1.0000 0 


0 


-1.6579 


0 1.0000 


0 


0.6842 


0 0 


1.0000 


0.7632 
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Compare this with: 

» rref (D) 
arts = 

1.0000 
0 
0 

(iii) 

» A = [1 2 -2 0 1-2 

24-10-4 -19 
-3 -6 12 2 -12 -8 

1 2-2-4 -5 -34]; 

» D = A; 

» A(2,:) = A(2,:> - 2*A(1,:); ’/, Subtract 2*R1 from R2. 

» A(3,:) = A(3,:) + 3*A(1,:); '/, Subtract -3*R1 from R3. 

» A(4, : ) = A(4, : ) - 1*A(1,:) '/. Subtract R1 from R4. 

A = 

1 2-2 0 1-2 
0 0 3 0 -6 -15 

0 0 6 2 -9 -14 

0 0 0 -4 -6 -32 

» A(2, : ) = A(2, : ) / (3); */, Normalize R2. 

» A(l,:) = A(l,:) + 2*A(2,:); '/. Subtract -2*R2 from R1 . 

» A(3,:) = A(3,:) - 6*A(2,:) '/, Subtract 6*R2 from R3. 

A = 

1 2 0 0 -3 -12 

0 0 1 0-2-5 

0 0 0 2 3 16 

0 0 0 -4 -6 -32 

» A(3 , : ) = A(3 ,:) / (2); */, Normalize R3. 

» A(4, : ) = A(4, : ) + 4*A(3,:) */. Subtract -4*R3 from R4. 

A = 



1.0000 


2.0000 


0 


0 


-3.0000 


-12.0000 


0 


0 


1.0000 


0 


-2.0000 


-5.0000 


0 


0 


0 


1.0000 


1.5000 


8.0000 


0 


0 


0 


0 


0 


0 



Compare this with: 



ref (D) 












1.0000 


2.0000 


0 


0 


-3.0000 


-12.0000 


0 


0 


1.0000 


0 


-2.0000 


-5.0000 


0 


0 


0 


1.0000 


1.5000 


8.0000 


0 


0 


0 


0 


0 


0 



0 0 -1.6579 

1 . 0000 0 0 . 6842 

0 1.0000 0.7632 
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6. (a) First enter A and b and then let C be the augmented matrix. 

» A =[ 1 2 -2 0; 2 4 -1 0; -3 -6 12 2; 12 -2 -4] ; 

» b = [1; -4; -12; 3] ; 

» C = [A b] 

C = 

1 2-2 0 1 

2 4-1 0-4 

-3 -6 12 2 -12 

1 2-2-4 3 

» rref (C) '/, Reduce the augmented matrix, 

ans = 

1 2 0 0 0 

0 0 10 0 

0 0 0 1 0 

0 0 0 0 1 

Since the bottom row represents the equation Oxi + 0 x 2 + OX 3 + OX 4 = 1, which has no solution, 
the system has no solutions. 

(b) 

» b = 2*A( : , 1) + A( : ,2) + 3*A(:,3) - 4*A(:,4) 
b = 

-2 

5 

16 

14 

» rref ( [A b] ) 
ans = 

1 2 0 0 4 

0 0 10 3 

0 0 0 1 -4 

0 0 0 0 0 

The solution is x 4 = -4, X 3 = 3, x 2 is arbitrary, and xi = 4 - 2x 2 . 

(c) For any choice of coefficients, there will be a solution. 

(d) No, it is not possible. The conjecture is true. For any choice of a; if 

b = aiA(:, 1) + a 2 .<4(:, 2) + asA(:, 3) + a 4 yl(:, 4), 

then Ax = b will have at least the one solution x,- = a,-. 

(e) Repeating the experiment will generate a random singular matrix A. In each case, if b is a sum 
of multiples of columns of A, the system [A b] will have a solution. 

7. (a) 

» A = [ 1 1 1 ; 234; -20 3 3; 

» b= [4; 9; -7 ] ; 

» c = [4; 16; 11 ]; 

» Aug = [Abe] 

Aug = 

1114 4 

2 3 4 9 16 

-2 0 3 -7 11 

» rref (Aug) 
ans = 

1 0 0 2 -1 

0 10 3 2 

001-13 
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(b) 



The solution to the first system is x = (2, 3, —1) and to the second is x = (—1, 2, 3). 



» A = [ 2 3 -4; 1 2 -3; -1 5 -11]; 

» b = [ 1; 0; -7]; 

» c = [- 1 ; - 1 ; - 63 ; 

» d = [1; 2; -7]; 

» Aug = [Abed] 

Aug = 

2 3-4 1-1 1 

1 2-3 0-1 2 

-1 5 -11 -7 -6 -7 

» rref(Aug) 
ans = 

10 12 10 
0 1 -2 -1 -1 0 

0 0 0 0 0 1 



For the first system, X 3 may be chosen arbitrarily, xi = 2 — X3, ® 2 = — 1 + 2^3. For the second system, 
*3 may be chosen arbitrarily, x\ — 1 — £ 3 , x 2 = — 1 + 2 x 3 . The third system is inconsistent, and has 
no solutions. 

(c) Any columns will generate solutions using the same method as above. 

(d) (i) No. Since the number of variables that may be chosen arbitrarily is not determined by the 

right hand side, a system cannot have a unique solution with one right hand side and in- 
finitely many solutions with another. 

(ii) No. A system will have no solutions only when it reduces to a system with zeros on the left 
hand side, and nonzeros on the right hand side, so one of the columns on the left will not 
have a pivot. However, if a system will have a unique solution if every column on the left has 
a pivot. Both cases cannot happen for the same left hand side. 

(iii) Yes. In (b) above, there were infinitely many solutions for the first two systems, but no solu- 
tion for the third. This happened because there was a column on the left hand side without 
a pivot. This missing pivot may or may not appear on the right hand side, causing the sys- 
tem to be inconsistent or consistent. 

8 . (a) If we consider each node in numerical order we find the equations: 



7i =(100 + T 2 + T 4 + 50)/4, or 
T 2 =(100 + T s + T 8 +Ti)/4, or 
T 3 =(100 + 50 + Tg + T 2 )/4, or 
T 4 = (Ti + T 5 + TV + 50)/4, or 
T 5 = (T 2 +T 6 + T 8 + T 4 )/4, or 
Tg = (T3 + 50 + T9 + T5)/ 4, or 
T 7 = (T 4 + T 8 + 0 + 50)/ 4, or 

T 8 = (T 4 + X9 + 0 + TV)/ 4, or 

T9 = (Tg + 50 + 0 + T 8 )/4, or 



4T\ - T 2 - T 4 =150 
— TV + 4T 2 -T3-T3 =100 
-T 2 + 4T 3 - Tg =150 
-T\ + 4T 4 - T 5 - T> = 50 
-TV +4T 5 -Tg -T 8 = 0. 
-TV - TV + 4Tg -TV = 50 
-T 4 + 4TV - T s = 50 
-TV - TV + 4TV - TV = 0 
-Tg - Tg + 4T 9 = 0 
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To express the equations on the right as AT = bwe see that we can form the coefficient matrix 
A and the righthand side b as follows: 

» A = 4*eye(9); '/, The diagonal terms axe all 4 and the non-zero 

» A ( 1 , [2 4]) = -Cl 1]; A(2, [1 3 5]) = -[1 1 1] ; */. off diagonals are -1 

» A(3, [2 6]) = -Cl 1]; A(4,Cl 5 7]) = -Cl 1 1]; 

» A (5, C2 6 8]) = -Cl 1 1]; A(6,C3 5 9]) = -Cl 1 1] ; 

» A(7,C4 8]) = -Cl 1]; A(8,C5 7 9]) = -Cl 1 1]; 

» A(9,C6 8]) = -Cl 1]; 

» b=C 150; 100; 150; 50; 0; 50; 50; 0; 50]; 

» CA b] '/. Here is the augmented matrix for the system, 
ans = 



4 


-1 


0 


-1 


0 


0 


0 


0 


0 


150 


-1 


4 


-1 


0 


-1 


0 


0 


0 


0 


100 


0 


-1 


4 


0 


0 


-1 


0 


0 


0 


150 


-1 


0 


0 


4 


-1 


0 


-1 


0 


0 


50 


0 


-1 


0 


0 


4 


-1 


0 


-1 


0 


0 


0 


0 


-1 


0 


-1 


4 


0 


0 


-1 


50 


0 


0 


0 


-1 


0 


0 


4 


-1 


0 


50 


0 


0 


0 


0 


-1 


0 


-1 


4 


-1 


0 


0 


0 


0 


0 


0 


-1 


0 


-1 


4 


50 



Notice that the non-zero terms are (relatively) near the diagonal; specifically there is a diagonal 
band about the main diagonal containing all non-zero entries. 

» R = rref ( CA b] ) ; 

» R( : , 10) ’ 
ans = 

Columns 1 through 7 
65.0794 65.8730 65.0794 

Columns 8 through 9 
23.0159 29.3651 

i.e. Tj = 65.0794, T 2 = 65.8730, T 3 = 65.0794, T 4 = 44.4444, T 5 = 33.3333, T 6 = 44.4444, T r = 
29.3651, T s = 23.0159, T 9 = 29.3651. 



» 

» y = (A\b) ; y’ 
ans = 

Columns 1 through 7 

65.0794 65.8730 65.0794 44.4444 33.3333 44.4444 29.3651 

Columns 8 through 9 
23.0159 29.3651 



V, The initial 9x9 in rref(CA b] ) is I so 
'/, The solution is just the 10 ’th column 

44.4444 33.3333 44.4444 29.3651 
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9. (a) First set A to be the matrix, and b to be the right hand side: 



» A = [ 1-.2 -.5 -.15; -.4 l-.l -.3; -.25 -.5 1-.15; ] 

A = 

0.8000 -0.5000 -0.1500 

-0.4000 0.9000 -0.3000 

-0.2500 -0.5000 0.8500 

» b = [10; 25; 20] 

b = 

10 

25 

20 

» A\b 



ans = 

110.3058 

118.7429 

125.8211 



(b) (i) The value 032 = .05 tells us that industry 2 needs .05 units of output from industry 3 in order 
to manufacture one unit. The value <133 = 0 tells us that industry 3 needs none of its own 
output. 

(ii) The augmented matrix will be: 



» A = [ 1-.2 
-.15 
-.1 

A = 

1.0e+05 * 
0.0000 0 
0.0000 0 
0 . 0000 0 



-.1 -.3 

1-.25 -.25 

-.05 1-0 



0000 0 . 0000 

0000 0.0000 

0000 0.0000 



300000 

200000 

200000] 



3.0000 

2.0000 
2.0000 



MATLAB has printed only the most significant digits, .8 is very small compared to 300,000 
so it is rounded off to 0. It does, however, keep the smaller numbers in memory: 



» A( : , 1:3) 
ans = 

0.8000 -0.1000 -0.3000 
-0.1500 0.7500 -0.2500 
-0.1000 -0.0500 1.0000 



Also, A is printed using “scientific notation”. The “1.0e+05 *” tells us that we must multi- 
ply every number in A by 100,000. 
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(iii) 



» R = rref(A) */. This reduces A to row echelon form. 
R = 



1.0e+05 * 

0.0000 0 

0 0.0000 

0 0 



0 5.3720 

0 4.6645 

0.0000 2.7704 



» R(:,l:3) */, This is the reduced echelon form of the coefficient part: 



ans = 

1 0 

0 1 

0 0 

» x = R( : ,4) 



0 

0 

1 

'/, The solution vector. 



x = 

1.0e+05 * 
5 . 3720 
4.6645 
2 . 7704 



In order to balance supply and demand, industry 1 should make 537,200 units, industry 2 
should make 466,450 units, and industry 3 should make 277,040 units, to 5 significant digits. 



» format long 
» x 

X = 

1.0e+05 * 
5.37197626654496 
4.66453674121406 
2.77042446371520 
» format 



10. (a) The equation for each intersection will be: 



(The negation of each is also a valid equation). 



(b) 



at [1] xi —x 3 


+x 5 = 200 


at [2] — x\ +x 2 


- 0 


at [3] -x 2 +Z 3 


O 

O 

i“H 

1 

II 

1 


at [4] 


x 4 —x 5 = —100 



» A = [1 0-1 0 1 200 

-1 1 0 0 0 0 
0-1 1-1 0 -100 
000 1-1 -100 ] ; 

» rref(A) 
ans = 

1 0-1 0 1 200 

0 1-1 0 1 200 

0 0 0 1 -1 -100 

0 0 0 0 0 0 



We may choose x 3 and x 5 arbitrarily, then xi = x 3 — x 5 +200, x 2 = x 3 — 2 : 5 +200 and x 4 = x 5 — 100. 
(c) If we set x 5 = 0, then x 4 = —100, i.e. the traffic from [3] to [4] would have to be reversed. The 
smallest x$ can be chosen is 100 , in order to keep all of the other numbers nonnegative. 
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11. (a) As in the example, the equations to solve will be Ax = b, with: 

» A = [ 1‘2 1 1; 3"2 3 1; 4~2 4 1] 

A = 

1 1 1 

9 3 1 

16 4 1 

» b = [-1; 3; -2]; 

» x = A\b 
x = 

-2.3333 

11.3333 

- 10.0000 



The parabola will be y = —2.3333a: 2 + 11.3333* — 10. 

» x = [1; 3; 4]; 

» V = vander(x) 

V = 

1 1 1 

9 3 1 

16 4 1 



V is the same matrix as A. 

(b) This is similar to (a), except we now need a third degree polynomial to fit four points: 



» A = C 0'3 0'2 0 1; 

A = 

0 0 0 

111 
27 9 3 

64 16 4 

» b = [5; -2; 3; -2]; 
» x = A\b 



1*3 1'2 1 1; 3“3 3-2 3 1; 4*3 4~2 4 1] 

1 

1 

1 

1 



-1.4167 

8.8333 

-14.4167 

5.0000 



The cubic polynomial will be y = — 1.4167a; 3 + 8.8333a: 2 — 14.4167a: + 5. 



» x = [0; 1; 3; 4]; 
» V = vander(x) 

V = 



0 0 0 

1 1 1 

27 9 3 

64 16 4 



1 

1 

1 

1 



V is the same matrix as A. 
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(c) For part (a): 



» x = Cl; 3; 4]; 

» y = [-1; 3; -2]; 

» V = vander(x); 

» c— V\y 
c = 

-2.3333 

11.3333 

- 10.0000 

» s= ntin(x) : .01:max(x); 
» yy = polyval(c,s) ; 

» plot(x,y, , s,yy) 



For part (b): 



» x = [0; 1; 3; 4] ; 

» y = [5; -2; 3; -2] 

» V = vander(x); 

» c= V\y 
c = 

-1.4167 
8.8333 
-14.4167 
5.0000 

» s= min(x) : .01:max(x); 

» yy = polyval(c.s) ; 

» plot(x,y, ,s,yy) 

(d) This will generate a 7th degree polynomial, passing through each of the seven points. 

» x = rand(7,l) 
x = 

0.0668 

0.4175 

0.6868 

0.5890 

0.9304 

0.8462 

0.5269 

» y = rand (7,1) 

y = 

0.0920 

0.6539 

0.4160 

0.7012 

0.9103 

0.7622 

0.2625 
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» V = vander(x); 

» c= V\y 
c = 

1.0e+04 * 

-0.9529 

3.3259 

-4.6043 

3.1855 

-1.1279 

0.1808 

-0.0079 

» s= min(x) : .01:max(x) ; 
» yy = polyval(c,s) ; 

» plot(x,y, ,s,yy) 




Notice that even though the y coordinates of the original points were all between 0 and 1, the result- 
ing polynomial can oscillate dramatically. 
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Note: Any variables appearing in a solution can take arbitrary values. 

L (3 1 0) (0 11/2 0) (J 1 2 )' Solution: 



Solution: (0,0). 



1 -5 0 
-15 0 



o) _+ (o "o 0 )' SoIution: ( 5 x 2 ,* 2 )- 



11-10 
2-430 
3 7-10 



11-10 
0-650 
0 4 2 0 



10 1/60 
0 1 -5/6 0 
0 0 16/3 0 



1 0 0 0\ 

0 10 0 I . Solution: (0, 0, 0). 
0 0 10 / 



11-10 
2-430 
-1-7 6 0 



11-10' 
0-650 
0 -6 5 0, 



'1 0 -1/6 O' 
0 1 -5/6 0 
0 0 0 0 , 



Solution: ( 2 : 3 / 6 , 5 ^ 3 / 6 , * 3 ). 



/ 1 1-10 
5. 2 -4 3 0 

\ -5 13 -10 0, 



'11-10 
0-6 5 0 

0 18 -15 0 



1 0 - 1/6 0 \ 

0 1 -5/6 0 1. Solution: ( 2 : 3 / 6 , 5 * 3 / 6 , 2 : 3 ). 
0 0 0 0 / 



2 3-10 
6-570 



1 3/2 -1/2 O' 
0 -14 10 0 



(0 1 - 5/7 o)' Solution: (- 4 * 3 / 7 , 52 : 3 / 7 , 2 : 3 ). 



4-10 
7 3 0 
-8 6 0 



1 -1/4 0 
0 19/4 0 
0 4 0 



1 0 0\ 

0 1 0 1. Solution: (0,0). 
0 0 0 / 



fl -1 7 -1 0\ fl -1 7 -1 O') 

V2 3 -8 1 Oj \ 2 5 -22 3 0 / 

Solution: ((— 13*3 + 2*.j)/5, (22*3 — 3*4)/5, *3, *4). 



1 0 13/5 -2/5 O' 
0 1 -22/5 3/5 0, 



1-2110' 
302-20 
0 4-1-1 0 
503-10/ 



'1-2 1 1 O' 
0 6-1-5 0 

0 4-1-1 0 

.0 10 -2 -6 0 > 



Solution : ( — 4*4 , 2*4 , 7*4 , *4) . 



'10 2/3 -2/3 O' 

0 1 -1/6 -5/6 0 

0 0 -1/3 7/3 0 

,0 0 -1/3 7/3 0/ 



1 0 0 4 O' 
010-2 0 
001-7 0 
0 0 0 0 0 / 



-2 0 0 7 0\ 
12-14 0 
3 0 -1 5 0 
42 30 0/ 

/I 0 0 -7/2 O' 

010 -40 

0 0 1 -31/2 0 

\0 0 0 155/2 0, 



'100 -7/2 On 

0 -2 1 -15/2 0 
00-1 31/2 0 
,0 2 3 28/2 0y 
/I 0 0 0 O' 
[ 0 1 0 0 0 
~ + 0 0 10 0 
Vo 0 0 1 0/ 



( 1 0 0 -7/2 O' 

0 -2 -1/2 15/4 0 
00-1 31/2 0 
00 4 13/2 0/ 

. Solution (0, 0, 0, 0). 
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/ 2 -1 0\ /I -1/2 0\ 

„ | 3 5 0 1 | 0 13/2 0 1 

1L 7 -3 0 P 0 1/2 0 

V -2 3 0/ \0 2 0/ 



/I -3 0\ /I -3 0\ 

12. I -2 6 0 1 — + I 0 0 0 J . Solution: (3x2) * 2 )- 

\ 4 -12 0/ \0 0 0/ 



/ 1 1 -1 0 \ /I 1-1 0 \ /I 1 -1 0 \ 

„ 4 -1 5 0 [ 0 -5 9 0 [0-1-30 

-2 1-2 0 "" 0 3 -4 0 03-40 

V 3 2 -6 0/ Vo -1 -3 0/ Vo -5 9 0/ 

Solution: (0,0,0). 

14. If an = a2i = 0, then xi is arbitrary and therefore infinitely many solutions and 011022 — 012021 = 0. 
If either an or a 2 i is non-zero, (say an ^ 0), then 

/ an an 0 \ ^ / 1 a 12 /a n 0\ 

V°21 a 22 0/ Vo ( a ll a 22 — a 12«2l)/ail 0/ 

There will be an infinite number of solutions when an ^ 0 if and only if (011022 — ai2 a 2i/)/an = 0. 
This is true if and only if ana 2 2 ~ «i2«2i = 0. Similarly o 2 2 / 0, get infinite solutions if and only if 
On022 — U12O21 = 0. 

/ 2 -3 5 0\ /I -3/2 5/2 0\ /l 0 32/11 0\ 

15. I -1 7-1 0 I — »- J 0 11/2 3/2 0 J — *■ I 0 1 3/11 0 ]. In order to have a non-trivial 

\4 -11 k 0/ V° - 5 k ' 10 0/ \° 0 k-95/11 0/ 

solution, we need k — 95/11 = 0. Therefore, k = 95/11. 

16. Repeat the solution to Problem 43 in Section 1.3 to see that for a unique solution we need 011022033 — 
011023032 — 012021033 + 012023031 + 013021032 — 013022031 ^ 0. 



( 10 -4 0\ 
0 1 3 0 

0 0 -13 0 
0 0 24 0/ 



/I 0 0 0\ 

0 10 0 
0 0 10 ' 
Vo 0 0 0/ 
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CALCULATOR SOLUTIONS 1.4 

Refer to the CALCULATOR SOLUTIONS 1.3 to review the conventions followed in presenting TI-85 solutions. If 
the input for one of these homogeneous systems can be obtained by simple modifications to the (saved) input for a 
previous problem, the solution will outline how to do that. Sometimes the appropriate augmented matrix can be 
formed by extracting the coefficient submatrix (of a previous augmented matrix) and augmenting with a new right 
hand side, consisting of all zeros. At other times an entire new row will be added to an existing (augmented) matrix 
by first copying the previous matrix to a new variable and then editing the new matrix using the [2nd] [MATRX 1 
<EDIT> menu entry, which allows changing dimensions and specifying specific elements of the matrix being edited. 
Also recall that the TI-85 recognizes both rref and RREF as a name for the HATRX ops rref reduced row 
echelon form function. (It is slightly easier to key in the all upper case version of most function names.) 

17. This is the homogeneous system associated to Section 1.3, Problem 54. To solve it we form A1417 by (2nd] 

(MATRX 1 <ops> (MOREl (fT) <aug>A13 54 ( 1, 1, 2 , 3 ) ,[0,0]) [STO ] A1417. Then from the equiv- 
alent system derived from RREF A1 4 1 7 : 

[[10 -1.66206896552 0 ] 

[ 0 1 -.002298850575 0 ] ] 

we find that the solutions are (1.66206896552x 3 , 2.29885057472e-3x 3 , x 3 ) with x 3 arbitrary. 

18. This is the homogeneous system associated to Problem 55 in Section 1.3. We use 

RREF (AUG (A1 3 55 (1,1, 3,4), [0,0,0])) (ENTER) to produce the reduced echelon form: 

[[10 -1.27469748219 0 ] 

[ 0 1 .403399919808 0 ] 

[ 0 0 0 0 ] ] 

from which we find the solutions are (1.27469748219x 3 , ,403399919808x 3 , x 3 ) with x 3 arbitrary. 

19. We input the augmented matrix A1419 by [[25, -16, 13, 33, -57, 0] [-16, 3, 1, 0, 12, 0] [0, -8, 0, 16, -26, 0]] 

[STO 1 A14 19, being careful to get the zero’s in the correct places to represent the missing variables. Then 
from RREF A1419: 

[[100 -.330472103004 -.146995708155 0 ] 

[010-2 3.25 0 ] 

[001 .712446351931 -.101931330472 0 ]] 

we find both x 4 and x 5 arbitrary and the solutions are ( .330472103004x 4 +.146995708155x 5 , 2x 4 -3.25x 5 , 
-,712446351931x 4 +. 101931330472x5, x 4 , x 5 ). 

20. The first three equations have the same coefficients as Section 53, Problem 57, so A1420 can be formed by 
A1357 (STO* 1 A1 4 20, and using (MATRX | <EDIT>A1420 (ENTER] to edit this new augmented matrix. 
First we change the number of rows to 4 and leave the number of columns at 6 by 4 (ENTER) (ENTER) and 
then we use the arrow keys and <col>> to edit the bottom row to contain -1, 1 1, -9, 13, -20 and 0. Finally we 
make the last (6’th) column all zeros. Now from RREF A1420: 






Homogeneous Systems of Equations Calculator 1 .4 39 



[[1000 -.288096195186 0 ] 
[0100 - 1.00777740881 0 ] 
[0010 - 1.28332488596 0 ] 
[0001 - 1.59634374399 0 ]] 



we see the solutions are (,288096195186x 5> 1.00777740881x s , 1.28332488596x s , 1.59634374399x s , x s ) with 
x 5 arbitrary. 
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MATLAB 1.4 

1. (a) One example: (your answer will differ). 



» A = rand(3,4) 
A = 




0.0331 0.9554 0.8907 


0.1598 


0.5344 0.7483 0.6248 


0.2128 


0.4985 0.5546 0.8420 


0.7147 


0>) 




» rref(A) 




ans = 




1.0000 0 


0 0.3685 


0 1.0000 


0 -1.1232 



0 0 1.0000 1.3704 

(c) The solution of this system has x\ arbitrary, xi = — .3685x4, X 2 = 1.1232x4, and X 3 = —1.3704, 
since Ax = 0 is equilvalent to ans*x=0. The associated homogeneous equation has more un- 
knowns than equations, which gives a non-pivot column. Since a variable can be chosen arbitrar- 
ily, there are an infinite number of solutions, as predicted in Theorem 1. 

2. Most matrices with more rows than columns will have only one solution to the homogeneous system. 
However, it is not true that all of them do. The matrix in (ii), for example, does not. 

(i) 

» A = [ 1 2 3 0; -1 4 5 -1; 0 2 -6 2; 1 1 1 3; 0 2 0 1] 

A = 

12 3 0 

-1 4 5-1 

0 2-62 

1113 
0 2 0 1 
» rref(A) 
ans = 

10 0 0 

0 10 0 

0 0 10 

0 0 0 1 

0 0 0 0 

There is a unique solution, x = 0, to Ax = 0. 

(ii) 

» A =[ 1 -1 3; 2 1 3; 0 2 -2; 4 4 4] 

A = 

1-13 
2 13 

0 2 -2 

4 4 4 

» rrei(A) 
axis = 

10 2 

0 1 -1 

0 0 0 

0 0 0 

Since column 3 has no pivot, the solution will have x 3 arbitrary, xi = — 2 * 3 , and x 3 = 1 x 3 . There 
are an infinite number of solutions to the homogeneous equation. 
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3. ’’Balancing” will lead to a homogeneous system because any solution, x will lead to another solution 
by scaling rx. In order to get a unique solution, we must require that x is made of positive integers 
with no common divisor. 

(a) The solution for the example will be: 

» A = [ 1 0 -6 0 

21 - 6-2 
0 2 -12 0 ]; 

» R = rref(A) 

R = 

1.0000 0 0 -1.0000 

0 1.0000 0 -1.0000 

0 0 1.0000 -0.1667 

» format rat '/, Use this in version 4.0 in order to view 
'/, output as rational numbers. 

» z = R( : ,4) 

2 = 

-1 

-1 

- 1/6 

» format '/, This returns to the standard output format. 

From this, we can see that if £4 is chosen to be 6, then £3 = 1, £2 = 6 and xi = 6. With this 
choice, there will be no common divisors, and all of the variables are positive integers. 

(b) First we set up the equations, £1 and £2 will correspond to the compounds on the left and £3 
through xe for those on the right. Once we subtract the right from the left, all of the coefficients 
from the right will be negative: 

» format rat '/, as above, use rational numbers in the output. 



» A = [ 


1 


0 


-0 


-0 


-3 


0 




'/, Those nith 


Pb. 


3*2 


0 


-0 


-0 


-0 


-1 




'/, Those with 


N. 




0 


1 


-2 


-0 


-0 


-0 




/, Those with 


Cr. 




0 


2 


-0 


-1 


-0 


-0 




'/, Those with Mn. 




0 


4*2 


-3 


-2 


-4 


-1 


] 


V, Those with 


0 . 


A = 




















1 




0 






0 




0 


-3 


0 


6 




0 






0 




0 


0 


-1 


0 




1 




- 


-2 




0 


0 


0 


0 




2 






0 




-1 


0 


0 


0 




8 




- 


-3 




-2 


-4 


-1 


» rref(A) 




















ans = 




















1 




0 






0 




0 


0 


-1/6 


0 




1 






0 




0 


0 


-22/45 


0 




0 






1 




0 


0 


-11/45 


0 




0 






0 




1 


0 


-44/45 


0 




0 






0 




0 


1 


-1/18 



» format */, As above, it is a good idea to return to the default format. 

From the reduced echelon form of A, we can see that in order to make all of the variables inte- 
gers, xe must be chosen to be the least common multiple of 45, 6, and 18, which is 90. With this 
we get the answer: £1 = 15, £2 = 44, £3 = 22, £4 = 88, £5 = 5, £5 == 90. 
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11. (3, — 1, 4, 2) + (—2, 3, 1,5)= (1,2,5, 7) 12. (6, 0, -1, 4) - (3, -1, 4, 2) = (3, 1, -5, 2) 

13. 4(— 2, 3, 1,5)= (-8,12,4,20) 14. (-12,0, 2, -8) 15. (6, —2, 8, 4)— (—2, 3, 1, 5) = (8, —5, 7, — 1) 

16. (24, 0, -4, 16) - (21, -7, 28, 14) = (3, 7, -32, 2) 

17. (7,2,4,11) 18. (-2,1,10,5) 19. (-11,9,18,18) 

20. (3a + 6/? — 27, —a + 37, 4a — p + 7, 2a + 4/3 + 57) 




/ — 2 0 \ / 00 \ 

25. 0 14 = 00 

\ — 7 5/ \00/ 
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/ 13\ / — 2 0\ / —7 —21 \ { -6 0\ / — 13 —21 \ 

26.-7 2 5 + 3 14 = -14-35 + 3 12 = -11 -23 

V— 1 2 / V- 75 / V 7-14/ \— 21 15/ V— 14 1/ 



27. 




-2 4\ 
7 15 

-15 10 / 



28. 




-2 

-3 

-4 



29. 2 





+ 4 




2 6\ / -6 0\ / -4 4\ 

4 10 - 3 12 + 16 24 

-2 4/ \ —21 15/ \ —28 12/ 



4 10\ 
17 22 



30. 7 



( 



-1 1 
4 6 
-7 3 



-2 0 
1 4 
-7 5 




( ~ 1 7 \ /~ 20 \ ( 2 6 \ 

= 28 42 - 1 4 + 4 10 

V —49 21 / \— 7 5/ \ —2 4 / 



31. 2A + B — D = O, says D = 2A + B = 



0 6 \ 
5 14 
-9 9/ 



32. A + 2B — 3C + E = O says £ = 3C - 25 - A = 




/ 1 — 1 2\ /0 2 1\ / 1 — 1 2\ / 0 4 2 

33. 3 4 5 -2 3 05 = 3 4 5 - 6 0 10 

\0 1—1/ \ 7 — 6 0 / \0 l-l/ \ 14 -12 0 



1 -5 0\ 

-3 4-5 
-14 13 -1 ) 



/l-l 2\ / 0 0 2\ / 3 -3 6\ /0 0 2\ /3 -3 4\ 

34.3 3 4 5 - 3 10 = 9 12 15 - 3 1 0 = 6 11 15 

\0 1 - 1 / \0-24/ V° 3 “ 3 / V° —2 4/ \o 5 - 7 / 



/I -1 2 \ /0 2 1 \ /0 0 2 \ /I 1 5\ 

3 4 5 + 3 05 + 3 10 = 9 5 10 

Vo l-l/ \7-60) V° — 2 4 / V 7 “ 7 3 / 



35 . 




44 Chapter 1 Systems of Linear Equations and Matrices 



Instructor’s Manual 





/0 0 2 \ 

+ 2 3 10 

\0 -2 4 / 

0 2 l\ / 0 0 4\ l 2— 47\ 
3 05 + 6 20 = 9 10 5 

7-60/ \0-44/ \ —7 4 2/ 




- 1-1 - 1 \ 
-3 -3 -10 
-7 3 5/ 



/0 0 2\ /0 2 l\ / 1 — 1 2\ 

38. 4 3 10 -2 3 05 +3 3 4 5 

\0-24 / V 7 — 6 0/ \0 1-1/ 



/ 0 0 8\ / 0 4 2\ / 3 — 3 6 

12 4 0 1- 6 0 10 + 9 12 15 

\ 0 -8 16/ \ 14 -12 0/ V° 3 ~ 3 



3 -7 12 \ 
15 16 5 
-14 7 13/ 



39. D — —A — B — C 




40. A + 25 - 3C - 4E 1 = O says 417 = 3(7 - 2J3 + 8A Divide by 4 to get 17 = | [3(7 - 2B + 8A] 




0 4 2\ / 8-8 16\ 

6 0 10 + 24 32 40 

14 -12 0/ V 0 8 -8/ 



/ 2-3 5\ 

27/4 35/4 15/2 
\ -7/2 7/2 1 / 



an a i2 • • ■ ajn 



41. 0 




O an O-an 0-a ln \ /0 O' 



= O 



s a m i a m 2 • • • a mn / \ 0 • a m i 0 • a m 2 • • ■ 0 • a mn / \ 0 • • • 0 / 

'0 ••• 0\ / an a i 2 ••• ai„\ / 0 + an 0 + ai 2 0 + ai„' 

v 0 • ■ • 0 / \ a m i a m 2 • • • a m „ / V 0 + a m i 0 + a m2 • • • 0 + a„ 

an ai2 ■ • • ai„ \ 

a '/ 

^ml ^m2 ' ‘ ’ Q-mn ' 

1 • an 1 • ai2 ■ • • l-ai„\ / an ai2 • • • ai„ ' 



an ai2 • • • ain 



< a m 1 a m2 * * * a mn 




1 * ^ml 1 ‘ ^m2 * 1 




= A 



a ml a m2 ’ ‘ ‘ a mn > 
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42 . ( A + B) + C = 



an a 12 ••• a i„\ / bn 612 ■■■ b i„' 

= ) + = 

y. a m i dm2 ‘ * ‘ a mn / V ^ml ^m2 * * * ^mn - 



+ 



C 11 C 12 • • • Ci n ' 



^ c ml c m2 ' ' ' c n 



( (an + 6n) + cn (ai2 + 612) + C12 • • • (ain + bi n ) + ci„ ^ 

\ ( a ml. + b m 1 ) -f“ C m l (®m2 4" ^mT.) 4" Cm2 * * * (®mn 4“ ^mn) 4* C m n / 
( an 4- (611 4- Cn) ai2 4- (612 4- C12) • • ■ a\ n + (6i„ -f ci„) \ 

\ a m i 4- (6mi 4- Cmi) «m 2 4- ( 6 m 2 4" Cm 2 ) • ■ • a m « 4- (6 mn ~f" C mn ) y 
= A + (B + C) 



43 . If A = ( aij ) and B = (6<y), then a(A + B) = a((a,y) + (6 fj -)) = a(a,j 4- 6,y) = (a(a ij - + 6,y)). 
a /1 + aB = a (aij) + a(6, ; ) = (aa i; ) + (a6,j) = (aa,j + a6,j) = (a(a ij - + &;_,)), and 
Therefore a(.A + B) — aA + aB. Similarly, 

(a + p)A = (a 4- /?)(a,j) = ((a + /?)a ti ) = (aa.y + /?a (j ). 
ayl + /?T = a(aij) + /?(a tJ ) = (aa.y) + (/?a<j) = (qa,,) + /?a, ; ) 

Therefore (a + / 3 )A = a A 4 - PA. 



/0 1 1 0\ 
110 10 
110 1 
\0 0 1 0/ 



45 . 



/0 1 0 1 0\ 
10 110 
0 10 0 1 
110 0 0 
\0 0 1 0 0/ 



46 . The entries of d 4- e represent the demand for all four raw materials if each factory is to produce 1 
unit. 2d gives the total raw material needs for factory 1 to produce 2 units. 
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1. (a) 

» A = [ 1 2-2 0 1 

24-10-4 
-3 -6 12 2 -12 

1 2-2-4 -5 ] 

A = 

1 2-20 1 

2 4-1 0-4 

-3 -6 12 2 -12 

1 2-2-4 -5 

» '/, Choose A (1,1) as the first pivot. 

» c = -A(2, 1)/A(1 , 1) ; A(2, : ) = A(2,:) + c*A(l,:); 

» c = -A(3 , 1)/A(1 , 1) ; A(3, : ) = A(3, : ) + c*A(l,:); 

» c = -A(4,l)/A(l,l); A(4, : ) = A(4,:) + c*A(l,:) 

A = 

1 2-20 1 

0 0 3 0 -6 

0 0 6 2 -9 

0 0 0 -4 -6 

» */. How let A(2,3) be the pivot. 

» c = -A(3,3)/A(2,3) ; A(3,:) = A(3,:) + c*A(2,:) 

A = 

1 2-2 0 1 

0 0 3 0 -6 

0 0 0 2 3 

0 0 0 -4 -6 

» */, Now let A (3, 4) be the pivot. 

» c = -A(4,4)/A(3,4) ; A(4,:) = A(4,:) + c*A(3,:) 

A = 

1 2-20 1 

0 0 3 0 -6 

0 0 0 2 3 

0 0 0 0 0 

This matrix is now in echelon form. 

(b) Recall that rand generates a random matrix, so your answer will be different. 

» A = rand(4,5) ; 

» A ( : ,3) = 2*A( : , 1) + 4*A(:,2) 

A = 

0.5269 0.7012 3.8586 0.7564 0.9826 

0.0920 0.9103 3.8252 0.9910 0.7227 

0.6539 0.7622 4.3566 0.3653 0.7534 

0.4160 0.2625 1.8818 0.2470 0.6515 

» !, Choose A (1,1) as the first pivot. 

» c = -A(2, 1)/A(1 , 1) ; A(2, : ) = A(2,:) + c*A(l,:); 

» c = -A(3 , 1)/A(1 , 1) ; A(3, : ) = A(3,:) + c*A(l,:); 

» c = -A(4,l)/A(l,l); A(4, : ) = A(4,:) + c*A(i,:) 
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A = 



0.5269 


0.7012 


3 . 8586 


0.7564 


0.9826 


0.0000 


0.7879 


3.1518 


0.8590 


0.5512 


0 


-0.1080 


-0.4319 


-0.5734 


-0.4660 


0 


-0.2911 


-1.1645 


-0.3501 


-0.1242 


i Now let A(2,2) be 


the pivot 






; = -A(3,2)/A(2,2) ; 


A(3, : ) = 


A ( 3 , : ) + 


c*A(2, :); 


; = -A(4,2)/A(2,2); 


A(4, : ) = 


A(4, : ) + 


c*A(2, : ) 


0.5269 


0.7012 


3.8586 


0.7564 


0.9826 


0.0000 


0.7879 


3.1518 


0.8590 


0.5512 


0.0000 


0 


0.0000 


-0.4556 


-0.3905 


0.0000 


0 


0 


-0.0327 


0.0795 


1 Now choose A(3,4) 


as the pivot. 




= -A(4,4)/A(3,4); 


A(4, : ) = 


A(4, : ) + 


c*A(3, : ) 


0.5269 


0.7012 


3.8586 


0.7564 


0.9826 


0.0000 


0.7879 


3.1518 


0.8590 


0.5512 


0.0000 


0 


0.0000 


-0.4556 


-0.3905 


0.0000 


0 


0.0000 


0 


0.1075 



It is worth noticing that we have introduced some small round-off error. For example: 

» A(2, 1) 
ans = 

1 . 3878e-17 



This is not exactly zero. However, up to the 14 significant digits that MATLAB can print out, 
the matrix we have above is in row echelon form. 



2. (a) 



» a = zeros (5); 

» a(l , [2 4]) = [1 1]; 

» a(2, [1 3 4]) = [1 1 1]; 
» a(3, [2 5]) = [1 1]; 

» a(4,[l 2]) = [1 1]; 

» a(5,3) = 1 



0 10 1 
10 11 
0 10 0 
110 0 
0 0 10 



0 

0 

1 

0 

0 



(b) A will have 5 rows, since there are 5 nodes, and 8 columns, since there are 8 edges. 



» A = zeros(5,8) ; 

» A ( [1 2], 1) = [-1; 1] 

» A([2 4], 2) = [-1; 1] 

» A( [4 1], 3) = [-1; 1] 

» A( [1 3], 4) = [-1; 1] 

» A ( [3 5], 5) = [-1; 1] 

» A ([5 1], 6) = [-1; 1] 

» A([3 4], 7) = [-1; 1] 

» A( [5 4], 8) = [-1; 1] 
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-l o l-l 

l-ioo 

0 0 0 1 

0 1-10 

oooo 

3. (a) 

» A = [ 1 2 3; 4 5 6]; 7.2 x 3 
» B = [ 1 2; 3 4; 5 6] ; */.3 x 2 
» A + B 

??? Error using ==> + 

Matrix dimensions must agree. 

(b) In general, if A and B are the same size, 

sA + sB = s(yl + B ). 

We can check this on random matrices as follows: 



» A = rand(3) 




A = 


0.9866 0.0907 


0.5007 


0.4940 0.9478 


0.3841 


0.2661 0.0737 

» B = rand(3) 


0.2771 


B = 


0.9138 0.9410 


0.7702 


0.5297 0.0501 


0.8278 


0.4644 0.7615 

» s = rand(l) 


0.1254 


s = 

0.0159 


» C = s*A + s*B 




C = 


0.0302 0.0164 


0.0202 


0.0162 0.0158 


0.0192 


0.0116 0.0133 

» D = s* (A+B) 


0.0064 


D = 


0.0302 0.0164 


0.0202 


0.0162 0.0158 


0.0192 


0.0116 0.0133 


0.0064 


» C - D 
ans = 


1.0e-17 * 




0.3469 0 


0 


0 0 


0 


0 0 


0 



We see that C — D is not exactly zero, but the difference between C and D can be accounted for 
by round off error in the computer. This process can be repeated with any matrices A and B, 
and the same results will occur. 



0 10 0 

0 0 0 0 
-1 0-1 0 

0 0 11 

1-1 0-1 




Section 1.6 

1. 2-3 + 3-0 + (-5) -4 = -14 
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2. 1 • 3 + 2 • (— 7 ) + (— 1 ) *4 + 0 • (—2) = — 15 
3 . 5 • 3 + 7 • (—2) = 1 4 . 8 - 7 + 3 -(- 4 ) + 1-3 = 47 5 . ac+bd 

6. xy + yz+zx 7 . (— l) 2 + (— 3) 2 + 4 2 + 5 2 = 51 

8. Since a 2 > 0 then a • a = a 2 + a§ + 1- a 2 > 0 

9. If a = 0, then a • a = 0. Conversely, if a y 0, then a 2 > 0 for some i, and hence a • a > 0. Thus, a = 0 
if and only if a • a = 0. 




( 8 + 0 2+ 18 \ _ ( 8 20 \ 
V— 4 + 0 -1+ 12 ) ~ 11 ) 



f-15-2 18 — 6 \ _ / —17 12 \ 
V -5 + 4 6 + 12; _ v -1 is; 



17. 



(-1 -2 0-3\ _ f-3 —3 \ 
1+20 + 3; _ V 1 3; 



f - 15 + 6 10 + 24\ _ f-9 34 \ 

V 3 + 3-2 + 12;“^ 6 10 y 



/— 12 + 25 + 0 4 + 30 + 1 —4 + 20 + 2\ _ / 13 35 18\ 
V 0 + 20 + 0 0 + 24 + 2 0+ 16 + 4^ _ ^20 26 20; 



f 7+ 0-8 42+4+ 12\ _ /— 1 58N 
^2 + 0-10 12 -12 + 15 J ~ V-8 15 ) 



( 7 + 12 1- 18 4 + 30\ / 19 —17 34 \ 

21. 0 + 8 0- 12 0 + 20 = 8 -12 20 22. not defined 

\ — 14 + 6 —2 — 9 — 8 + 15 / \ —8 —11 1 ) 



/ 2 + 4+ 12 -3 + 0 + 18 5 + 24 + 6 \ / 18 15 35 \ 

23. -4 + 3+ 10 6 + 0 + 15 -10 + 18 + 5 = 9 21 13 

\ 2 + 0 + 8-3 + 0 + 12 5 + 0 + 4; \10 9 9/ 
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2 + 6 + 5 8-9 + 0 12- 15 + 20 
1 + 0 + 64 + 0 + 0 6 + 0 + 24 

2-6 + 18 + 9 + 0 12 + 15 + 4 



13 -1 17 
7 4 30 
-3 17 31 



25. (3 + 8 + 0 - 4 



-6 + 16 + 0 + 6 ) = (7 16) 



3 + 8 + 0 — 4 
-6 + 16 + 0 + 6 



)-U) 



3 -2 1 

4 0 6 

5 1 9 



/ 3 -2 1 
28. 4 0 6 

\ 5 19 



( a*l-f 6*0 + c *0 fl*04"6’l*4*C"0 ct • 0 -f - • 0 — }~ c • 1 
d • 1 -{■ c • 0 -f- / •Od’O'fC’l'f 
^•l + A'O + i- O.flfO + A'l + i’O flfO + A'O + i*! 



orj Txr / 2a 4* 6 3a - 1 - 26 

30. We want 2e + d3c + 2J 

2a -}- 6 — 1 2c -f* d — 0 
3a -j" 26 = 0 3c 2d — 1 



j r ^ o j ) — ( n ? ) • This gives two systems of equations: 
a 6c 2d J yuly 

+ 2d = r SoIvingforaand6wefind (3 2 |o) _+ (o 1/2 | -3/2) 

Hence, a = 2 and 6 = —3. Similarly, one can show c = — 1 and d = 2 



1 0 2 
0 1 -3 



31. We want | 
equations: 



A Oll&ll 
\ 021&11 



+ 012^12 au &12 + 012^22 

+ Ct 22^21 <*21^12 + 022^22 



- (”)■ 



ail&ll + <112^21 = 1 
021 ill + a 22^21 = 0 



. As in problem 30, this gives two systems of 

«11^12 + “12^22 = 0 
021&12 + U22&22 = 1 



Since 011022 — ai2«2i i 1 0, then either an or a 2 i is nonzero. Without loss of generality, we may assume 



an 0. Solving for in and 621 we obtain 



/ an 


012 




\ o 2 i 


022 


0/ 



0 022 — a 2 iai 2 /an — 021/011 



Solving for fe 2 i gives 621 = — 021/(^11022 — 021012)- Use back substitution to find in = 
022/(011022 — 021012). Similarly, solving the second set of equations for 612, 622, gives 
612 = — 012/(011022 — 021012), and 622 = on/(ana 2 2 ~ 021012). 



33. (a) There are 3 people in group 1, 4 in group 2, and 5 in group 3. 

/ 2 1 1 0 0\ 

(b) AB =(ll010J 

V 1 0 2 0 1 / 
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34. (a) There are 2 people in group 1, 5 in group 2, and 7 in group 3. 



(b) 



AB = 



/2 1 0 1 2 1 3\ 
V0202101/ 



35. a • b = 2 • 3 + (-3) -2 = 0 orthogonal 

36. a • b = 2 • (-3) + (-3) • 2 = -12 ^ 0 not orthogonal. 

37. a* b = 1 • 2 + 4 • 3 + (-7) -2 = 0 orthogonal. 

38. a-b=l -0 + 0- l+ l*0 + 0- l = 0 orthogonal. 

39. a-b = a- 0 + 0 - d+ 6*0 + 0 - e + c- 0 = 0 orthogonal. 

40. We want (1, -2, 3, 5) • (-4, a, 6 , -1) = -4 - 2a + 18 - 5 = 0. Hence, a = 9/2. 



41. We want 




= 4 — 5a — 4/? + 21 = 0. Let /? be arbitrary. Then a = 5 — 



42. (i) a-0 = ai • 0 + a 2 • 0 H |-a„-0 = 0 (ii) in text 

(iii) a • (b + c) = a • (6i + c\ , b 2 + C 2 , • • • , b n + c„) 

= ai(5i + ci) + 02(62 + C2) + • • • + o n (6 n + c n ) 

= ai6i + 0262 + ' • • + a„6„ + OlCi + 02 c 2 + ■ • • + dn c n 



= a- b + a- c 

(iv) (aa) • b = (aai, 002, • ■ aa n ) • b 

= aa 1 6 1 + aa 2 b 2 + h aa n b n 

— a(ai6i + 0262 + • ■ • + a n b n ) 
= a(a • b) 



43. (a) (2, 3, 5,1) 




(c) total hours = 2 + 4.5 + 2.5 + 2 = 11 



/ 0.055 \ 
1.80 
0.20 
0.001 
V 0.40/ 



44. (a) (1000,20,100,5000,50) (b) 



(c) $136.00 



Oi|^ 
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45. 




80,000 45,000 40, 00(A 
50 20 10 ) 





255, 000 \ 
120 ) 



46. Express the sales of each item in each month as the 4x3 matrix A = 



/4 2 20' 

5 3 12 | ' E x P ress the unit 

\8 2.5 20, 



3.5 1.5 

profit and unit taxes as the 3x2 matrix B — I 2.75 2 I . Then AB = 



1.5 0.6, 



/ 49.5 22 \ 

37.25 16.4 
43.75 20.7 
V 64.875 29/ 



shows the 



total profits and taxes in each of the four months. 



47. 



/ 4 — 4 -2-6 
1^8 + 24 -4 + 36 




49. A 2 





f 11 38 \ 
V57 106 ) 



50. A 2 = 



( 0 0 1 0 \ 
0 0 0 1 
0 0 0 0 
0 0 0 0 / 



( 0 0 0 1 \ 
0 0 0 0 
0 0 0 0 
0 0 0 0 / 



A 4 = A 5 = 0 



/° 


0 


1 


0 






/° 


0 


0 


1 






(° 


0 


0 


0 


!\ 


0 


0 


0 


1 


0 




0 


0 


0 


0 


1 




0 


0 


0 


0 


0 


0 


0 


0 


0 


1 


yl 3 = 


0 


0 


0 


0 


0 


A 4 = 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


Vo 


0 


0 


0 


0/ 




Vo 


0 


0 


0 


0/ 




Vo 


0 


0 


0 


0/ 



52. Let a,ij be an element of A. 
AB — C . C — 0 implies c,i 



Let B be an n x n matrix with bj\ = 1 and 0’s everywhere else. Let 
= a.ij = 0. Since a,j was arbitrary, then + is the zero matrix. 



53. PQ = 



/ ii 

/ 90 

11 



120 120 60 



V I 




Each component is > 0 and the sum of the elements in each row is 1. 



54. Clearly, the elements of P 2 are positive. If P is n x n, let v be a column n-vector with 1 as its ele- 
ments. Note that if A is an n x n matrix, then the sum of the elements in each row of A is 1 if and 
only if Av = v. We have P 2 v = P(Pv) = P\ = v. Hence, P 2 is a probability matrix. 



55. Since every entry of P and of Q is positive, the entries of PQ are all positive. If p and Q are n x n 
matrices, let v be the column n-vector with every element 1; note that if A is an n x n matrix, the 
sum of its entries along any row is 1 if and only if j4v = v. Since ( PQ)v = P(Qv) = Pv = v, we see 
that PQ is a probability matrix. 
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56. ABCD = AB(CD) = A(B(CD)) = A{BC)D = (AB)CD = ((. AB)C)D = ( AB)(CD ) 



57. R 2 = 



0 0 11 



2 0 1 0 Si = 1 + i(l + 1) = 2 S 2 = 1 + 1 + |(2 + 1) = 3.5 



110 0 



0 100 S 3 = 1 + |(1)= 1.5 5 4 = 1 + 1 + |(1 + 1) = 3 



(a) S 2 > S 4 > Si > S 3 (b) The score is the number of games won by player i plus half the 
number of games won by players who player i beat. 

n 

58. OA = 0\. The ij th component of OA, c,j , can be written Cij — ^ onaij . Since each element of O is 

k = 1 

0, then Cij — 0. Hence, 0 1 is the m x p zero matrix. 



5,.( S + C , = W -(«“)=^C=(«») + (3») 
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65. AB = 



( I 0\ ( I 0\ _ ( I 2 + OD I-0 + 0-l\ _ ( I 2 0\ 
\C I ) \D I J ~ \CI + ID C-O + I 2 J ~ \C + D I 2 J 



( I 0\ ( I 0\ _ f I 2 + OC i-o + o-i\ 
\D I J \C I J ~ \DI + IC D-O + I 2 ) 



1 2 O 
D + C I 2 



I 2 0 
C + D I 2 
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74. 2° + 2 1 + 2 2 + 2 3 + 2 4 = ^2* 

k=zO 



5 

75. (-3)° + (— 3) 1 + (— 3) 2 + (— 3) 3 + (— 3) 4 + (-3) 5 = ^(-3)*' 

1=0 



n n 7 9 4 

76 . ^T^fc/(fc + l) 77 . X/ 1/f 78 . XV* 79 . ^(-l) i+1 o-*' as = 8 + 27 + 

ifc = 2 1 = 1 *=0 i=0 7 = 2 

64 = 99 . 

8 

80 . (2 • 1 - 1)(2 • 1 + 1 ) + • • • + (2 • 8 - 1)(2 • 8 + 1 ) = ] T )( 2 * - l )( 2 Jb + 1 ) 

k = l 



7 7 

81. 2 2 (2 • 2) + 3 2 (2 • 3) + • • • + 7 2 (2 • 7) = X> 2 ( 2 *) = X> 3 

Jfc = 2 i = 2 



3 7 3 2 3 2 

82. + a 2 j) = X^XXij 83. X^aij + a 2 j + a 3j ) = ^ Xl a u 

j— l 1=17=1 i=i i=i j=i 



4 4 4 5 

84. ^(a 2 j + a 3j + a 4j ) = ^ ^2 a ij 85 ‘ X! a3 * 6 * 2 

7=1 1=2 7=1 i=l 



4 3 4 

86. ^](a2l6l7'C7'5 + a22&27 c 7'5 + <*23637 C75) = X] X/ j2l '^i C J 5 

7 = 1 1 = 1 7 = 1 

N 

87. 2(°* + 6 *) = ajv/ + 6 at + aM+i + 6 m+i + ■ ■ • + ajv + 6w 
k=M 



— “M + “M+l + • • ' + + 6 m + 6 m + 1 + • • • + 677 



JV AT 



= Yl ak ~ ^2 bk 

k=M k=M 



88. X^ (a/c — bk) — X^ ( a * + (~1)M ^ X/ ajk X/ ( — 0^t ^ ak ~ X] ^ 

i=M k=M k=M k=M k=M k=M 
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N 

89. ^ ak = a\f + a.M + i + • • • + 

t=JU 

= a\f + a,M + 1 -f • • • + cim-i + d m + a m +i + 1- a;v 

m N 

= J2 ak+ J2 ak 

k—M k=m + l 
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CALCULATOR SOLUTIONS 1.6 

In this section many of the problems have two input matrices, and we append the problem number to the matrix 
name; for instance we use A16nn and B16nn to refer to the left and right factors in the products in problem nn, 
nn=90,91,92. We assume that the matrices have been entered and we will only show how to obtain the solutions 
from these input matrices. 



[ [ -12.0704 -26.8412 ] 
90. A1690 (x) B1690 (ENTER) gives the product: [ 44.4207 -45.0695 ] 

[ 91.7485 4.2242 ]] 



91. A1691 (x) B1691 | ENTER) gives the product: 



t 34192 38621 ] 

[ 50408 44115 ] 

[ 62661 71731 ] 

[ 59190 55046 J] 



[f -.6557 -3.3655 ] 

92. A1692 (x) B1692 (ENTER) gives the product: [ .6907 3.4072 ] 

[ -1.4255 5.5459 ] ] 



93. (a) To see PI 69 3 and Q1693 are probability matrices compute their row sums by computing their products 
with 

the column of all ones: [[1][1][1][1]] (STO) ONES 

[[ 1 ] [[ 1 ] 

P1693 (x) ONES (ENTER) : and Q1693 (x) ONES (ENTER) : . 

[ 1 ]] [ 1 ]] 



The product P1693 (x) Q1693 (ENTER) : ^ 

[ 

can be saved as PQ1693 via (2nd) (ANS) (STO ) 



by showing its row sums are all 1 by PQ 1 6 9 3 (x) 



.31118 .18444 .14174 .36264 ] 

.32625 .27585 .08454 .31336 ] 

.17955 .22651 .19619 .39775 ] 

.30047 .15251 .33558 .21144 ]] 

PQ1693. Then you see the product is a probability matrix 

[[ 1 ] 

ONES (ENTER] which yields: j . 

[1] ] 



94. Entering A1 69 4 Q n, n = 2, 5, 10, 50, 100 gives: 



A1694 2 

is 



A1694 5 

is 



A1694 10 

is 



[[1 93] [[1 3069] 

[0 32]] [0 1024]] 





A1694 50 




A1694 100 




is 




is 


[ [1 


3 . 378e15] 


t [1 


3 . 803e3 0 ] 


[0 


1. 126e15] ] 


to 


1. 268e30] ] 



[ [1 9] 

[0 4] ] 
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95. From the A 2 , A 5 , A 10 results in Problem 94, where the diagonal components are {l 2 , 2 2 }, {l 5 , 2 5 } and 
{ 1 1 0 , 2 1 0 } , it appears that the diagonal components of A n are just a n , b n , c n . 
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MATLAB 1.6 



» A = rand(3,4) 

A = 

0.6885 0.7362 0.8886 0.3510 

0.8682 0.7254 0.2332 0.5133 

0.6295 0.9995 0.3063 0.5911 

» B = rand(4,2) 

B = 

0.8460 0.4154 

0.4121 0.5373 

0.8415 0.4679 

0.2693 0.2872 

» A*B 
ans = 

1.7281 1.1982 

1.3679 1.0070 

1.3614 1.1116 

» B*A 

??? Error using ==> * 

Inner matrix dimensions must agree. 

The product of a 3 x 4 with a 4 x 2 will be a 3 x 2 matrix, so AB is a 3 x 2 matrix. However, the 
product of a 4 x 2 with a 3 x 4 is not defined since 2, the number of columns of left factor, is not 3, 
the number of rows of the right factor. 



» A = round(10*(2*rand(3)-l) ) 

A = 

-9 -10 -2 

1-2 4 

3-9 2 

» B = round(10*(2*rand(3)-l)) 

B = 

9-8 4 

7 3 8 

1-2 5 

» A*B 
ans = 

-153 46 -126 

-1 -22 8 

-34 -55 -50 

» B*A 
ans = 

-77 -110 -42 

-36 -148 14 

4 -51 0 

The probability that for two random matrices, AB = BA is very small. 
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3. (a) 



» A = [ 2 9 -23 0 
0 4 -12 4 
7 5-11 
7 8 -10 4 ] ; 



» b = [ 34; 24; 15; 33 ] ; 

» z = C -2; 3; 1; 0] ; 

» x = [ -5; 10; 2; 2] ; 

» A*x '/. This will be b. 
ans = 

34 

24 

15 

33 

» A*z '/, This will be zero, 
ans = 

0 

0 

0 

0 



(b) For any scalar s, 



yl(x + sz) = .Ax + sAz = Ax. + 0 = Ax. 



This can be tested by computing the following: 



»s = rand(l) 
s = 

0.7529 

» A * (x+s*z) */, This will be b 

ans = 

34.0000 

24.0000 

15.0000 

33.0000 



4. 



»'/, Generate Random matrices: 

» A = round ( 10* ( 2*rand(2,4)-l) ) 
A = 

-3 3 0 4 

-5 4 5 9 



» B = round( 10* ( 
B = 

-1 0 10 

9 3-7 

-4 -2 3 

-1 2 -3 



2*rand(4,5)-l) ) 

1 -5 

0 3 

7 1 

-10 8 




Vector and Matrix Products MATLAB 1.6 



61 



» B( : ,2) 


= B( 


:,3) 






B = 


-1 


10 


10 


1 


-5 


9 


-7 


-7 


0 


3 


-4 


3 


3 


7 


1 


-1 

» '/, part 


-3 

(b): 


-3 


-10 


8 


» A*B 
ans = 


26 


-63 


-63 


-43 


56 


12 


-90 


-90 


-60 


114 



Since columns 2 and 3 were the same in B, they will also be the same in AB. 

(c) The above can be repeated several times. 

(d) Proof: Assume that the columns m and n are the same in B. This implies that 6,- m = for 
any i. If C = AB, what we wish to show is that c jm = cj n for any j, i.e., that the m and nth 
columns of C are the same. To do this, we will write cj m using J2 notation. Assume that the 
width of A is p, then: 

p p 

Cjm = ^ ] Qjjbjm = 'y ] Cljjbjn = Cj n . 
i = 1 i=l 

Which concludes the proof. 



» A = 
A = 


round ( 


10*( 


2*rand(5, 


6) - 


- 1) ) 


-5 


-1 


4 


4 


0 


5 


9 


6 


-1 


-4 


-8 


-1 


-5 


-10 


2 


5 


-7 


-10 


7 


-7 


5 


-7 


8 


2 


4 


4 


10 


-10 


6 


0 


II 

X 

A 

A 


round ( 


10*( 


2*rand(6 


,D 


- 1) ) 



x = 

-9 

1 

-2 

-6 

-2 

3 

» A*x - ( x(l)*A(:,l) + x(2)*A(:,2) + x(3)*A(:,3) + ... 
x(4)*A(:,4) + x(5)*A( : ,5) + x(6)*A(:,6) ) 

ans = 

0 

0 

0 

0 

0 

This is essentially the same as expression (10) in the text. The expression inside the parentheses is 
the definition of matrix multiplication. 
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6. (a) We can write AB = BA as AB — BA = 0 . Then, with this choice of A and B, 



0 = AB- BA 

_ f ax i + bx 3 ax 2 + bx 4 
y cx 1 + dx 3 cx 2 + dx 4 



\ / aari + cx 2 + dx 2 \ 

y y 0x3 + CX4 6x3 -f c?X4 y 






_L U 



-CX2 + bx 3 —6xj + (a — d)x 2 + 6x4 

\ I . 



If we use each of the 4 entries in this matrix as one equation in our system, we will get a 4 x 4 
system with coefficient matrix R and variables x. 

(b) 0) 



» a 

» R 



R = 



1 ; b = — 1 ; c=5;d=-4; 
[ 0 -c b 0 

-b a-d 0 b 

c 0 d-a -c 

0 c -b 0 ] 



0-5-1 0 

15 0-1 

5 0-5-5 

0 5 10 



» rref(R) 
sms = 

1.0000 

0 

0 

0 



0 


- 1.0000 


- 1.0000 


1.0000 


0 . 2000 


0 


0 


0 


0 


0 


0 


0 



So Xi = x 3 + x 4 , x 2 = |a:3 or B = ^ X3 + ■'■4 5^3 'j _ ^ ^1 1 / 5 ^ + ^ ^1 0 ^ jy ot i ce 

that the matrices *4 = x 4 I commute with all matrices, and so there will always be 

infinitely many solutions for any A. 



» format rat */, Use rat(rref(R), ’s’) in Hatlab 3.5 
» rref(R) 
ans = 



1 

0 

0 

0 



0 

1 

0 

0 



-1 

1/5 

0 

0 



-1 

0 

0 

0 



If we choose X3 = 5 , and X4 = 1 we will get Xi = 6, and x 2 = — 1 . You may choose any other 
integers, as long as X3 is divisible by 5 . 

(iii) The matrix B will be 



» B = [6 -1; 5 1] 
B = 



6 

5 



-1 

1 
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» A = [a b; c d] 

A = 

1 -1 

5 -4 

» A*B - B*A 
axis = 

0 0 

0 0 

Since AB — BA is zero, we have verified that AB = BA. 

(iv) This can be repeated for any other choice of x 3 and x 4 . 

(c) We can repeat the above, using the new matrix A: 

» a = 1; b = 2; c = 3; d = 4; 

» R = C 0 -c b 0 

-b a-d 0 b 

c 0 d-a -c 
0 c -b 0] 

R = 

0-320 
-2-3 0 2 

3 0 3 -3 

0 3-20 

» format rat */. This gives rational numbers in the output, for Matlab 4.0. 

» rref(R) '/, Use rat(rref(R), l’l) in Matlab 3.5. 
ans = 

10 1-1 

0 1 -2/3 0 

0 0 0 0 

0 0 0 0 

As above, we may choose x 3 and x 4 arbitrarily and B = x 3 

possible B’s.) If we choose x 3 = 3, and x 4 — 1 we will get x 4 
will be 

» B = [-2 2; 3 1] 

B = 

-2 2 

3 1 

» A = [a b; c d] 

A = 

1 2 

3 4 

» A*B - B*A 
ans = 

0 0 

0 0 



i oj + X4 [0 1 j ' ( A S ain X/tI are 

= —2, and £2 = 2. The matrix B 



Since AB — BA is zero, we have verified that AB = BA. (d) The above may be repeated using 
any matrix for A. To avoid round off error, you should use an integer matrix. 
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» A = round( 10*( 2*rand(2,2) -1) ) 

A = 

-6 4 

-9 4 

» B = round( 10*( 2*rand(2,2) -1) ) 

B = 

9 0 

-2 7 

» C = (A+B) "2 
C = 

-35 56 

-154 77 

» D = A'2 + 2*A*B + B'2 
D = 

-43 48 

-192 85 

» C-D 
ans = 

8 8 

38 -8 

In general, it is not true that C — D. However, they will be equal when we use A and B from prob- 
lem 6: 

» A = [1-1; 5 -4] ; 

» B = [6 -1; 5 1]; 

» C = (A+B) “2 
C = 

29 -8 

40 -11 

» D = A'2 + 2*A*B + B~2 

D = 

29 -8 

40 -11 

» C-D 
ans = 

0 0 

0 0 

When this is repeated with the matrices from 6(c), C will again be the same as D. In fact the state- 
ment 

{A + B) 2 = A 2 + 2 AB + B 2 

if and only if 

AB = BA. 

Proof: We may expand ( A + B) 2 as follows: 

{A + B) 2 = (A + B){A + B) 

= A(A + B) + B(A + B ) 

= AA + AB + BA + BB 

= A 2 + AB + BA + B 2 . 

If we subtract this from A 2 + 2 AB + B 2 , we get AB — BA, which is zero whenever AB = BA. Thus 
we may say that ( A + B ) 2 is A 2 + 2 AB + B 2 exactly when AB is BA. 
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8. (a) 



» A = round( 10*(2*rand(6,5)-l)) 
A = 



4 




3 




-9 


-3 


-9 


2 




-2 




5 


-5 


3 


9 




4 




-3 


10 


8 


7 




8 




3 


4 


-5 


1 




5 




5 


5 


-1 


-8 




-5 




10 


3 


5 


» E = 
E = 


[1 


0 0 


0 


0 0] 






1 




0 




0 


0 


0 


» E*A 














ans = 














4 




3 




-9 


-3 


-9 


» E = 
E = 


[0 


0 1 


0 


0 0] 






0 




0 




1 


0 


0 


» E*A 














ans = 














9 




4 




-3 


10 


8 



In the first case, EA was the first row of A, in the second case, EA was the third row of A. In 
general if E is all zeros except a 1 in the ith column, EA will be the ith row of A. 



» E = [2 0 0 0 0 0]; 

» E*A 
ans = 

8 6 -18 -6 -18 

» E = [002000]; 

» E*A 
ans = 

18 8 -6 20 16 

As above, EA will be the ith row of A, but this time it will be multiplied by 2. 

(c) (0 

»E = [1 01000]; 

» E*A 
ans = 

13 7 -12 7 -1 

(ii) Here EA is the sum of the first and third rows. 

» E = [201000], 

» E*A 
ans = 

17 10 -21 4 -10 
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Here, EA is the twice the first row plus the third row. In general EA will be made up by 
multiplying the ith row of A by the zth element of E, and then adding these rows together. 

(d) In general, if E is zero except for a p in the fcth entry, then EA will be the kth row of A multi- 
plied by p. To test this, 

» A = round( 10*(2*rand(3,5)-l)) 

A = 

-5 -9 -4 -1 5 

-8 0 8 9 5 

9-2 1-9 7 

» E = [ 0 3 0]; '/, A 3 in the 2nd entry . 

» E*A '/. This will be the 2nd row of A, multiplied by 3. 
ans = 

-24 0 24 27 15 

This may be repeated with a different choice of E and A. 

(e) In general, if E has a p in the fcth entry and a q in the j th entry, and zeros elsewhere, EA will be 
the Ath row of A times p plus the yth row of A times q. 

» E = [0 3 1]; */. 3 in the 2nd entry . 1 in the 3rd entry . 

» E*A '/, This will be 3*(2nd row) + l*(3rd row), 
ans = 

-15 -2 25 18 22 

This may be repeated with a different choice of E and A. 

(f) You should find the same results as in (d) and (e), except using columns instead of rows: 

» F = [0; 0; 2; 0; 0] '/, A 2 in the 3rd entry. 

F = 

0 

0 

2 

0 

0 

» A*F '/, This will be 2 times the 3rd column of A. 
ans = 

-8 

16 

2 

This may be repeated with different A’s and F’s. 

9. (a) 

» A = round( 10*(2*rand(3)-l)) 

A = 

-7 7 5 

-10 3 10 

4 5 8 

» B = round( 10*(2*rand(3)-l)) 

B = 

-5 0 -2 

-4 2 7 

-3 7 -5 
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» UA = triu(A) '/ The upper triangular part of A. 

UA = 

-7 7 S 

0 3 10 

0 0 8 

» = triu(B) */, The upper triangular part of B. 

UB = 

-5 0 -2 

0 2 7 

0 0-5 

» UA*UB 
ans = 

35 14 38 

0 6 -29 

0 0 -40 

This has the property that UA*UB is also upper triangular. When this is repeated for larger matrices, 
UA*UB will still be upper triangular. 

(b) If A and B are upper triangular matrices, then C = AB is also an upper triangular matrix. 
Proof: The matrix A is upper triangular when a,ij = 0 for any i > j. Since B is also upper trian- 
gular, bij = 0 when i > j. What we wish to show is that c t j = 0 for i > j. Assume that i > j. If 
we use summation notation, 

n 

c i j — ^ ^ a ik^kj • 

kzzl 

Split this sum into two sums: k < i and k > i: 

i— 1 n 

Cij = ^ a ik b k j + ^2 a ikb k j. 
k=l k=i 

In the first sum, we have i > k so a,-* = 0. In the second sum, we have k > i > j , so b k j = 0. In 
either sum, we are just adding up 0, so Cij = 0. Which implies that C is upper triangular. 

(c) The product of two lower triangular matrices is also lower triangular. To test this, generate two 
random lower triangular matrices.: 

» A = tril( round( 10*(2*rand(3)-l)) ) 

A = 

-2 0 0 

0 7 0 

-7 2 -7 

» B = tril( round( 10*(2*rand(3)-l)) ) 

B = 

10 0 0 

-2 -5 0 

-7 0 -7 

» A*B '/, This is also loner triangular, 
ans = 

-20 0 0 

-14 -35 0 

-25 -10 49 
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10. (a) 



» A = 


round ( 


10*(2*rand(5) 


-D) 


A = 










-2 


-7 


-8 


-7 


6 


-6 


9 


-7 


6 


9 


-9 


-2 


-9 


-1 


3 


8 


-7 


-3 


-3 


-6 


-1 


8 


-5 


-1 


4 


» B = 


triu(A,l) 






B = 










0 


-7 


-8 


-7 


6 


0 


0 


-7 


6 


9 


0 


0 


0 


-1 


3 


0 


0 


0 


0 


-6 


0 


0 


0 


0 


0 



B is the matrix made from A with nonzero elements above the main diagonal. If we call the di- 
agonal above the main diagonal the first diagonal and the one above that the second, and so on, 
B has zeros below the first diagonal. Type help triu for more information about triu. 



» B'2 
ans = 

0 0 

0 0 

0 0 

0 0 

0 0 



49 -34 -45 

0 7 -57 

0 0 6 

0 0 0 

0 0 0 



(b) 



B 2 has zeros below the second diagonal. Similarly, B 3 will have zeros below the third diagonal. 
B 5 = 0, so B is nilpotent with index 5. 



» B = 


triu(A,2) 








B = 










0 


0 


-8 


-7 


6 


0 


0 


0 


6 


9 


0 


0 


0 


0 


3 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



This time, B has zeros below the second diagonal. 



» B'2 
ans = 

0 

0 

0 

0 

0 



0 0 
0 0 
0 0 
0 0 
0 0 



0 -24 

0 0 

0 0 

0 0 

0 0 



B 2 has zeros below the fourth diagonal. B 3 will be zero, so B is nilpotent with index 3. 

(c) If we repeat (a) with a 7 x 7 matrix, we will find that B k will have zeros below the kth diagonal. 
This means that B will be nilpotent with index 7. If we repeat (b) with a 7 x 7 matrix, we will 
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find that B k 


will have 


zeros below the 2Hh diagonal. This means that B will be nilpotent with 


index 4. 












(d) If we choose B = triu(A,j), then B k will have zeros below the j kth diagonal. For a 6x6 matrix, 


B will be nilpotent with index 3 when we choose j to be the smallest number so that j ■ 3 > 6, 


which is 2. 












» A = 


round ( 


10*(2*rand(6)-l) ) 




A = 












7 


-6 


-5 


-4 


6 


-7 


-8 


-8 


4 


-6 


-7 


10 


-8 


-2 


-4 


-10 


-2 


-5 


5 


9 


6 


-6 


2 


-5 


3 


9 


6 


10 


-6 


-8 


-6 


-2 


-2 


-5 


7 


-6 


» C = 


triu(A,2) 








c = 












0 


0 


-5 


-4 


6 


-7 


0 


0 


0 


-6 


-7 


10 


0 


0 


0 


0 


-2 


-5 


0 


0 


0 


0 


0 


-S 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


» C~2 


*/. This 


is not 


zero. 




ans = 












0 


0 


0 


0 


10 


45 


0 


0 


0 


0 


0 


30 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


» C~3 


'/. This is zero 


. So 


index o t nilpotency of C is 3. 


axis = 












0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


11. First generate the eight matrices: 








» A = round( 6*( 2*rand(2) - 


D) 




A = 












3 


■2 










6 


3 










» B = round( 6*( 2*rand(2) - 


l)) 




B = 












-5 


5 










2 


3 
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» C = round ( 6*( 2*rand(2) - 1)) 
C = 

-1 0 

3 -3 

» D = round( 6*( 2*rand(2) - 1)) 
D = 

-3 -4 

-2 0 

» E = round( 6*( 2*rand(2) - 1)) 
E = 

5 -5 

5 5 

» F = round ( 6*( 2*rand(2) - 1)) 
F = 

0 -2 

0 6 

» G = round( 6*( 2*rand(2) - 1)) 
G = 

0 -5 

-3 5 

» H = round( 6*( 2*rand(2) - 1)) 



H = 


-5 

0 

» AA = 


-1 

-3 

C A B; 


C D] 


!, The block matrix. 


AA = 


3 


-2 


-5 


5 


6 


-3 


2 


-3 


-1 


0 


-3 


-4 


3 


-3 


-2 


0 


» BB = 


C E F; 


G H] 


'/, Another block matrix. 


BB = 


5 


-5 


0 


-2 


5 


5 


0 


6 


0 


-5 


-5 


-1 


-3 


5 


0 


-3 


» AA*BB 
ans = 


-10 


25 


25 


-28 


24 


-70 


-10 


-23 


7 


0 


15 


17 


0 


-20 


10 


-22 


» K = 


[ A*E+B*G 


A*F+B*H ; C*E+D*G C*F+D*H] 


K = 


-10 


25 


25 


-28 


24 


-70 


-10 


-23 


7 


0 


15 


17 


0 


-20 


10 


-22 
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» AA*BB 


- K 






ans = 








0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


fact, AA*BB will j 


always be the same as K 



12 . 



» A = 


round ( 


10*(2*rand(3,4)-l)) 


A = 










8 


9 


5 


-10 




1 


-9 


7 


4 




-1 


5 


-7 


7 




» B = 


round ( 


10*(2*rand(4,5)-l)) 


B = 










7 


-1 


1 


0 


8 


-5 


-4 


6 


9 


2 


-2 


-6 


-9 


5 


7 


1 


-7 


1 


1 


-7 


» D = 


A( : , 1)*B(1 , : 


) + A( : 


,2)*B(2, 




A ( : ,3)*B(3, : 


) + A( : 


,4)*B(4, 


D = 










-9 


-4 


7 


96 


187 


42 


-35 


-112 


-42 


11 


-11 


-26 


99 


17 


-96 


» D - 


A*B 








ans = 










0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



In general, D will be AB. To see this, notice that the kth term in this sum is the product of the kth 
column of A and the ibth row of B. The ijth entry in this matrix will be: au-b^j . When we add to- 
gether all of the matrices, the ijth entry will be aikbkj , which is the ijth entry of AB. 

13. (a) 



» AB = zeros (3, 5); 

» AB(1 , [1 2]) = [1 1]; 

» AB(2, [2 3]) = [1 1]; 

» AB(3, [1 4 5]) = [1 1 1] 

AB = 

110 0 0 
0 110 0 
10 0 11 
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» BC = zeros(5,8); 

» BC(1,[1 3 5]) = [1 1 1]; 

» BC(2, [3 4 7]) = [1 1 1]; 

» BC(3, [1 5 6 8]) = [1 1 1 1]; 

» BC(4,8) = 1; 

» BC(5, [5 6 7]) = [1 1 1] 

BC = 

1010100 0 
0 0110010 
10001101 
00000001 
00001110 

» CD = zeros (8, 10); 

» CD(1 , [1 2 3]) = [1 1 1]; 

» CD(2, [3 4 6]) = [1 1 1] ; 

» CD (3, [8 9 10]) = [111]; 

» CD (4, [4 5 7]) = d 1 1]; 

» CD(5, [1 4 6 8]) = [1 1 1 1] ; 

» CD (6, [2 4]) = d 1]; 

» CD(7,[1 5 9]) = d 1 1]; 

» CD (8, [1 2 4 6 7 9 10]) = [1 1 1 1 1 1 1] 

CD = 

1110000000 

0011010000 

0000000111 

0001101000 

1001010100 

0101000000 

1000100010 

1101011011 

(b) Person i in group 1 will have contact through person j in group 3 through person k in group 2 
if ABik — 1 and BCkj = 1, or in other words ABikBCkj = 1. If we sum over k, this will tell 
us how many indirect contacts person i has with person j. So the indirect contact matrix will be 
AB ■ BC. Similarly, to get from group 1 to group 4 via groups 2 and 3, we will multiply all three 
matrices together: 

» AD = AB*BC*CD 
AD = 

3112211332 

4314222232 

5314131332 

None of the entries are zero. This signifies that everybody in group 1 has some indirect contact 
with each person in group 4. Since the (1,5) entry is 2, there are 2 different paths that connect- 
ing person 1 in group 1 with person 5 in group 4. Similarly, since the (2,1) entry is 4, there are 4 
different paths connecting person 2 in group 1 with person 1 in group 4. 

(c) In order to find the total number of contacts a person in group 4 has, we add the rows of the in- 
direct contact matrix. From problem 12, we see that a simple way to do this is to multiply by 
[111]. 

» ones (1,3) * AD 
ans = 

12 73 10 564896 
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The person with the most indirect contacts is person 1, with 12 contacts. Person 3 has 3 con- 
tacts, which is the least. Similarly, we can add the columns together to find out how many indi- 
rect contacts each person in group 1 has. 

» AD * ones (10,1) 
ans = 

19 

25 

26 

Person 3 has the most indirect contacts with people in group 4, and so is the most dangerous. 

14. (a) Column one means that of the households using product 1, after one month 80/switch to prod- 
uct 3. Column two means that of the households using product 2, after one month 75/switch to prod- 
uct 3. Column three means that of the households using product 3, after one month 90/switch to 
product 2. 

(b) Since x is the distribution of households using the products after 0 months, Px will be the distri- 
bution of households using the products after 1 month. Similarly, P k x will be how many house- 
holds use each product after k months. 

( c ) 



» 


x = [10000 


; 10000; 10000]; 


» 


P = [ .8 


.2 


.05 




.05 


.75 


.05 




.15 


.05 


.9 ] 


p = 


= 








0.8000 


0.2000 0.0500 




0.0500 


0 . 


7500 0.0500 




0.1500 


0 . 


0500 0.9000 



It may be convenient to round off your answers. One way to do this is to set the output format 
to bank. This rounds numbers to the nearest hundredth. 

» format bank 
» P'5 * x 
ans = 

10275.83 
5840.35 

13883.83 

» P'10* x 
ans = 

9477.30 

5141.24 

15381.46 

» P~15 * x 
ans = 

9142.60 

5023.74 

15833.66 

» P*20 * x 
ans = 

9038.77 

5003.99 

15957.24 
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» P'25 * x 
ans = 

9010.03 
5000.67 

15989.30 

» P'30 * x 
ans = 

9002.52 

5000.11 

15997.37 

» P*35 * x 
ans = 

9000.62 

5000.02 

15999.36 

» P*40 * x 
ans = 

9000.15 

5000.00 

15999.85 

» P'45 * x 
ans = 

9000.04 

5000.00 
15999.96 

» P*50 * x 
ans = 

9000.01 
5000.00 

15999.99 

As n gets larger, P n x tends to (900,500, 1600). This may be interpreted by saying that eventu- 
ally every month the same number of households switch to product i, as switch from product i, 
for i = 1, 2, 3. 



» x = [0; 30000; 0]; 

» P'5 * x 
ans = 

12056.86 

9201.75 

8741.39 

» P'50 * x 
ans = 

9000.04 

5000.00 

15999.96 

Although for small n, P n x is different from that in (c), for large n, P n x. tends to the same value: 
(9000,5000,16000). 

(e) For any choice of x, P n x will tend to the same vector: P n x — >• (9000, 5000, 16000). 
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(0 

» P"50 
ans = 

0.30 0.30 0.30 

0.17 0.17 0.17 

0.53 0.53 0.53 

» 30000 * P'50 
axis = 

9000.00 9000.04 8999.99 

5000.00 5000.00 5000.00 

16000.00 15999.96 16000.01 

When n is large, the columns of 30000P" are very close to each other, and to the vector 
(9000, 5000, 15000). So that when 300000P" is multiplied by any vector x, it will be close to 
(xi + *2 + £ 3 ) times this vector. 

(g) 

» P = [.8 .1 .1; .05 .75 .1; .15 .15 .8]; 

» n = 50; '/, Try this for n= 5,10,15,25... 

» 1000 * P"n 
ans = 

333.33 333.33 333.33 

238.10 238.10 238.10 

428.57 428.57 428.57 

» format */, Don’t forget to return to normal output format at 
» '/, the end of this problem. 

In the limit, the columns of 1000P" tend to the vector (333.33,238.1,428.57). This will be the 
long term distribution of cars no matter what the starting distribution is. A car rental agency 
could use this information by planning to have a larger parking lot at office 3 than at office 1, or 
by planning to hire more mechanics at office 3 than at office 1. 

15. (a) Column 1 says that 40% of the fish in group 1 survive to belong to group 2 the next year. Col- 
umn 2 says that 20% of the fish in group 2 stay in group 2, and 50% of the fish in group 2 survive to 
be in group 3 the next year. This would happen if group 2 covers fish in an age range of more than 1 
year. Column 3 says that each fish in group 3 has 2 babies, and that 20% survive and stay in group 3, 
and 50% survive and enter group 4. Column 4 says that each fish in group 4 has 2 babies, and then 
20% survive to stay in group 4 and 40% survive to enter group 5. Column 5 says that each fish in 
group 5 has a 10% chance of survival to stay in group 5. 

(b) If x is the distribution of fish at this moment, then y = Sx will be the number of fish after one 
year. Since y is also a distribution of fish, Sy = S • Sx = S 2 x will be the number of fish one year 
later, or two years from now. 

(c) 

»S=[0 0 2 2 0 

.4 .2 0 0 0 

0 .5 .2 0 0 

0 0 .5 .2 0 

0 0 0 .4 .1 ]; 
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» x = [5000; 10000; 20000; 20000; 5000]; 

» n = 10; floor (S"n*x) 
ans = 

41016 

21666 

12949 

7754 

3709 

» '/, repeat for n=20,30, . . . 

» n = 50; floor (S*n*x) 
ans = 

49063 

24412 

15183 

9443 

4179 

After about 10 years, the population grows steadily. The growth rate will be exponential. 

(d) 



» S(1 ,3) = 1; 

» n = 10; floor(S~n*x) 
ans = 

15774 

8573 

5572 

4100 

2105 

» '/, repeat for n=20,30, ... 

» n = 50; floor (S~n*x) 
ans = 

550 

305 

212 

147 

71 

This time, the population decays exponentially. Not enough new fish are being born to keep the 
population steady. 

» S(l,3) = 2; S(3 ,2) = .3; 

» n = 10; floor (S~n*x) 
ans = 

13280 

8007 

3334 

2438 

1321 

» '/, repeat for n=20,30, . . . 

» n = 50; floor (S“n*x) 

100 
58 
25 
18 
9 



ans 
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Again, the total population decays when the survival rate from group 2 into group 3 is dropped 
from 50% down to 30%. Not enough fish survive to be group 3 in order to create new fish. 

(e) At the end of the first year, the number of fish that have survived will be Sx. If we harvest h 
of these fish, we will end up with u = Sx — h fish. Similarly, Su will be the number of these 
remaining fish that survive to the end of the second year. If we harvest h of these, we will have 
Su — h fish remaining. 



» S(1 ,3) = 2; S(3,2) = .5; 

» h = [0 ; 0 ; 0 ; 0 ; 2000] ; 

» u = S*x - h 
u = 

80000 

4000 

9000 

14000 

6500 

» u = S*u - h '/. This command should be repeated several times, 
u = 

46000 

32800 

3800 

7300 

4250 

After two more iterations of the last command, the vector u will have a negative number in the 
last entry. This means that we may harvest this many fish for 3 years, and at the end of the 
fourth year, there will be less than 2000 mature fish to be harvested. 

(g) After repeating the following experiment, for different values of n, you will find that u begins to 
drop for the first few years, and then begins to increase after about 5 years. If n is chosen to be 
1530 or smaller, u will never be negative. 

» n = 1530; 

» h = [0 ; 0 ; 0 ; 0 ; n] ; 

S*x -h 

80000 
4000 
9000 
14000 
6970 

S*u -h '/, This step should be repeated several times, until 
'/, u begins to rise again. 

46000 
32800 
3800 
7300 
4767 

(h) If the above experiment is repeated using nonzero values in the last two entries of h, the total 
harvest can be improved. For example, if h = (0,0,0, 1500,790), the total harvest will be 2290, 
and the population of fish will not become negative over 15 years. These values can be improved 
by slightly changing h. 



» u = 
u = 



» u = 
» 
u = 




78 Chapter 1 Systems of Linear Equations and Matrices 



Instructor’s Manual 



Section 1.7 



‘■CD ft)-© 



3 6-7 
2-1 3 



/ 0 1 -1 
5. 1 0 1 

V 3 2 0 




x\ + x 2 - x 3 = 7 
7. 4xi — X 2 + 5x 3 = 4 
6xi + x 2 + 3x 3 = 20 

2xi + 3*2 + x 3 = 2 
4x 2 + x 3 = 3 

6xi + 2 x 2 + x 3 = 2 
13. — 2xi -f 3x2 + x 3 = 4 
0 = 2 

/2 0 0\ /xA 

16. = I 0 4 0 , x = I x 2 , 

\0 0 — 5/ \x 3 J 

n ( 1 - 3 2 W* - 3 2N ) 

u - \-2 6 - a ) V ° 0 °7 

X = (2,0) + x 2 (3, 1) 



1 -1 3\ /xi 

4 1 -1 x 2 

2-1 3/ \x 3> 

4 -1 1 - 1 \ / 

3 1-5 6 

2-1 1 0 / . 



2 3-1 
-4 2 1 

7 3-9 





X2 = 2 
xi =3 

xi =2 

x 2 =3 
x 3 = -5 

X4 = 6 

3xi + x 2 + 5 x 3 = 6 
2xi + 3x2 + 2 x 3 = 4 



2xi + x 3 = 2 

9. — 3xi + 4x2 = 3 

5x 2 + 5x 3 = 5 

2xi + 3x2 + x 3 = 0 
12. 4xi - x 2 + 5 x 3 = 0 
3xi + 6x 2 — 7 x 3 = 0 

7xi + 2 x 2 = 1 
15. 3xi + x 2 = 2 
6xi + 9x2 = 3 



5 , xx = 3/2, x 2 = 5/4, x 3 = —2/5 

2 / 



x p — (2, 0); f j 6 o 



1 -3 0 
0 0 0 



,x h = x 2 (3, 1), 



1-11 6 
3-33 18 / 

1-11 0\ 
3 -3 3 0 / 



A -1 1 6 
A 0 0 0 

1-11 0> 
0 0 0 0 ; 



x = (6, 0, 0) + (x 2 — x 3 , x 2 , x 3 ) 



x p = (6,0,0) 



Xh = (x 2 - x 3 ,x 2 ,x 3 ) 



1-1-1 2 
2 12 4 

1-4-5 2 

1-1-1 0 
2 12 0 
1-4-5 0 



1-1-1 2 
0 3 4 0 

0- 3-4 0 

1- 1-1 0 
0 3 4 0 
0-3-4 0 



1 0 1/3 2' 
014/3 0 
0 0 0 0 , 

1 0 1/3 O' 
014/3 0 
0 0 0 0 , 



Xp = (2,0,0) 



x„ = x 3 (-l/3, -4/3,1) 



x = (2,0,0) + x 3 (-l/3, -4/3,1) 
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/I -1 -1 2\ /I -1 -1 2\ / 1 0 1/3 2\ 

20. 2 1 2 4 J — I 0 3 4 0 — 0 1 4/3 0 x p = (2,0,0) 

\1 -4 -5 2/ \° ' 3 - 4 0/ \0 0 0 0/ 

/I -1 -1 0\ /I -1 -1 0\ /I 0 1/3 0\ 

2 1 2 0 — 0 3 4 0 — 0 1 4/3 0 x h = a? 3 (— 1/3, -4/3, 1) 

\ 1 -4-5 0/ \0 -3 -4 0 / \0 0 0 0/ 

x = (2,0,0) + * 3 (-l/3, -4/3,1) 



A 1 -1 2 f 11-12 3\ A 0 3 -5 -A 

V 3 2 1-1 5/ VO -1 4-7 -4/ V° 1 -4 7 V 

A 1 01 2 3\ A 1 -1 2 3\ A 0 3 -5 -l\ 

\3 2 1 -1 5 ) VO -1 4 -7 -4/ ~ + V0 1 -4 7 4/ 

x* = (-3x 3 + 5x 4 , 4 x 3 - 7x4, x 3 , x 4 ) 
x = (-1,4, 0,0) + (— 3 x 3 -(- 5x 4 , 4x 3 - 7x 4 , x 3 , x 4 ) 



x p = (-1,4, 0,0) 



/ 1-1 1-1 -2\ /l -1 1 -1 -2\ /I 0 2-1 -1\ /I 0 0 -1/2 5\ 

22. -2 3 -1 2 5 — 0 1 1 0 1 — 0 1 1 0 1 — 0 1 0 1/4 4 

V 4 -2 2 -3 6/ \° 2 2 1 14 / \0 0-4 1 12/ \0 0 1 -1/4 -3/ 

x p = (5,4, -3,0) 

/ 1 -1 1-1 0\ /I -1 1-1 0\ /I 0 2 -1 0\ /I 0 0 -1/2 0\ 

-2 3 -1 2 0 — 0 11 0 0 — 01 1 0 0 — 010 1/4 0 

V 4 -2 2 -3 0/ \0 2 2 1 0 / \0 0 -4 1 0/ \0 0 1 -1/4 0/ 

x k = x 4 (l/2, -1/4, 1/4,1) 
x = (5, 4, —3, 0) + x(l/2, —1/4, 1/4, 1) 



23. Plugging y = c^yi + c 2 j /2 into the left side of the differential equation gives, since ( cy )" - cy" and 
(yi + Vi)" = y" + y'l 

Ciy” + C 2V2 + a ( x )( c iyl + c 2 y' 2 ) + 6( x )( c iyi + c 2 y 2 ) =-ci(yi' + a(x)y[ + b(x)yi) + c 2 (y'^ + a(x)y' 2 + b(x)y 2 ) 

— c i(0) + c 2 (0) = 0 



24. y" - y” + a(x)(j/p - y' ) + 6(x)(y p - y q ) 

= ( Vp + a(x)x/ p + b(x)y p ) - (y" + a(x)j/ s + b{x)y q ) 
= /( x ) - /( x ) = 0 

Thus, y p (x) - y q (x) solves y"(x) + a(x)j/(x) -f 6(x)y(x) = 0. 
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l. (a) 

» A = round(10*( 2*rand(3)-l) ) 

A = 

-6 4 0 

-9 9 7 

4 -2 -9 

» b = round(10*(2*rand(3, 1)-1) ) 
b = 

-9 
1 
3 

» R = rref ( [A b] ) 

R = 

1.0000 0 

0 1.0000 

0 0 

» x = R(: .4) 
x = 

3.8053 
3.4579 
0.5895 

» A*x '/. First find A*x 
ans = 

-9.0000 

1.0000 

3.0000 

» A*x - b / Compare Ax with b. 
ans = 

1.0e-14 * 

0 

-0.2665 

0.0888 

In theory, Ax — b should be zero, but in practice, the computer will have some round-off error. 
Here the error is on the order of 10~ 14 . 

» y = x(l)*A( : , 1) + x(2)*A( : ,2) + x(3)*A(:,3) 

y = 

-9.0000 

1.0000 

3.0000 

» y-b 
ans = 

1.0e-14 * 

0 

-0.2665 

0.0888 

Again, y is the same as A * x, so it will be b up to some round-off error. 



0 3.8053 

0 3.4579 

1.0000 0.5895 
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(b) (i) 

» A = [4 9 17 5 

2 15-1 
5 9 19 4 

9 5 23 -4 ] 

» b = [11; 9; 16; 40] ; 

» R = rref ( [A b] ) 

R = 

10 2-15 

0111-1 
0 0 0 0 0 

0 0 0 0 0 

The general solution will have x 3 and x 4 arbitrary, x x = — 2ar 3 + X 4 + 5, and ar 2 = — x 3 — x 4 — 1. 
We may pick, for example, X 3 = 1 and x 4 = 0, to get the particular solution: 

» x = [3 ; - 2 ; 1 ; 0 ] ; 

(ii) 

» A*x */. This should be the same as b. 
ans = 

11 

9 

16 

40 

» y = x(l)*A(: ,1) + x(2)*A( : ,2) + x(3)*A(:,3) + x(4)*A( : ,4) 

y = 

11 

9 

16 

40 

As in (a), if x is a solution of the system with [A b] as augmented matrix, Ax = b, then Ax 
and x x A(\, 1) H (- x 4 A(:, 4) = y are both b. 

(iii) If this is repeated with other choices of X 3 and x 4 , the same results will occur: Ax = b and 
y = b. For example: 

» x = [ 4; -3; 1; 1]; 

» A*x 
ans = 

11 

9 

16 

40 

(iv) 

» y = x(l)*A(:,l) + x(2)*A( : ,2) + x(3)*A(:,3) + x(4)*A(:,4) 

y = 

11 

9 

16 

40 

(c) The solution of a system of equations represented by [A b] is the same as the solution of the ma- 
trix equation Ax = b. Also, multiplying a matrix by a single column is equivalent to adding 
multiples of the columns of this matrix. 
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2. (a) The ij entry of Ax is the inner product of the ith row of A with the jth column of x. Since x has 
only one column, j must always be 1, and the first column of x is just x itself. Since Ax = 0, the 
inner product of the ith row of A with x is zero. This means that the ith row is orthogonal to x. 

(b) This can be done by solving Ax = 0 where A is the matrix whose rows are the given vectors. 

» A = [1 2 -3 0 4; 5-5201]; 

» rref(A) 
ans = 

1.0000 0 -0.7333 0 1.4667 

0 1.0000 -1.1333 0 1.2667 

For a solution, we may choose x 3 , x 4 , and x$ arbitrarily, and then set x\ = . 7333 x 3 — 1.4667x5 
and x 2 = 1.1333x3— 1.2667x5. 

3. (a) The matrix x solves the nonhomogeneous system Ax = b. The matrix z was a solution of the 
homogeneous system Az = 0. If we set y = x + sz, then we found that y was also a solution of the 
nonhomogeneous system Ay = b. The corollary tells us the converse, i.e. that any such solution can 
be written as y = x + sh where H is a solution of Ah = 0. 

(b) (i) Refer to the solution of 1(b) above. The matrix R is the reduced echelon form of [A b]. Since 
two variables can be chosen arbitrarily in the solution of Ax = b, there are infinitely many solu- 
tions. 

(ii) 

» x = A\b 

Warning: Matrix is singular to working precision, 
x = 

Inf 

Inf 

Inf 

Inf 

» '/, Since A is singular, it has no inverse. Instead, enter a solution 
» '/, from the answer to 1(b). 

» x = [ 5; -1; 0; 0] ; 

(hi) 



» rref(A) 
ans = 

10 2-1 
0 111 
0 0 0 0 

0 0 0 0 

Solutions of Ax = 0 are of the form X 3 and X 4 are arbitrary, xi = — 2x 3 + x 4 , and x 2 = 
— X 3 — x 4 . 

Pick one solution, with x 3 = 1, x 4 = 0 

» z = [-2; -1; 1; 0] ; 

» A*(x+z) '/. This should yield b. 

11 
9 

16 
40 



ans 
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Another solution, with £3 = 0 , £4 = 1 . 

» z = [1; -1; 0; 1]; 

» A*(x+z) Again, this should be b. 
ans = 

11 

9 

16 

40 

This can be repeated two more times by choosing other values for X3 and X4. 

4 . (a) 



» A = [5580 
4 5 8 7 

3 9 8 9 

9 116]; 

» rref(A) 
ans = 

10 0 0 
0 10 0 

0 0 10 

0 0 0 1 

If we were to reduce [ A b] for any b, we would still have the same left hand side. Since there 
is a pivot in every row, the system will be consistent, and there will be a unique solution. 

(b) If we solve Ax = b, for x, then we may write b as 

b = CjXj + C2Z2 + C3X3 + 040:4 

where c , is the ith column of A. For example: 

» b = round( 10*(2*rand(4, 1 ) — 1 ) ) 
b = 

-6 
-9 

4 
4 

» R = rref ( [A b] ) 

R = 

1.0000 0 

0 1.0000 
0 0 

0 0 

» x = R(: ,5) 

0 . 6476 
3.5799 
-3.3922 
-0.3361 



0 


0 


0.6476 


0 


0 


3.5799 


1.0000 


0 


-3.3922 


0 


1.0000 


-0.3361 



x 
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» */, Now we may write b as a combination of columns of A: 

» y = A(:,l)*x(l) + A(:,2)*x(2) + A(:,3)*x(3) +A(:,4)*x(4) 

y = 

- 6.0000 

-9.0001 

3.9999 

3.9999 

Up to round off error, y is the same as b. This can be repeated two more times. 



» A = [ 5 5 -50 

4 5 -6 7 

3 9 -15 9 
9 1 7 63; 

» rref(A) 
ans = 

10 10 
0 1-20 
0 0 0 1 

0 0 0 0 

It is possible that there is a b for which the reduced form of [A b] does not have a zero in the 
(4,5) entry. Since the left hand side has all zeros in the fourth row, this would mean the system 
is inconsistent. By experimenting with several possible b you can find that if b = (1,0, 0, 0 )*, 
there will be no solution: 

» b = [1; 0; 0; 0] ; 

» rref ( [A b] ) 
ans = 

10 10 0 
0 1-200 
0 0 0 1 0 

0 0 0 0 1 

Notice that the above solution is inconsistent. 

(d) If you start with a vector b which is a combination of columns of A, then there will always be a 
solution of .Ax = b. 

» k = round( 10*(2*rand(4, 1 ) — 1 ) ) */. generate 4 random numbers, 
k = 

9 

-2 

0 

7 

» '/. Write b as a combination of columns of A: 

» b = A(: , 1 ) *k ( 1 ) + A(:,2)*k(2) + A(:,3)*k(3) +A(:,4)*k(4) 
b = 

35.0000 

23.6475 

5.9754 

76.9836 
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» rref([ A b] ) 








ans = 

1.0000 


0 


1 . 0000 


0 


9.0000 


0 


1.0000 


-2.0000 


0 


-2.0000 


0 


0 


0 


1.0000 


-0.3361 


0 


0 


0 


0 


0 



This system is consistent. 

(e) To see that this system will always have a solution, we only need to show that it has at least one 
solution. However, writing b as a combination of columns of A using the scalars is equivalent 
to the matrix multiplication .Ak = b. This means that the vector k will be a solution of ylx = b. 
Since the system has a solution, it is consistent. 

5. (a) (i) For A from 4(c): 

» A = C 5 5 -50 

4 5 -6 7 

3 9 -15 9 

9 1 7 6]; 

» rref(A) 
ans = 

10 10 
0 1-20 
0 0 0 1 

0 0 0 0 

(ii) The solutions of this homogeneous system have x 3 arbitrary, and x\ = — x 3 , x 2 = 2x 3 , and 
x 4 = 0. 

(iii) If we set x 3 = 1, then we have x\ = —1 and x 2 = 2. This corresponds to 

0 = ylx = — lei + 2 c 2 + lc 3 + 0c 4 , 
where c,- is the ith column of A. Which can be rewritten as 

c 3 = lei - 2 c 2 . 

To check this: 

» 1*A(:,1) - 2*A(:,2) '/, This should be the same as A(:,3). 

ans = 

-5 

-6 

-15 

7 

» A(:,3) '/, This is the third column, 

ans = 

-5 

-6 

-15 

7 

(iv) This system only allows one arbitrary variable. 

(v) If we let x be the third column of rref (A) then Ax will be the third column of A. 
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(b) 



» A = [4 9 17 5 

2 1 5-1 

5 9 19 4 

9 5 23 -4 ]; 

» rref(A) 
ans = 

10 2-1 
0 111 
0 0 0 0 

0 0 0 0 

The solution of this system will have r 3 and X4 arbitrary, and x\ — — 2X3 + X4, and r 2 — —X3 — X4. 
If X 3 = 1 and X 4 = 0, then we have = —2 and X 2 = — 1. This corresponds to 



c 3 = 2cj + lc 2 . 



Similarly, if x 3 = 0 and X 4 = 1, then we have ajj = 1 and z 2 = — 1. This corresponds to 



c 4 = — lcj + lc 2 . 



As above, if x is the third column of rref (A) then Ax will be the third column of A. Similarly, 
if x is fourth column of rref (A) then Ax will be the fourth column of A. To check this: 

» R = rref (A); 

» x = R(; ,3) 
x = 

2 

1 

0 

0 

» A*x */, This will be the 3rd column of A. 
ans = 

17 

5 

19 

23 

» x = R(: ,4) 
x = 

-1 

1 

0 

0 

» A*x '/, This is the 4th column of A. 
ans = 

5 

-1 

4 

-4 
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(c) First we generate a random matrix and modify it as directed. 

» A = round( 10*(2*rand(6,6)-l)) ; 

» A( : ,3) = 2*A( : , 2) - 3*A(:,1); 



V 

V 

> 


,5) 


= -A( 


:.l) 


+ 2*A( 


:,2) - 


3*A( 


» A( : , 
A = 


,6) 


= A( : 


,2) + 


4*A( : 


,4) 




10 




8 


-14 


-9 


33 


-28 


4 




-5 


-22 


8 


-38 


27 


5 




-1 


-17 


0 


-7 


-1 


3 




5 


1 


0 


7 


5 


-9 




0 


27 


-4 


21 


-16 


3 




-5 


-19 


10 


-43 


35 


» R = 


rref(A) 










R = 














1 




0 


-3 


0 


-1 


0 


0 




1 


2 


0 


2 


1 


0 




0 


0 


1 


-3 


4 


0 




0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


0 



This will be the same, for almost every random matrix that we start with. The solution of this 
system has £3, X5, and xq arbitrary and xi = 3x 3 -fx5, X2 = — 2x 3 — 2 x 5 — X6, and x 4 = 3x 5 — 4 x 6 . 
As in (b) above, if we set one of the three arbitrary variables to 1 and the others to 0, and write 
out x we will get: 



(~*\ 




/-!\ 




/°\ 


2 




2 




1 


0 




0 




0 


0 


, X — 


-3 


, or x — 


4 


0 




0 




0 


\ 0/ 




\ 0/ 




W 



Writing out Ax = 0 as a linear combination of columns of A, we will get the original equations: 

» A( : ,3) = 2*A(:,2) - 3*A(:,1); 

» A( : ,5) = -A( : , 1) + 2*A(:,2) - 3*A(:,4); 

» A(: ,6) = A(: ,2) + 4*A(: ,4) 
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1. Since det A — 2 ■ 2 — 3 • 1 = 1 ^ 0, .A 1 exists. A 



'- K-ry 



1 


1 


i (2 


1 


1 0\ (2 


1 


1 


0 


2 


o l ; 


i-U 


1 

2 


=T 1 / U 


1 


-3 


2 



t -? ) - ( 


1 0 


2 -1 \ 


. So A- 1 = ( \ 


-3 2) V 


0 1 


-3 2) 


\-3 2) 



2. Since det A = (—1) • (—12) — 1 • 12 = 0, A is not invertible. 



3. det A 



=o-i=-i#o ^-■=(-i)(_;-j) = (;;) 



4. Since det A = 3 — 3 = 0, A 1 does not exist. 



5. A is not invertible since det A = ab — ba = 0. 



Ill 100 
0 2 3 0 1 0 

5 5 1 0 0 1 



1 1 1 
0 1 3/2 
0 0-4 



1 0 0 13/8 -1/2 -1/8 

0 1 0 -15/8 1/2 3/8 

001 5/4 0 -1/4 



1 0 0 
0 1/2 0 
-5 0 1 



1 0 - 1/2 
0 1 3/2 
0 0 1 



/ 13/8 -1/2 -1/8 
A~ l = -15/8 1/2 3/8 

\ 5/4 0 -1/4 



1 - 1/2 0 
0 1/2 0 
5/4 0 -1/4 



3 2 1 10 0 

0 2 2 0 1 0 

0 0 -1 0 0 1 



1 2/3 1/3 1/3 0 0 

011 01/20 

0 0 1 0 0 1 



100 1/3 -1/3 -1/3 

010 01/2 1 
0 0 1 0 0 -1 



A -1 = 



1 0 -1/3 1/3 -1/3 0 

01 1 01/20 
0 0 1 0 0 -1 



1/3 -1/3 -1/3 
0 1/2 1 
0 0-1 



/111 100 
8. 0 1 1 0 10 
VO 0 1 001 



10 0 1-10 
0 10 0 1-1 

0 0 1 0 0 1 



/ 1 — 1 0 
A~ l = j 0 1-1 

Vo o i 



/1 6 2 10 0 

9. -2 3 5 0 1 0 

\ 7 12 -4 0 0 1 

A is not invertible. 



1 6 2 10 0 

0 15 9 2 1 0 

0 -30 -18 -7 0 1 



1 6 2 

0 1 0.6 

0 0 0 



1 0 0 
0.1333 0.0667 0 
-3 2 1 



3 10 10 0 



1 -1 2 

1 1 1 



0 1 0 
0 0 1 



1-12 0 10 

0 4-6 1-3 0 

02-1 0-11 



100 3/8 1/8 -1/4 

0 1 0 -1/8 -3/8 3/4 

0 0 1 -1/4 1/4 1/2 



10 1/2 I 1/4 1/4 0 



—01 -3/2 

\0 0 2 

3/8 1/8 -1/4 

-1/8 -3/8 3/4 
-1/4 1/4 1/2 



1/4 -3/4 0 
- 1/2 1/2 1 
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/ 2 -1 


4 


1 


0 


°\ 


11 . -1 0 


5 


0 


1 


° 


\ 19 -7 


3 


0 


0 


1 / 







A is not invertible. 




12. 










13 . 



/! 


1 


1 


1 


1 


0 


0 


°\ 






1 


1 


1 


1 


0 


0 


°V 




Z 1 


0 


3 


0 


2 


-1 


0 


°\ 


1 


2 


-1 


2 


0 


1 


0 


0 




0 


1 


-2 


1 


-1 


1 


0 


0 




0 


1 


-2 


1 


-1 


1 


0 


0 


1 


-1 


2 


1 


0 


0 


1 


0 




0 


-2 


1 


0 


-1 


0 


1 


0 




0 


0 


-3 


2 


-3 


2 


1 


0 


Vi 


3 


3 


2 


0 


0 


0 


1 / 




Vo 


2 


2 


1 


-1 


0 


0 


lJ 




Vo 


0 


6 


-1 


1 


-2 


0 


1 / 



'o 



0 0 2 
1 0 - 1/3 
0 0 1 - 2/3 

Vo 0 0 3 



-1 

1 

1 

-5 



1 



1 0 



- 1/3 - 2/3 0 
- 2/3 - 1/3 0 
2 2 1 




A~ l = 



/ 7 / 3 - 1 / 3 - 1 / 3 - 2 / 3 ' 
4/9 - 1/9 - 4/9 1/9 

- 1/9 - 2/9 1/9 2/9 

V — 5/3 2/3 2/3 1 / 3 / 



10 0 0 
0 10 0 
0 0 10 
0 0 0 1 



7/3 - 1/3 - 1/3 - 2/3 \ 
4/9 - 1/9 - 4/9 1/9 
- 1/9 - 2/9 1/9 2/9 
- 5/3 2/3 2/3 1 / 3 / 



14 . 



A is not invertible. 



/ 1 


0 


2 


3 


1 


0 


0 


°\ 




Z 1 


0 


2 


3 


1 


0 


0 








0 


2 


3 


1 


0 


0 


°v 


-1 


1 


0 


4 


0 


1 


0 


0 




0 


1 


2 


7 


1 


1 


0 


0 




' 0 


1 


2 


7 


1 


1 


0 


0 


2 


1 


-1 


3 


0 


0 


1 


0 




0 


1 


-5 


-3 


-2 


0 


1 


0 




1 0 


0 


-7 


-10 


-3 


-1 


1 


0 


V-l 


0 


5 


7 


0 


0 


0 


l) 




Vo 


0 


7 


10 


1 


0 


0 


l) 




Vo 


0 


0 


0 


-2 


-1 


1 


lJ 



15 . 



( 1 


-3 


0 


2 


1 


0 


0 


°\ 




/I 


-3 


0 


-2 


1 


0 


0 


°\ 


3 


-12 


-2 


-6 


0 


1 


0 


0 




0 


-3 


-2 


0 


-3 


1 


0 


0 


-2 


10 


2 


5 


0 


0 


1 


0 




0 


4 


2 


1 


2 


0 


1 


0 


V-l 


6 


1 


3 


0 


0 


0 


\) 




Vo 


3 


1 


1 


1 


0 


0 


1 / 



/! 


0 


3/2 


- 5/4 


5/2 


0 


3/4 


°\ 


/I 


0 


0 


1 


-2 


3 


3 


°\ 


0 


1 


1/2 


1/4 


1/2 


0 


1/4 


° - 


0 


1 


0 


1 


-1 


1 


1 


0 


0 


0 


- 1/2 


3/4 


- 3/2 


1 


3/4 


° 


0 


0 


1 


- 3/2 


3 


-2 


- 3/2 


0 


Vo 


0 


- 1/2 


1/4 


- 1/2 


0 


- 3/4 


1 / 


Vo 


0 


0 


- 1/2 


1 


-1 


- 3/2 


1 / 



yr 1 = 



/! 


0 


0 


0 


0 


1 


0 


2 \ 


0 


1 


0 


0 


1 


-1 


-2 


2 


0 


0 


1 


0 


0 


1 


3 


-3 


Vo 


0 


0 


1 


-2 


2 


3 


- 2 / 


/ 


0 


1 


0 


2 ' 


\ 




_ 




1 


-1 


-2 


2 










0 


1 


3 


-3 










-2 


2 


3 


-2 


/ 
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16. ABCC-'B-'A- 1 = ABIB~ 1 A~ 1 = ABB~ l A~ x = AIA~ X = AA~ X = I. Hence, by theorem 8, ABC 
is invertible and ( ABC)~ l = A~ 1 B~ 1 C~ 1 . 

17. Show A\A 2 ■ ■ ■ A m A^ ■ ■ ■ A^ 1 A \ 1 = I. Then by theorem 8, A\A 2 ■ ■ -A m is invertible with inverse 
,4-1 . . . 4 _1 4-1 

A m ^2 -'m • 



18. 



3 4 
-2 -3 



3 4 
-2 -3 



f 9 — 812 — 12\ 
V— 6 + 6 -8 + 9; 



19. 



If A = ±1 then A 2 = I. If an = —022 and 021012 = 1 — a? 1 then ( 022 ° 12 ^ = 

V fl21 «22 J 



J 22 



+ ai 2 02i — 022012 + 012022 \ _ j 



— 022021 + <J22 a 21 021 a 12 + a 22 ) 



( 4/5 -1/5 0\ 

20. I-A= -2/5 3/5 -3/5 

\ —1/5 —1/10 3/5 / 

/ 96.4\ 

x = (I-A)- 1 e& 235.7 

V 138.1/ 



/ 1.7857 0.7143 0.7143 \ 
(/-i4) _1 = 2.1492 2.8571 2.8571 

\ 0.9524 0.7143 2.3810/ 



21. Bx = 0 gives m equations in n unknowns. Since m < n, there are an infinite number of solutions. In 
particular, there exists a nonzero solution. Hence, there exists a nonzero vector x such that ABx = 0. 
By theorem 6, AB is not invertible. 



22. (a) det A — —i 2 — 2 = — 1 




(b) 



(c) 



det .4 = l-i 2 =2 




( 



0 

1 — i 



((1 + 0/2 0 \ 

V 0(1 - 0 / 2 J 



/l 


i 


0 


1 


0 


°\ 


( _i 


0 


1 


0 


1 


0 


\0 


1+i 


l-i 


0 


0 


1/ 



i 0 


1 


0 


°\ 


(l 


0 


i 


0 


i 


1 -1 


-i 


-1 


0 - i 


0 


1 


-1 


-i 


-1 


1+i l-i 


0 


0 


1/ 


\0 


0 


2 


-1+i 


1+i 



1 0 0 
0 1 0 
0 0 1 



(1+0/2 (1+0/2 -i/2 \ 
-(l+i)/2 (-1+0/2 1/2 
(-1+0/2 (1+0/2 1/2/ 



/ (1 + 0/2 (1 + 0/2 -i/2 \ 
A~ x = -(1 + 0/2 (-1 + 0/2 1/2 

V (- 1 + 0/2 (1 + 0/2 1 / 2 / 



0 

0 

1 



( sin 6 cos 9 0 \ 2 

cosd — sin 0 0 I = I ; this matrix is its own inverse. You can discover this by trying to do row 

0 01 / 

elimination or by finding the inverse of the upper left 2x2 block using (12) in Theorem 4. 

/ 1/2 0 0 \ 

24. A~ l — I 0 1/3 0 , from (A | /) -♦ (/ | +" 1 )- 

\ 0 0 1/4/ 
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25. Let D be a diagonal matrix. Suppose D is invertible with inverse A. Since AD = J, then aada = 1 
for each i. Hence, the diagonal components of D are nonzero. Conversely, suppose da / 0 for each i. 
Then the only solution to Dx = 0 is the trivial solution. By theorem 6, D is invertible. Or you can 
write down D~ x directly as in Problem 26. 



26. A- 1 = 



hi 
0 a 



0 

-l 

22 



V 0 0 



0 \ 

0 ' 

:i/ 



, from {A | /)-*(/ | A -1 ) 



/ 2 


1 


-1 


1 


0 


°\ 


0 


3 


4 


0 


1 


0 


\0 


0 


5 


0 


0 


1 / 



/I 0 -7/6 
0 1 4/3 
\0 0 5 



1/2 - 1/6 0 \ 
0 1/3 0 
0 0 1 / 



/ 1/2 —1/6 7/30 \ 

A- 1 = 0 1/3 -4/15 

\ 0 0 1/5/ 



/I 0 0 
0 10 
\0 0 1 



1/2 -1/6 7/30 \ 
0 1/3 -4/15 
0 0 1/5/ 



f 1 


0 


0 


1 


0 


°\ 


-2 


0 


0 


0 


1 


0 


V 4 


6 


1 


0 


0 


1 / 



(l 


0 


0 


1 


0 


°\ 




4 


6 


1 


0 


0 


1 


A does not have an inverse. 


\0 


0 


0 


2 


1 


0/ 





29. If U is upper triangular with all diagonals nonzero, then divide each row by its diagonal entry. The 
result is an echelon form of U with n pivots. So by Theorem 6.v, U is invertible. Conversely, if U has 
some diagonal zero, let j be such that ujj — 0 is the first zero on the diagonal. Then form x = (x*), 
with x k = 0 for k > j and Xj — 1, but the x k ,k < j unknown. Then the nonzero equations Ux — 0 
formed from the first j — 1 rows of U and this x can be backsolved to get the remaining x k , since they 
have the form 

akkXk H h a kj = 0, a kk ^ 0, k < j. 

Thus Ux — 0 has a nonzero solution and so U is not invertible. 



30. Row reducing ( U \ I) to (7 | U~ l ) requires only dividing through by the diagonals of U and then 
adding multiples of lower rows to higher rows, i.e. only backsolving is needed as U is already essen- 
tially in echelon form. But both these types of row operations only change the elements on or above 
the diagonal of the right hand block. Thus when this reduction is done A = U~ 1 will be upper tri- 
angular. (You can also solve this using partitioned matrices. To get an idea look at the solution to 
Problem 49.) 



31. 




1 - 1/2 
0 0 




then Ax = 0. For example, 




is one such vector. 



32. 



/ 1-1 3 
0 4-2 
\ 2 -6 8 




/ 1 0 5/2 
0 1 - 1/2 
\0 0 0 





then Ax = 0. For example, 




is such a 



vector. 
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33. Let c be the number of chairs and t the number of tables produced each day. We have 812 = 96 labor 
hours per day in the machine shop. Hence, for the machine shop we must have yy ■ c + ■ t = 96. 

Similarly, yy • c + ^jy ■ t = 8 • 20 = 160 for the assembly and finishing division. Write this system 



as Ax = b where A = 



x = A~ 1 b = 



289 



130,560 



384 


240 \ 


17 


17 


480 


640 


17 


17 ) 


640 


240 


17 


17 


480 


384 


17 


17 



, x = ( , and b 



Vieo; 



Since detA = — ^ 0, then 



96 

160 



-il) 



Hence, 3 chairs and 2 tables can be produced 



each day. 



34. Let / be the amount of love potion and c be the amount of cold remedy needed. The witch wants to 



find x = ( ' ! such that Ax = b where A = 
c 



V si ' 1 

13 13 

„ 2 10 
2 — 10 — 
13 13 



and b 



- Si 



Since detA = y^(40 



OOf) 1 9 

140 - 28 • 70) = # 0, then x = A~'b = — 



( 



\ 



140 _70 
13 13 

28 40 

"13 13 



f 10\ _ / 3/2 

V 14 / V 1 



Hence, 1- batches 



of love potion and 1 batch of cold remedy are needed. 



35. The farmer needs Ax = b where A — 008^ ’ X ” (fc)’ an< * ^ = 0 ) ^* nce A = — 0.01 ^ 

| . Thus, 4 units of type A and 5 units of type 



0 then x = A *b = —100 
B are needed. 



0.08 — 0.12 \ (l 

-o.i5 o.io y v i 



36. (a) 0.293 (b) 200,000 0.293 = 58,600 (c) 0 

(d) 50, 000 0.044 = 2,200 



/ 0.293 0 0\ 

37. (a) technology matrix A = I 0.014 0.207 0.017 1 

\0.044 0.010 0.216/ 

/ 0.707 0 0\ 

Leontief matrix = I — A = I —0.014 0.793 —0.017 I 

\ -0.044 -0.010 0.784/ 



/ 1.414 0 0\ / 13, 213 \ / 18, 689\ 

0.027 1.261 0.027 x = (/-A)~ 1 17,597 = 22,598 

\ 0.080 0.016 1.276/ \ 1,786/ \ 3,615/ 



(b) (/-A)-i = 
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It would require 18,689 pounds of agricultural products, 22,598 pounds of manufactured goods, and 
3,615 pounds of energy. 



38 . 



1 0 
0 1 



invertible 






invertible 



40. 



1 1 
3 3 



(ii) 



not invertible 




1 2/3 1/3' 

0 1 1 I invertible 

0 0 1 




1 6 2 \ / 1 6 2 ' 

0 15 9 J — >■ 10 1 3/5 | not invertible 

0-30-18/ V 0 0 0, 



43. 




- 1/2 2 \ 

1 —14 ] not invertible 

0 0 / 



44. 



f 1 


1 


1 A 




( ll 


1 


A 




( 11 


1 


A 


1 


2 - 


-1 2 1 




0 1 


-2 


1 




0 1 


-2 


1 


1 


-1 


2 1 i 


* 


0 0 


-3 


2 




0 0 


1 -2/3 


ll 


3 


3 2/ 




\0 0 


6 






loo 


0 


1/ 



invertible 



45. 



/ 1 0 2 3\ 


/ 1 0 2 3\ 


(-11 0 4 


(012 7 


2 1-13 


0 0 1 10/7 


\ — 1 0 5 7/ 


looo 0 / 



not invertible 



46. Since ana 22 — <Ji2 a 2i 7^ 0, then either an or a 12 is nonzero. We may assume without loss of generality 
that an ^ 0. 

f an “12 

\a2i <*22 

10 a 22 /det.A —021/661.4^ 
y0 1 — a 2 i/detyl an/det-Ay 



1 0 N / 1 ai 2 /«n 1/ «i 1 (A 

0 l) \0 a22 ~ ai2«2i/an — «2i/ a n 1/ 
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47. (i) Suppose A is invertible. Then Ax = 0 implies x = A~ 1 Q = 0. When reducing the augmented 
matrix (A|0) to reduced row echelon form, we must have (A|0) — > (J|0). Using the same elementary 
row operations will give A — ► I. Conversely, suppose A is row equivalent to I. Write (yl|/) and row 
reduce A to I. Then we will have (A|7) — ► ( I\B ). Hence, AB = I. We want to show BA = I. It will 
suffice to show BAx — x for every n-vector x. Note that B is row equivalent to I. Hence, for every x, 
we can find a y such that By = x. Thus, BAx = BA(By) = B(AB)y = Bly = By = x. Therefore, 
A is invertible. 

(ii) Suppose A is invertible. Suppose Ax = b = .Ay. Multiplying by A -1 we obtain A -1 (Ax) = 
A -1 (Ay). It follows that x = y. Conversely, suppose the system Ax = b has a unique solution for 
every n-vector b. This implies A is row equivalent to I. Therefore, by part (i), A is invertible. 

(iv) Suppose A is invertible. By (i) A is row equivalent to the identity matrix I n . I n is in row ech- 
elon form and has n-pivots. Conversely, suppose that an echelon form of A has n pivots. Then the 
reduced echelon form of A is the identity matrix /„ . Thus A is row equivalent to the identity matrix 
By (i) A is invertible. 



(^q^J ^ 12 ^ = (^q O'j we obtain the following equations: 



An + AA 2 i — I, A 12 + AA 2 2 — 0, A 2 i — 0 and A 22 = I. 
Solving for An and Ai 2 we get An = I and A 22 = —A. 






49. From ( P | ( o 11 q 12 J = ( L ® ) we obtain the following equations: AnBn = 7, An7?i 2 = 

V A 2 i -<4 22 J \ B 21 B 22 ) \U 1 J 

O , A 2 i5n + A 22 B 2 1 = O and A 2 ii?i 2 -f A 22 7? 22 = 7. Solving for the Bij we get 
Bn = AJ^ 1 , Bn = O, B 21 — A^i—AziA^) and 7? 22 = Aj/ . 

Thus / ^ 1 1 _ A j? °). 

\A 2 1 A 22 y V^22 ( — ^2lA u ) A 22 J 
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CALCULATOR SOLUTIONS 1.8 

The solution for Problem nn assumes the data has been entered into the matrix A18nn. 



50. The inverse of A1850, calculated by A1850 [2nd] FI [ENTER] is 

[[ .099858156028 -.076501182033 -.076312056738 ] 

[ -.101843971631 .272151300236 -.192056737589 ] . 
[ .143262411348 -.067139479905 .075177304965 ]] 

51. The inverse A1851, calculated by A1851 [2nd] FI [ENTER] is 

[[ -1.40754039497 .456014362657 .303411131059 ] 

[ .657091561939 -.166965888689 -.154398563734 ] . 

[ .439856373429 -.26750448833 -.03231597846 ]] 



52. The inverse of A1852, calculated by A1852 [2nd] FI [ENTER] is 



[[ -1.69701053654 1.60958317276 2.30794647641 

[ -.793941492883 1.67786019382 .951949895801 

t 1.95621257025 -.200637557471 -.462876267722 

[ -.654423076042 .641726876009 -.249248482807 



1.41376091667 

.261501231662 

.357416829406 

.536103553562 



] 

] 

] 

]] 



53. The inverse A1853, calculated by A1 8 53 [2nd] FI [ENTER] is 

[[ .03984485542 .00954069806 .035197499863 .010590297574 ] 

[ -.003683920834 .00460111628 4 . 2629849506E-4 -.005608453127 ] 

[ .018345358489 .009413418129 .008529325192 .00297300599 ] 

[ .019410170095 .007025671643 .025503332841 .015762697722 ]] 



54. The inverse of A1854, calculated by A1854 [2nd] [F] [ENTER] is 

[[ .333333333333 -.208333333333 1.675 -1.42142857143 ] 

[ 0 .125 -.325 .578571428571 ] 

[0 0 .2 -.114285714286 ] 

[0 0 0 -.142857142857 ] ] 

which has zeros below the diagonal. 



55. The inverse of A1855, calculated by A1855 [2nd] FI [ENTER] is 



.04329004329 




.125690401552 


-.196440007897 


.126871293496 


.203412258426 ] 


0 




.068965517241 


-.067111605488 


.035660968327 


.113274788336 ] 


0 


0 




-.02688172043 


.019100169779 


.007849426165 ] 


0 


0 




0 


.010964912281 


- . 003226349456 ] 


0 


0 




0 


0 


.02132196162 ] 



which has zeros below the diagonal. 



56. The results in Problem 54 and 55 suggest that the inverse of an upper triangular matrix is upper triangular. 
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MATLAB 1.8 



1. (a) (i) 



» A = [12 3; 254; 1 -1 10] ; 
» R = [A eye(3)] 

R = 



1 

2 

1 



2 

5 

-1 



3 

4 
10 



1 0 
0 1 
0 0 



0 

0 

1 



» rref(R) 
ans = 

1 0 

0 1 

0 0 



0 

0 

1 



» S = ans ( : , [4 : 6] ) 
S = 

54 -23 -7 

-16 7 2 

-7 3 1 

(ii) 



54 

-16 

-7 



-23 

7 

3 



-7 

2 

1 



» S*A 
ans = 

1 

0 

0 



» A*S 
ans = 

1 

0 

0 



(iii) 



Both SA and AS are the identity matrix. Hence, S is the inverse of A. 



» inv(A) 
ans = 

54.0000 -23.0000 

-16.0000 7.0000 

-7.0000 3.0000 



-7.0000 

2.0000 

1.0000 



This seems to be the same as S although S-inv(A) may not be exactly zero due to round off. 
The command inv(A) computes the inverse of the matrix A. 



» A = 2*rand(5)-l 
A = 



0.8206 


-0.3435 


-0.5059 


-0 . 8546 


0.5330 


0.5244 


0.2653 


0.9651 


0.2633 


-0.0445 


-0.4751 


0.5128 


0.4453 


0.7694 


-0.5245 


-0.9051 


0.9821 


0.5067 


-0.4546 


-0.4502 


0.4722 


-0.2693 


0.3030 


-0.1272 


-0.2815 
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» R = [A eye(5)] ; 
» rref(R) 



ans = 



Columns 1 


through 7 




1.0000 


0 


0 


0 


1.0000 


0 


0 


0 


1.0000 


0 


0 


0 


0 


0 


0 


Columns 8 


through 10 




1.7333 


-0.2542 


0.3383 


2 . 2202 


0.2963 


-0.8329 


-1.8634 


0.2153 


0.1091 


0.8678 


-0.6154 


-0.5434 


-1.6150 


-0 . 2000 


-1 . 8254 


» S = ans ( : 


,[6:10]) 




S = 






1.6785 


0.1073 


1.7333 


2.0062 


0.1332 


2.2202 


-1.5142 


0.9549 


-1.8634 


0.0574 


0.1234 


0.8678 


-0.7602 


1 . 0249 


-1.6150 


» A*S 






ans = 






1.0000 


0 


0.0000 


0.0000 


1.0000 


0.0000 


0.0000 


0.0000 


1.0000 


0 . 0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


» S*A 






ans = 






1 . 0000 


0.0000 


0.0000 


0.0000 


1.0000 


0.0000 


0.0000 


0.0000 


1.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


» inv(A) - 


S */. This ! 


should be 


ans = 






1.0e-15 * 






-0.4441 


0.1943 


0 


0 


-0.0278 


0 


-0.2220 


-0.1110 


0 


0.0278 


-0.0416 


0.1110 


0.1110 


-0.2220 


0.2220 



0 


0 


1.6785 


0 


0 


2.0062 


0 


0 


-1.5142 


1 . 0000 


0 


0.0574 


0 


1.0000 


-0.7602 



-0.2542 


0.3383 


0.2963 


-0.8329 


0.2153 


0.1091 


-0.6154 


-0.5434 


-0.2000 


-1.8254 



0.0000 


0.0000 


0.0000 


0.0000 


0 . 0000 


0.0000 


1.0000 


0.0000 


0.0000 


1.0000 



0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


1 . 0000 


0.0000 


0.0000 


1.0000 


zero up to 


round off error 



0.3886 


-0.2220 


0.4996 


-0.1110 


-0.1943 


0.1943 


0 


0 


-0.1388 


0.2220 



0.1073 

0.1332 

0.9549 

0.1234 

1.0249 
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2. (i) 

» A = (1/13)* [275; 098; 740]; 

» rref(A) */. part (a) 
ans = 

10 0 
0 1 0 

0 0 1 

» B = inv(A) '/, This will exist since there are no zero rows above 
B = 

-32.0000 20.0000 11.0000 

56.0000 -35.0000 -16.0000 

-63.0000 41.0000 18.0000 

» '/, For part (c): 

» A*B '/, This should be I . 
ans = 

1.0000 0.0000 0.0000 

0 1.0000 0.0000 

0 0.0000 1.0000 

» B*A This should also be I. 
ans = 

1.0000 0.0000 0 

0 1.0000 0.0000 

0 0.0000 1.0000 

» b = 2*rand(3,l) - 1 '/. Choose a random b, with 3 rows, 

b = 

0.0090 

0.0326 

-0.3619 



» rref ( [A b] ) */. Solve Ax=b. 



ans = 

1.0000 


0 


0 


-3.6191 


0 


1 . 0000 


0 


5.1571 


0 


0 


1 . 0000 


-5.7488 


» x = ans(: 


,4); */. Set 


x equal 


to the solution 



» y = inv(A)*b '/, Solve Ax=b using inverses. 

y = 

-3.6191 

5.1571 

-5.7488 

» x-y '/. This should be zero up to round off error, 

ans = 

l.Oe-15 * 

0.4441 

0.8882 

-0.8882 
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(ii) 



» A = [2 -4 5; 0 0 8; 7 -14 0]; 

» rref(A) '/, part (a) 
ans = 

1-2 0 
0 0 1 

0 0 0 

» B = inv(A) V, This will not exist since there are zero rows above. 
Warning: Matrix is singular to working precision. 

B = 

Inf Inf Inf 

Inf Inf Inf 

Inf Inf Inf 

» '/, For part (b) : "singular" means that the matrix is not invertible. 



» A = [1 4 -2 1; 5 1 9 7; 7 4 10 4; 0 7 -7 7]; 
» rref(A) '/, part (a) 
ans = 

10 2 0 
0 1-10 
0 0 0 1 

0 0 0 0 



» B = inv(A) */, This will not exist since there are zero rows above. 
Warning: Matrix is close to singular or badly scaled. 

Results may be inaccurate. RC0ND = 7.178166e-18 

B = 



1.0e+15 * 

-3.6029 -1.8014 

1.8014 0.9007 

1.8014 0.9007 

0.0000 0 



1.8014 1.2867 

-0.9007 -0.6434 

-0.9007 -0.6434 

0.0000 0.0000 



For part (b): From the command rref , we see that A is actually singular. However, MATLAB 
gives us an answer since round off error during the computation makes A seem to be invertible. 
However, MATLAB gives a warning since it could tell “nonsingularity” might be due to round off 
error. 



(iv) 



» A = [1 4 6 1; 5 1 9 7; 7 4 8 4; 0 7 5 7]; 
» rref (A) I, part (a) 
ans = 



10 0 0 
0 10 0 
0 0 10 
0 0 0 1 



» B = inv(A) '/. This will exist since there are no zero rows. 
B = 



-0.1558 -0.0779 

0.0115 -0.1609 

0.2121 0.1061 

-0.1631 0.0851 



0.2208 -0.0260 
0.1133 0.0945 

-0.1061 -0.0758 

-0.0375 0.1025 
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» */. For part (c) : 

» A*B '/, This should be I. 
ans = 

1.0000 0.0000 0.0000 0.0000 

0 1.0000 0.0000 0.0000 

0.0000 0 1.0000 0.0000 

0.0000 0.0000 0 1.0000 

» B*A '/, This should also be I. 
ans = 

1.0000 0.0000 0.0000 0.0000 

0 1.0000 0.0000 0 

0 0 1.0000 0.0000 

0.0000 0.0000 0.0000 1.0000 

» b = 2*rand(4,l) - 1 '/, Choose a random b. 

b = 

0.9733 

- 0.0120 

-0.4677 

-0.8185 



» rref ([A b] ) 
ans = 

1.0000 


'/, Solve 
0 


Ax=b. 

0 


0 


-0.2327 


0 1 


.0000 


0 


0 


-0.1172 


0 


0 1 


.0000 


0 


0.3168 


0 


0 


0 


1.0000 


-0.2260 



» x = ans(:,5); '/. Set x equal to the solution. 

» y = inv(A)*b '/. Solve Ax=b using inverses. 

y = 

-0.2327 

-0.1172 

0.3168 

-0.2260 

» x-y V. This should be zero up to round off error, 

ans = 

1.0e-16 * 

0 

-0.5551 

0 

0.2776 



» A = (-1/56) * [1 2 3 4 5 

0-1 2-1 2 
1002-1 
11-111 
0 0 0 0 4 ]; 



(v) 
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» rref(A) */. part (a) 
ans = 



1 0 
0 1 
0 0 
0 0 
0 0 



0 0 
0 0 
1 0 
0 1 
0 0 



0 

0 

0 

0 

1 



» B = inv(A) 



*/, This will exist since there are no zero rows above. 



B = 


8.0000 - 


40.0000 


-8.0000 


-56.0000 


22.0000 


- 


12.0000 


4 . 0000 


40.0000 


-28.0000 


30.0000 




-8.0000 - 


16.0000 


8.0000 


0 


20.0000 




-4.0000 


20.0000 


-24.0000 


28 . 0000 


-18.0000 




0 


0 


0 


0 


-14.0000 


» 


*/, For part 


(c): 








» 


A*B 7, This should 


be I. 






ans 


= 












1.0000 


0 


0 


0 


0.0000 




0 


1 . 0000 


0.0000 


0 


0.0000 




0 


0 


1.0000 


0 


0 




0 


0.0000 


0 


1.0000 


0.0000 




0 


0 


0 


0 


1 . 0000 


» 


B*A '/, This should 


also be I. 






ans 


= 












1.0000 


0.0000 


0.0000 


0.0000 


0.0000 




0.0000 


1.0000 


0.0000 


0 . 0000 


0 




0 


0 


1.0000 


0 


0.0000 




0 


0.0000 


0.0000 


1 . 0000 


0.0000 




0 


0 


0 


0 


1.0000 



» b = 2*rand(5,l) - 1 '/. Choose a random b. 

b = 

0.5230 

0.5404 

0.6556 

-0.7493 

-0.9683 



» rref ( [A b] ) 


7, Solve 


Ax=b. 








ans = 












1.0000 


0 


0 


0 


0 


-2.0200 


0 1 


.0000 


0 


0 


0 


14.0423 


0 


0 


1.0000 


0 


0 


-26.9510 


0 


0 


0 


1 . 0000 


0 


-10.5699 


0 


0 


0 


0 


1.0000 


13.5557 



» x = ans(:,6); '/, Set x equal to the solution. 
» y = inv(A)*b */, Solve Ax=b using inverses. 



- 2.0200 

14.0423 

-26.9510 

-10.5699 

13.5557 




102 Chapter 1 Systems of Linear Equations and Matrices 



Instructor’s Manual 



» x-y */, This should be zero up to round off error, 

ans = 

1.0e-13 * 

-0.1066 

0.0355 

0 

0 

0.0178 

» A = [ 1 2 -1 7 5 '/, Matrix for (vi) 

0-1 2-3 2 

1031-1 
1114 1 

0 0 0 0 4 ]; 



» rref(A) 
ans = 

1 

0 

0 

0 

0 



'/, part (a) 

0 

1 

0 

0 

0 



3 

-2 

0 

0 

0 



1 

3 

0 

0 

0 



0 

0 

1 

0 

0 



» B = inv(A) % This will not exist since there are zero rows. 
Warning: Matrix is singular to working precision. 

B = 



Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 


Inf 



» '/, This is the same warning as in (ii). 



3. (a) 



» A = round(10*(2*rand(5)-l) ) 
A = 



4 


10 


0 


-5 


-6 


7 


8 


2 


-2 


-7 


3 


-5 


7 


1 


1 


5 


-4 


-2 


-1 


6 


5 


-3 


7 


-4 


-9 


» B = A; 


» B(3, : ) 


= ; 


3*B(1 , : ) 


+ 5*B(2 


1 > • 


B = 


1 


8 


4 


-2 


-6 


0 


2 


-7 


-9 


-4 


3 


34 


-23 


-51 


38 


5 


-7 


-5 


9 


3 


1 


-6 


-10 


-5 


-7 
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ref (B) 










1.0000 


0 


0 


0 


-6.8037 


0 


1.0000 


0 


0 


3.2398 


0 


0 


1.0000 


0 


-4.1014 


0 


0 


0 


1.0000 


4.3544 


0 


0 


0 


0 


0 



Since the bottom row of the reduced matrix is zero, B will not be invertible. 

(b) 

» B= A; 



' - - " 1 
» B(4, : ) 


= 2*B(2, : ) 


- B(1 








B = 












1 


8 4 


-2 


-6 






0 


2 -7 


-9 


-4 






9 


7 -8 


4 


8 






-1 


-4 -18 


-16 


-2 






1 


-6 -10 


-5 


-7 






» rref(B) 












ans = 












1.0000 


0 




0 


0 


-6.6780 


0 


1.0000 




0 


0 


3.0199 


0 


0 


1 . 


0000 


0 


-3.8251 


0 


0 




0 


1.0000 


4.0906 


0 


0 




0 


0 


0 



Again, the bottom row is zero, so B is hot invertible. 

(c) Assume that row i is a linear combination of the other rows. Since one of the valid operations 
in Gaussian Elimination is to add a multiple of one row to another, we may subtract this linear 
combination of the other rows from row i. This will leave a matrix with zeros in row i. By rear- 
ranging the rows, we may put this zero row at the bottom. After continuing the Gaussian Elimi- 
nation, this bottom row will still be zero, so the matrix will be singular. 



4. 



» A = round(10*(2*rand(7)-l) ) 
A = 



4 


2 


1 


-4 


0 


9 


-3 


-2 


9 


-5 


3 


-2 


-2 


-5 


-2 


1 


0 


-7 


-6 


-7 


-7 


0 


-7 


-1 


3 


-9 


8 


6 


-7 


10 


9 


2 


8 


-8 


-1 


2 


-2 


-7 


6 


-1 


-7 


-3 


7 


-7 


-6 


-5 


-7 


-9 


-1 


B = 


A; 












B( : 


,3) = 


2*B(: ,1) 


- B( : 


,2) 






4 


2 


6 


-4 


0 


9 


-3 


-2 


9 


-13 


3 


-2 


-2 


-5 


-2 


1 


-5 


-7 


-6 


-7 


-7 


0 


-7 


7 


3 


-9 


8 


6 


-7 


10 


-24 


2 


8 


-8 


-1 


2 


-2 


6 


6 


-1 


-7 


-3 


7 


-7 


21 


-5 


-7 


-9 


-1 
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» C = A; 



» 


C(: ,4) 


= 


C( : 


: , 1)+C( : ,2) 


-C( : ,3) 


} 




» 


C( : ,6) 


= 


3*C( : , 


2) 










c = 


4 


2 




1 




5 


0 


6 


-3 




-2 


9 




-5 




12 


-2 


27 


-5 




-2 


1 




0 




-1 


-6 


3 


-7 




0 


-7 




-1 




-6 


-9 


-21 


6 




-7 


10 




9 




-6 


8 


30 


-1 




2 


-2 




-7 




7 


-1 


-6 


-3 




7 


-7 




-6 




6 


-7 


-21 


-1 


» 


D = A; 


















» 


D( : ,2) 


= 


3*D( : , 


1) 


» 








» 


D(:,4) 


= 


2*D( : , 


1) 


-D(:, 


2)+4*D( 


: ,3) ; 




» 


D(:,5) 


= 


D(: 


:,2) 


- 


5*D( 


• >3) 






D = 


4 


12 




1 




0 


7 


9 


-3 




-2 


-6 




-5 




-18 


19 


-2 


-5 




-2 


-6 




0 




2 


-6 


-7 


-7 




0 


0 




-1 




-4 


5 


8 


6 




-7 


21 




9 




43 


-66 


-8 


-1 




2 


6 




-7 




-30 


41 


-7 


-3 




7 


21 




-6 




-31 


51 


-9 


-1 


» rref (B) 


















ans 


= 




















1 


0 




2 




0 


0 


0 


0 




0 


1 




-1 




0 


0 


0 


0 




0 


0 




0 




1 


0 


0 


0 




0 


0 




0 




0 


1 


0 


0 




0 


0 




0 




0 


0 


1 


0 




0 


0 




0 




0 


0 


0 


1 




0 


0 




0 




0 


0 


0 


0 


» 


rref (C) 
















ans 


= 




















1 


0 




0 




1 


0 


0 


0 




0 


1 




0 




1 


0 


3 


0 




0 


0 




1 




-1 


0 


0 


0 




0 


0 




0 




0 


1 


0 


0 




0 


0 




0 




0 


0 


0 


1 




0 


0 




0 




0 


0 


0 


0 




0 


0 




0 




0 


0 


0 


0 


» 


rref (D) 
















ans 


= 




















1 


3 




0 




-1 


3 


0 


0 




0 


0 




1 




4 


-5 


0 


0 




0 


0 




0 




0 


0 


1 


0 




0 


0 




0 




0 


0 


0 


1 




0 


0 




0 




0 


0 


0 


0 




0 


0 




0 




0 


0 


0 


0 




0 


0 




0 




0 


0 


0 


0 



Each rref has a row of zeros so the original matrices are not invertible. So we conjecture any 
matrix with some columns equal to linear combinations of other columns will not be invertible. 
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(b) Repeat with your own E. (The example E = A; E( : ,3) = 3*E( : ,2) ; 

E( : ,4) = 2*E( : , 1) + 4*E( : ,3) is interesting.) 

(c) If column j of A is a linear combination of columns preceeding it, then there will be no pivot in 
column j of rref (A), and the non-zero entries in column j of rref (A) may be the coefficients in 
the linear combination which represents the j th column. For instance column 4 in rref (D) has 
-1 in row 1, 4 in row 2 and £)(:, 4) = -1 £>(:, 1) + 4 D(:, 2). (However, for the E suggested in (b), 
rref (E) will not recover the (2,4) coefficients for column 4 since column 3 will not be a pivot 
column.) 

(d) This is the converse of Problem 5, Section 1.7. There we saw that if column j of rref (A) had 
no pivot, then column j of A is a linear combination of the preceding (pivot) columns of A with 
coefficients given by the entries in column j of rref (A). 

5. (a) (i) 



» A = triu( round(10*(2*rand(5)-i)) ); 

» A(2,2) = 0 
A = 

1 -1 -1 -7 5 

0 0 6 3 5 

0 0 -3 2 -4 

0 0 0 -10 -2 

0 0 0 0 4 

» rref (A) '/. part (i) 

ans = 

1-1 000 
0 0 10 0 
0 0 0 1 0 

0 0 0 0 1 

0 0 0 0 0 

Since the bottom row is zero, A is not invertible. This can be repeated 4 more times. In general, 
if A is upper triangular, and there is a zero on the diagonal, A will not be invertible. However, if 
there is not a zero, it will be invertible: 



» A = triu( 
A = 

-6 -2 

0 2 

0 0 

0 0 

0 0 

» rref (A) 
ans = 

1 0 

0 1 

0 0 

0 0 

0 0 



round(10*(2*rand(5)-l)) ) 

10 3 3 

-5 4 3 

-5 6 -9 

0 4 2 

0 0-5 



0 0 0 
0 0 0 
10 0 
0 10 
0 0 1 



This matrix is invertible. 
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(ii) 



» B = inv(A) 
B = 


'/. part 


(ii) 






-0.1667 


-0.1667 


-0.1667 


0.5417 


0.3167 


0 


0 . 5000 


-0.5000 


0.2500 


1.3000 


0 


0 


-0.2000 


0.3000 


0.4800 


0 


0 


0 


0.2500 


0.1000 


0 


0 


0 


0 


-0.2000 



The inverse of an upper triangular matrix is also upper triangular. Also, if the diagonal en- 
tries of A are an then the diagonal entries of A~ l will be 1/a,-,-. Since one of the diagonal 
entries in (i) was 0, the matrix will not be invertible since 1/0 is not defined. 

(iii) Part (iii): To reduce the matrix [A /] to echelon form, we would divide row one by a, row 
two by d, and row three by /. After this step, we would have 1/a, 1/d, and 1// on the diag- 
onal of the right hand side. We would then add multiples of row three to rows one and two, 
which will not change the diagonal on the right. Finally, we would add a multiple of row two 
to row one, which again will not change the diagonal on the right. This will leave the I on 
the left, and A~ x on the right. As predicted in (ii) the diagonal entries of A~ l are the in- 
verses of those in A. Also, as predicted in (i), the first step of this process will fail if a, d, or 
/ is zero. (See the solution to Problem 29.) 



» A = [1 2 3; 4 5 6; 7 8 9]; 

» rref(A) 
ans = 

1 0 -1 

0 12 

0 0 0 

A is not invertible. 

» B = Cl 2 3 4; 5 6 7 8;9 10 11 12; 13 14 15 16]; 

» rref(B) 
axis = 

1 0 - 1-2 

0 12 3 

0 0 0 0 

0 0 0 0 

B is also not invertible. In general, a matrix of this form will not be invertible. If C = (c,y) is an 
n x n matrix with Cij = j + (i — l)n, then c$j = 2 * C 2 j — 1 * cy . This means that the third row is 
2 times the second row minus the first row. 

(c) The assertion that there is a unique nth degree polynomial that fits n + 1 points is the same as 
saying that the coefficient matrix is invertible. 

» x = 2*rand(5, 1)-1 
x = 

-0.2330 

0.0388 

0.6619 

-0.9309 

-0.8931 
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» V= vander(x) 
V = 



0.0029 


-0.0126 


0.0543 


-0.2330 


1.0000 


0.0000 


0.0001 


0.0015 


0.0388 


1.0000 


0.1920 


0.2900 


0.4382 


0.6619 


1.0000 


0.7508 


-0.8066 


0.8665 


-0.9309 


1 . 0000 


0.6361 


-0.7123 


0.7976 


-0.8931 


1.0000 


» inv(V) 










ans = 










8.9238 


-6.5334 


0.7240 


24.5579 


-27.6723 


10.0230 


-9.1141 


1.4612 


10.4447 


-12.8147 


-3.7580 


0.6876 


0.8518 


-13.6375 


15.8561 


-4.7803 


4.1676 


0.1049 


-2.8701 


3.3779 


0.1907 


0.8377 


-0.0054 


0.1314 


-0.1543 



This may be repeated several times. As long as the points in x are distinct, the Vandermonde 
matrix V will be invertible. 

6. (a) Enter Al, A2, ... A5. Then 

» rref(Al) 
ans = 

1 0 0 0 0 

0 10 0 0 

0 0 10 0 

0 0 0 1 0 

0 0 0 0 1 

We see Al is invertible. The same result will come from rref (A3) and rref (A4), so both A3 
and A4 are invertible. 

» rref (A2) 
ans = 

10 3 10 

0 1-230 

0 0 0 0 1 

0 0 0 0 0 

0 0 0 0 0 

» rref(A5) 
ans = 

1-2001 
0 0 10-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 

Both A2 and A 5 are not invertible. 

» rref (A1*A2) 

10 3 10 

0 1-230 

0 0 0 0 1 

0 0 0 0 0 

0 0 0 0 0 



ans 
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» rref (A1*A3) 
ans = 

1 0 0 0 0 
0 10 0 0 
0 0 10 0 

0 0 0 1 0 

0 0 0 0 1 

A1 ■ A2 is not invertible, and A1 ■ A3 is invertible. From the list given A1 • A3, Al ■ A4, and A3 ■ A4 
are invertible, while the others are not: 

» rref (A1*A5) 
ans = 

1-2001 
0 0 10-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 

» rref(A2*A3) 
ans = 

1-2001 
0 0 10-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 

» rref (A2*A4) 
ans = 

1.0000 0 1.2857 -0.0621 0 

0 1.0000 0.2857 0.8509 0 

0000 1.0000 
0 0 0 0 0 

0 0 0 0 0 

» rref (A2*A5) 
ans = 

1-2001 
0 0 1 0-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 

» rref(A3*A5) 
ans = 

1-2001 
0 0 10-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 

» rref(A4*A5) 

1-2001 
0 0 10-2 

0 0 0 1 1 

0 0 0 0 0 

0 0 0 0 0 



ans 
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A product of two matrices will be invertible only if both of the original two matrices are invert- 
ible. In the above, since A2 and A5 were not invertible, any product involving either of these is 
not invertible. 

(b) For AI-A3: 

» inv(Al*A3) - inv(Ai)*inv(A3) 
ans = 



4.8095 - 


20.8571 


-3.3810 


-15.8571 


8.0857 


16.5833 


-7.8690 


16.5238 


-4.7976 


0.7357 


1.0238 


13.0238 


7.0476 


10.0238 


-5.7357 


-5.8452 


-4.7024 


-9.1190 


-3.9881 


2.9500 


0 


0 


0 


0 


0 


inv(Al*A3) 


- inv(A3)*inv(Al) 






1.0e-14 * 










-0 . 0444 


0 


0.2665 


0 


0.0888 


0 


0.1776 


-0.0333 


0.0888 


-0.1332 


0.1332 


-0.1776 


-0.2665 


-0.3553 


-0.2665 


-0.1554 


0.1776 


0.1332 


0.1776 


0.0888 


0 


0 


0 


0 


0 



The second result is 0 to within round off error. Similarly: 



» inv(Al*A4) - inv(Al)*inv(A4) 
ans = 



-16.3571 


13.5714 


3.8571 


10.4286 


6.4286 


3.0000 


9.5714 


9.5000 


18.3571 


-14.1786 


-9.7143 


9.2857 


6.7143 


8.5714 


0.7500 


1.5000 


-1.6429 


-1.5000 


-0.7143 


-1.3929 


0.7857 


-1.4286 


-0.2857 


-2.0000 


0.7857 



» inv(Al*A4) - 
ans = 

1.0e-13 * 
0.0089 -0 

-0.1066 -0 
-0.0488 -0 

0.0155 -0 

0.0006 0 

» inv(A3*A4) - 
ans = 

-74.6667 223 

-531.2333 320 

-435.0667 153 

158.9667 -20 

33.6000 -22 

» inv(A3*A4) - 
ans = 

1.0e-ll * 
-0.1478 0 

0.0284 -0 

-0.0114 -0 

0.0014 0 

0.0092 0 



inv(A4)*inv(Al) 



3730 


0.2087 


0355 


0.3908 


0355 


0.2309 


0644 


-0.0266 


0600 


-0.0377 


inv(A3)*inv(A4) 


3333 


-265.0000 


8333 


-435.5000 


6667 


-228.0000 


3333 


38.0000 


6667 


30.0000 


inv(A4)*inv(A3) 



0284 


-0.2160 


1364 


-0.1307 


0625 


-0.0966 


0135 


0.0199 


0025 


0.0181 



-0.6040 

-0.1243 

-0.1599 

-0.0822 

0.0977 



117.6667 

214.5000 

119.0000 

-22.1667 

-14.6667 



0.1180 

0.0796 

0.0554 

-0.0117 

-0.0107 



0.3730 

0.1599 

0.1332 

0.0389 

-0.0622 



-63.7333 

-48.5667 

-1.9333 

-16.2667 

4.0000 



-0.0142 

0.0092 

0.0014 

- 0.0001 

0.0006 
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From these, we may conjecture that 

(AB)- 1 =0 or {AB)~ l = B~ l A” 1 . 

In fact, this will be true: see theorem 3 in this section. 



» A = [123; 456; 7 8 9]; 
» rref(A) 
ans = 

10-1 
0 12 

0 0 0 



Since the bottom row is all zeros, A isn’t invertible, 
(a) 



(b) 



» format short e ’/, Use scientific notation. 
» f = l.e-5; C = A; C(3,3) = A(3,3) + f; 

» inv(C) 



ans = 

9.9998e+04 

-2.0000e+05 

1.0000e+05 



-2.0000e+05 1.0000e+05 

4.0000e+05 -2.0000e+05 

-2 . 0000e+05 1 . 0000e+05 



C = A; C(3,3) = A(3 ,3) + f; 



» f = l.e-7; 
» inv(C) 
ans = 

1 . 0000e+07 
-2 . 0000e+07 
1 . 0000e+07 



-2.0000e+07 
4. 0000e+07 
-2.0000e+07 



1 . 0000e+07 
-2 . 0000e+07 
1 . 0000e+07 



» f = l.e-10; C = A; C(3,3) = A(3,3) + f; 
» inv(C) 
ans = 

1 . 0000e+10 -2 . 0000e+10 1.0000e+10 

-2 . 0000e+10 4.0000e+10 -2.0000e+10 

1 . 0000e+10 -2 . 0000e+10 1.0000e+10 

(c) The entries in C~ l are roughly the same size as 1//. 

(d) 



» x = [1 ; 1 ; 1] ; 

» f = l.e-5; C = A; C(3,3) = A(3,3) + f; b = [6; 15; 24+f ] ; 
» y = inv(C) * b 

y = 

1.0000e+00 
1 . 0000e+00 
1 . 0000e+00 

» z = x-y 
z = 

0 

0 

4.6566e-10 
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» f = l.e-7; C = A; C(3,3) = A(3,3) + f; b = [6; 15; 24+f ] ; 

» y = inv(C) * b 

y = 

1 . 0000e+00 
1 . 0000e+00 
1 . 0000e+00 

» z = x-y 
z = 

0 

-5 . 9605e-08 
-5 . 9605e-08 

» f = l.e-10; C = A; C(3,3) = A(3,3) + f; b = [6; 15; 24+f]; 

» y = inv(C) * b 

y = 

9 . 9997e-01 
1 . 0001e+00 
9 . 9997e-01 

» z = x-y 
z = 

3 . 0518e-05 
-6 . 1035e-05 
3.0518e-05 

» format '/, Return to standard format for the next problem. 

As / gets smaller, the error term z becomes larger. This means that the closer C is to a nonin- 
vertible matrix, the more error there is in the computation of C~ l . In fact the sum of the ex- 
ponents in / and z is always -15. This is related to the fact that there are about 15 significant 
digits in MATLAB ’s internal computations. 

8. (a) Problem 37: 

» A = [ 0.293 0.014 0.044; 0 0.207 0.010; 0 0.017 0.216]; 

» L = eye(3) -A '/, This is the Leontief matrix. 

L = 

7 . 0700e-01 -1 . 4000e-02 -4.4000e-02 

0 7 . 9300e-01 -1.0000e-02 

0 -1 . 7000e-02 7.8400e-01 

» x = inv(L) * C 13216; 17597; 1786] 
x = 

1.0e+04 * 

1.9305 

2.2225 

0.2760 

Israel needs 19,305 pounds for Agriculture, 22,225 for Manufacturing, and 2,760 for Energy to 
export the given amounts. 
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(b) 

» A = [ .2 .1 .3; .15 .25 .25; .1 .05 0]; 

» L = eye(3) - A 
» format long 

» x = inv(L) * [300000; 200000; 200000] 
x = 

1.0e+05 * 

5.37197626654496 
4.66453674121406 
2.77042446371520 
» format 

This is the same answer as the one to 9(b) in Section 1.3. 

9. If we arrange the message as a sequence of rows, then we need to multiply each row by the encoding 
matrix. Since these are rows, and not columns, the encoding matrix must be on the right. 

Write M for the message and C for the coded message. If encoding the message is done by multiply- 
ing by A then C = M ■ A. Multiplying C by A -1 will decode the message because 

C ■ A- 1 = (M ■ A) -A- 1 = M -(A- A- 1 ) = M ■ I = M. 

» A = [1 2 -3 4 5; -2 -5 8 -8 -9; 12-279; 1106 12; 24-68 11] 

» C = [47 49 -19 257 487 
10 -9 63 137 236 
79 142 -184 372 536 
59 70 -40 332 588] ; 

» M = round(C * inv(A)) '/, use round to get rid of any small error term. 

M = 



1 


18 


5 


27 


25 


15 


21 


27 


8 


1 


22 


9 


14 


7 


27 


6 


21 


14 


27 


27 



»setstr(M + ’A’ - 1) '/. setstr(l:27 + >A’ - 1) is ’ABC...XYZ[> 

ans = '/, So given command prints the message 

ARE[Y '/. but with ’[’ instead of ’ ’ . 

0U[HA 

VING[ 

FUN[[ 

The message decodes to “ARE YOU HAVING FUN”. 
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Section 1.9 




5. 




6. 



1 2 
2 4 



3 -5 




7. 



10 \ 

0 1 

1 0 

0 1/ 



8. 



2 2 1 1 A 
-1465; 



( a d g\ 
b e h I 

cfj) 



10 . 



0 0\ 
0 0 
0 0 / 





/ On O12 • 


' ‘ a lm ^ 




/ 6n 612 ■ 


' * ^lm ^ 


11. Let A = 


021 “22 ’ 


‘ ’ a 2m 


and B = 


621 622 • 


^2m 




V a„i a„2 • 


‘ a nm / 




' 6«i 6„2 • 


bnm f 



Then ( A + B )‘ = 



/ an + bn 021+621 ••• <Jni+6 n i\ 

Ol2 + 612 022 + 622 ' ' ' 0„2 + 6„2 

V ai m + 6lm 02m + b^m ' 1 ' a nm + b nm ) 



= A* +B t 



12. a = 5; /?= 3 

13. aij = aji and 6 tJ - = for 1 < i < n and 1 < j < n. Then, a {j + bij = a j{ + 6j,-. Thus, A + B is 
symmetric. 

14. Since A is symmetric, aj k = a k j for 1 < j < n, 1 < k < n. And since B is symmetric, b ki = b ik for 

n n 

l<i<n,l<k<n. Therefore, ^2, a j k b ki = ^ ~y ik a k j . Thus, (AB) ( = BA. 

k = 1 i = l 

15. Suppose A is m x n. Then A* is n x m. Then AA* is defined and is an m x m matrix. Note that 

n n 

5>a„ = J>*o ijb . That is, the ij th component of A A 1 is equal to the ji th component of AA t . 
k = l k - 1 

Thus, AA l is symemtric. Another proof is that (AA*)* = (+*)*+* = AA\ so AA t is equal to its own 
transpose, hence is symmetric. 

16. If i = j then clearly a, ; - = aji. If i ± j then a,j = 0 and aji = 0. Thus, a {j = Qji . Thus, A is 
symmetric. 
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17. Let U = 



Then U 



fu n u 12 • ■ • «i„\ 

0«22 ' • ' «2n 

V 0 0 • • • 0 u„„ / 



( un 0 • • • 0 ^ 

u 12 u 22 0 - 0 



be upper triangular. 



t _ 



Vuin 



is lower triangular. 



«2 n 



- 0 / 



18. (a) No (b) Yes (c) No (d) Yes 



19. >1* = -A and B' = -5. Then (4 + 5)* = vT + B t = -4 - 5 = -(T + B). Thus, A + B is 
skew-symmetric . 



20. A i = —A ciij = —aji. Elements on the main diagonal are of the form a,,-. a,j = —an. It follows that 
a ii — 0. 

21. (ABy = B t A t = (—B)(—A) = BA. AB is symmetric if and only if ( AB )* = AB. But (AB)* = BA. 
Thus AB is symmetric if and only if A and B commute. 



22. The ij th component of (T + A)/2 = (aij+aji)/ 2 and the ji th component of (T + j4‘)/2 = (aji + aij)/ 2. 
Thus ( A + A *)/ 2 is symmetric. 



23. The ij th component of ( A — A f )/2 = (a,y — aji )/ 2 and the ji th component is (aj,- — a,y)/ 2 = — (a,y — 
aji)/2. Thus (4 — A i )/2 is skew-symmetric. 



24. Let A, B, and C be n x n matrices. Suppose A = B + C where B is symmetric and C is skew- 
symmetric. Then cq j = 6,y + Cij and aji = bji + cj ,• . But b^ = and cji = —Cij. 

Then b , j — |— Cij — aij 

bij — a j i 

Then bij = ( a t; - + aji)/2 and Cjj = (a;j — aji)/2. Note that these solutions are unique. Thus, any 
square matrix can be written in a unique way as the sum of the symmetric matrix (A + A*)/ 2 and the 
skew-symmetric matrix {A — a 1 )/ 2. (This uses the results of Problems 22 and 23.) 



25. A A* 



f an ai2 \ / an a 2 i A _ ( -f aj 2 a na2i + ai 2 a 2 2 \ _ f 1 0 

\ a 2 i a 22 ) \ a i2 a 2 2 / ~~ Vaiia 2 i + a 12 a 22 a\ x + a | 2 ) ~ \0 1 



Thus A is invertible 



and A 1 = A* . 



26. 





(-2 3/2 \ 



= (A- 1 )* 



27. 





(AT 1 = 



2 -3 
-1 2 



= (^ 1 ) < 
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/ 3 0 0 \ 
28. A* = 2 2 0 

VI 2 -1/ 



/ 1/3 -1/3 -1/3 \ / 1/3 0 0 

A- 1 = 0 1/2 1 = -1/3 1/2 0 

V 0 0-1/ \ —1/3 1-1 



/ 105\ 


/ 13/8 - 


-15/8 


5/4 \ 


29. A* = 1 2 5 


(yl- 1 )* = -1/2 


1/2 


0 = (A*)- 1 


V 1 3 1 / 


\ —1/8 


3/8 


-1/4 ) 



= w 1 
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» A = round( 5*(2*rand(4,3)-l) ) 

A = 

-3 4 -5 

-5 -1 -4 

2 0 0 

2 3 2 

» B = round( 5*(2*rand(3,2)-l)) 

B = 

-5 -1 

-1 2 

-4 1 

» (A*B) ’ - A’ * B’ 

??? Error using ==> * 

Inner matrix dimensions must agree. 

» (A*B) ’ - B' * A’ 
ans = 

0 0 0 0 

0 0 0 0 

It is possible that A t B t is not always defined. However (A5) ( is always the same as B t A t . 



» A = [12 3; 2 5 4; 1-1 10]; ’/, This was invertible. 

» inv(A’) 
ans = 

54.0000 -16.0000 -7.0000 

-23.0000 7.0000 3.0000 

-7.0000 2.0000 1.0000 

» inv(A) ’ 
ans = 

54.0000 -16.0000 -7.0000 

-23.0000 7.0000 3.0000 

-7.0000 2.0000 1.0000 

This should be repeated for each of the other matrices. In each case A* is invertible if and only if A is 
invertible. Also, (A*) -1 = (A -1 )*. 

3. (a) 

» A = round( 10*(2*rand(4) -1) ) 

A = 

9 3 5 -3 

7-2-5 3 

14-95 
-8 8 5 10 
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» B = A’ + A 
B = 



18 


10 


6 


-11 


10 


-4 


-1 


11 


6 


-1 


-18 


10 


-11 


11 


10 


20 



For any matrix A, B will be symmetric. 

(b) 



» C = A’ - A 
C = 



0 4 


-4 


-5 


-4 0 


9 


5 


4 -9 


0 


0 


5 -5 


0 


0 


For any A, C will be antisymmetric: C 1 = —C. 


» G = A*A’ 5 


For 


the matrix above. 


G = 


124 23 


-39 


-53 


23 87 


59 


-67 


-39 59 


123 


29 


-53 -67 


29 


253 


» A = round( 


10*(2*rand(3,4) -1) ) */, A non-square matrix 


A = 


-3 4 


-9 


-5 


-5 5 


3 


-1 


10 3 


8 


5 



» G = A*A’ 
G = 



131 


13 


-115 


13 


60 


-16 


-115 


-16 


198 



For any matrix A, G will be symmetric. 

(d) If the ij entry of A is a,j then the ij entry of A* is ciji. The ij entry of B is = a,ij + a,, which 
is the same as bji = a,-,- + o,y, so B is symmetric. 

The ij entry of C is a,y — a j ,• , while the ij entry of C* is a,-,- — a 8 j = — c,y, so C* — —C. 

The ij entry of G will be gij = Yk a >kajk ■ The ji entry of G will be Yk a jk<iik which is the same 
as gij . Hence G is symmetric. 

4. (a) In problem 2 from section 1.7, it was shown that the solution of Ax = 0 produces all vectors x 
that are perpendicular to the rows of A. Since the columns of A are the same as the rows of A t , the 
solution of yFx = 0 will produces all vectors x that are perpendicular to the colums of A. 

(b) (i) 



» A = [2 0 1; 0 2 1; 1 1 1; -1 1 1; 111]; 
» rref(A’) 
ans = 



1.0000 


0 


0.5000 


0 


0.5000 


0 


1.0000 


0.5000 


0 


0.5000 


0 


0 


0 


1.0000 


0 
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The solutions have X3 and X5 arbitrary, and x\ = —.5x3 — .5xs, x 2 — — .5x3 — .5x5, and x 4 = 0. 

(ii) 

» A = [2 4 5; 0 5 7; 7 8 0; 7 0 4; 9 1 1] ; 

» rref(A’) 



ans = 



1.0000 


0 0 


5.8462 


5.3846 


0 


1.0000 0 


-3.6044 


-3.7033 


0 


0 1.0000 


-0.6703 


-0.2527 


The solutions have 
3.7033x5, and X 3 = 


x 4 and X 5 arbitrary, and x x — —5.8462x4 — 5.3846x 5 , x 2 = 3.6044x 4 -f 
: .6703X4 + .2527X5. 



» A = rand(5,3) 



A = 


0.2661 


0.3841 


0.9410 




0.0907 


0.2771 


0.0501 




0.9478 


0.9138 


0.7615 




0.0737 


0.5297 


0.7702 




0.5007 


0.4644 


0.8278 




>> rrefCA 1 ) 
ans = 


1.0000 


0 


0 


1.0319 0.5701 


0 


1.0000 


0 


1.7249 -0.4801 


0 


0 


1.0000 


-0.3771 0.4142 


The solutions have 


x 4 and X5 


arbitrary, and x x = —1.0319x4 — 0.5701x 5 , x 2 = —1.7249x4-} 


0.4801x5, and X3 = 


: 0.3771X4 


- 0.4142x5. 





5. 



» A = 2*rand(4)-l 
A = 

-0.7493 0.2591 

-0.9683 0.4724 

0.3769 0.4508 

0.7365 0.9989 



0.7771 0.0265 
-0.5336 0.1822 
-0.3874 0.6920 
-0.2980 -0.1758 



» Q = orth(A) 



0.5071 0.2864 

0.6553 0.4839 

-0.2551 0.3366 

-0.4984 0.7553 



0.7859 

-0.5678 

-0.1922 

0.1515 



-0.2077 

0.1184 

-0.8858 

0.3977 



(a) 



» x = 2*rand(4, 1)-1 
x = 

0.6830 

-0.4614 

-0.1692 

0.0746 
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» y = 2*rand(4, 1)-1 

y = 

-0 . 0642 
-0.4256 
-0.6433 
-0.6926 

» s = x’ * y '/, The scalar product of x and y. 
s = 

0.2097 

» r = ( Q*x )’ * (Q*y) '/ The scalar product of Qx and Qy. 

s = 

0.2097 

» format short e 
» s-r 
ans = 

5 . 5511e-17 

Since Q is orthogonal, x • y is the same as (Qx) • (Qy). 

(b) If the above steps are repeated, even if A is complex, the inner product of x and y will always be 
the same as that of Qx and Qy. (For complex A this depends on the fact that A! — (Ay . Also 
you should reverse the x, y variables as x ■ y = y' * x for complex vectors.) 



» x = Q(:,l); y = Q(:,2); '/, Let x be the 1st column, and y the 2nd. 
» sqrt(x’ * x) */, the length of the 1st column, 
ans = 

1.0000 



» x’ * y '/ The inner product of the 1st and 2nd column, 
ans = 

5 . 5Slle-17 



(d) 



This is zero up to round off and the same results follow for all the other columns. 



» inv(Q) 
ans = 

0.5071 0.6553 

0.2864 0.4839 

0.7859 -0.5678 

-0.2078 0.1184 



-0.2551 -0.4984 

0.3366 0.7553 

-0.1923 0.1515 

-0.8858 0.3976 



Q 1 and Q' are the same for any orthogonal matrix. 

(e) To show that x • y = Qx • Qy, we will use Q -1 = Q‘. (If Q is complex replace * with '.) 

Qx- Qy = (Qx)*Qy = x*Q‘Qy = x‘Q _1 Qy = x*/y = x*y = x • y. 

The ith column of Q is Qe; where e, is a vector with a 1 in the zth position and zeros elsewhere. 
This means that the inner product of the ith column of Q with the yth column can be written as 
Qej ■ Qej. Using (b) we have that Qe,- ■ Qej = e,- • ej , which is I if i = j and is 0 if i ^ j. The 
statement in (c) follows immediately. 
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Section 1.10 

1. RiTlR 2 is an elementary matrix 2. f ?2 + Ri is an elementary matrix 

3. Since two operations RxT^-R-i, R 2 — * R 2 + R\ are needed, this is not an elementary matrix 

4. 2R 2 is an elementary matrix 

5. Two operations needed: R\ — *■ 3i?i, R 2 — ► 3J?2, so not an elementary matrix 

6. is an elementary matrix 

7. Two operations, R,\7±R 2 , R 2 J±R 3 , so not an elementary matrix 



8. Two operations needed, R 2 —* R 2 + 2f?i, R 3 —> R 3 + 3i?i, so not an elementary matrix 

9. i ?2 + 2/?i is an elementary matrix 10. R 4 + R 2 is an elementary matrix 



11. Two operations needed. R 2 — ► R 2 + R \ , R 4 — ► i ?4 + R$, so not an elementary matrix 

12. R\ — R 2 is an elementary matrix 



13. 



17. 



0 0 1 
0 1 0 
10 0 , 



21 . 



25. 



29. 



33. 



14. 



18. 




15. 



19. 




16. 



20. 




22. 



26. 



30. 



34. 



1 0 
-2 1 




0 1 0 
1 0 0 
0 0 1 



23. 



1 2 
0 1 



31. 



35. 




24. 



32. 



36. 



0 1 
1 0 



1 0 0\ 




(-1 0 0\ 




/I 


-2 0 


-3 10 


27. 


0 10 


28. 


0 


1 0 


0 0 1/ 




V 001 ; 




\o 


0 1 




38 . 



/I 0 -1 0\ 
0 1 0 0 
0 0 10 
Voo 01/ 



/I 0 0 


-5\ 




/I 0 0 0\ 


0 1 0 


0 


40. 


f 0100 


0 0 1 


0 


,0310 


Vo 0 0 


1 / 




\0 0 0 1/ 
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41 . 




R\ —*■ ^R\ 

R 2 — ► R 2 3i?i 




R 2 2R 2 
Ri — ► Ri — ^R 2 




A = 



f20UlO\ A 0\/ll/2\ 

\0 1/ 1/ \0 1/2/ \0 \) 




43. 



(lll\ 
0 2 3 
\ 5 5 1 / 



R$ — ► /?3 — 5Hi 
«2 - \ R 2 
R\ — > i?i — R 2 



/ 1 0 1/2 \ 
0 1 3/2 

\ 0 0 -4/ 



•R 3 — > -yH 3 
/?2 ^2 — 1 . 5 i ?3 

R\ — ► R\ + O.5H3 



(!!!) 

\001/ 



/100\ /100\ / 1 1 0\ /10 0\ / 1 0 — 1/2\ / 1 0 0\ 
A = I 0 1 0 ) J020J I 0 1 0 I [Ol 0 01 0 013/2 

\5 0 1/ \001/ V 001 / \0 0 — 4/ \00 1/ \00 l) 



R 3 —* —R 3 
fi 2 —* R 2 — R 3 
Ri -* Ri + hR 3 




/3 0 0\ /100\ / 1 2/3 0 \ / 1 0 0\ / 1 0 0\ / 1 0 — 1/3\ 
A= 010 020 0 1001 001101 0 

\001/\001/\0 0 1/ \o 0 -1/ \oo 17 \0 0 l ) 



45. 



( 



0 -1 
0 1 

1 0 




R\ ^Rz 
R 3 — ► R 3 + ^2 




#3 — ► — /?3 

^2 — ► /?2 + 2 R 3 

/?1 — + — /?3 



(!!!) 

Vo 0 iy 



/001\ /I 00\ /10 0\ / 1 0 1\ / 1 0 o\ 
A = I 0 1 0 I I 0 10 01 0 010 01 -1 

\l 0 0/ \o -1 1/ \0 0 -1/ \0 0 1/ \0 0 1/ 




rtl - 5 R 1 

R 3 -* R 3 - 3R1 I J J J 

R 3 R 3 + R 2 \ Q Q _4 



R3^-\R 3 . 

Ri — ► R\ — 2 R 3 



\0 0 1 



/2 0 0 \/l 00 \/l 00 \/l 0 0 \ / 1 0 0 \ / 1 0 2 \ 
> 1 = 010 010 0 10 01 0 011010 
\0 0 1 / V 3 0 1 / \0 -1 1 / \0 0 -4 / V 0 0 1 / Vo 0 1 / 
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/2 0 0 0\ 


/I 0 0 0\ 




"s* 

H— ‘ 

0 

0 

0 


/I 0 0 0\ 


[0100 


[ 0 3 0 0 




0 1 0 0 


0 10 0 


0 0 10 


0 0 10 




0 0-40 


0 0 10 


Vo 0 0 1 / 


Vo 0 0 1 / 




Voo 01 / 


Vo 0 0 5 / 



/2 1 0 0\ 


Ri -+ O. 5 R 1 


/iO 


-1/4 ON, 


R 3 — 0.5i ? 3 


0 2 10 


R 2 — ► 0.5R 2 


0 1 


1/2 0 


r 2 —*■ R 2 — O. 5 R 3 


0 0 2 1 


R\ — ► R\ — 0.5i?2 


0 0 


2 1 J 


Ri —+ R 1 + \R 3 


Vo 0 0 2/ 


Vo 0 


0 2/ 





( 10 0 1/8 \ 
0 1 0 -1/4 | 
0 0 1 1/2 1 
0 0 0 2 / 



/?4 — ► 0 . 5 i ?4 
R3 — ► i?3 — O.5H4 
R 2 — ► ^2 4" 0.25/?4 
iij — ► Hi — 



( 1 0 0 0\ 
0 10 0 
0 0 10 
0 0 0 1 / 



/2 0 0 0 \ 


/I 0 0 0 \ 


/I 1/2 0 0 \ 


/I 0 0 0 \ 




/!0 


0 0 \ 


0 10 0 


0200 ] 


0 


10 0 


0 10 0 




0 11/2 0 


0 0 10 


0 0 10 


0 


0 10 


0 0 2 0 




0 0 


10 


Vo 0 0 1 / 


Vo 0 0 1 / 


Vo 


0 0 1 / 


Vooo 1 / 




Vo 0 


0 1 / 



/I 0 -1/4 0\ 

0 1 0 0 

0 0 10 

Vo 0 01 / 



A = 



/I 0 0 0\ 
0 10 0 
0 0 10 
Vo 0 0 2 / 



/1 0 0 0\ 
0 10 0 
0 0 11/2 
Vooo 1 / 



/1 0 0 0\ 
0 1 0 -1/4 
0 0 1 0 
Vooo 1 / 



/I 0 0 l/8\ 
0 10 0 1 
0 0 1 0 I 
Vooo 1 / 



49. ac ^ 0 implies a ^ 0 and c ^ 0. Row reducing A to / 2 by back elimination via 




c _1 it 2 

— bR 2 




a 1 




shows / 2 = (a 1 /?i)(/?i — 6i? 2 )(c 1 j R 2 )-A- So solving for A gives A = (ci? 2 )(Ri + 6/2 2 )(a7?i ) , or 



A = 




50. adf 0 implies a, d and / are nonzero. Row reducing A to I 3 by back elimination gives 

f-'Ri 



\ 



( a b c 
Ode 

00 f) 






R 2 

Ri 



R? — el?3 



/a b 0\ 
OdO 

Rl ~ cR3 V 0 0 1 / 



Ri 



d-'Ri 

— ► R\ — bR 2 



(a 0 0\ 
010 
\° 0 1 / 



l Ri 
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So taking all the inverse operations 



applied to la (in the opposite order) shows 



A = 






/I 0 0\ /I b 0\ fa 0 0\ 
( 0 <f 0 I 010 010 

\00 1 / \00 1 / \00 1 / 



51. Let U = («J“ 1 1 i?i)(uJ 2 1 ^ 2 ) • • '( u Zr%Rn)A. Then U is an n x n matrix with l’s down the diagonal 
and 0’s below it. Note that U is row equivalent to I and hence, by theorem 4, can be written as a 
product of elementary matrices. If U — E\ E- 2 ■ ■ ■ Ek, where each Ei is an elementary matrix, then 
A = ( u nn R n ) • ■ (u 22 ^ 2 )(«ii^i)£'i -^2 • • • Ek ■ By theorem 4, A is invertible. 



52. By problem 51, A is invertible. As in problem 51, let U = (u nfliX u 22^) ' ' ' ( u nnRn)A. Since U 

has l’s down the diagonal and 0’s below it, when row reducing U to I we need only add multiples 
of a row to those rows above it. Hence, we can write U as a product of upper triangular elementary 
matrices: U = E\E 2 ■ ■ ■ Ek- Note that E~ l is upper triangular for each i. Show that the product 
of two upper triangular matrices is upper triangular. Then deduce A -1 = E k l • • • W^i) 

(u 22 1 f? 2 ) ■ • -( u nnRn) is upper triangular. 

53. A* is upper triangular, so (A*) -1 is upper triangular by the result of problem 52. But (A*) -1 = (A -1 ) 4 , 
so (A -1 )* is upper triangular, which means that A -1 = [(A -1 ) 4 ]* is lower triangular. 



n n 

54. PijA = C where c rs = If r = i, then Ci , = = aj s since pik is 1 if k — j and 0 

*= 1 i=l 

otherwise. If r = j, then Cj s = a, 3 since pjk is 1 if k = i and 0 otherwise. If r ^ i and r ^ j, then 
c TS = a rs . Hence, PijA is the matrix obtained by permuting the i th and j th rows. 



n 

55. AijA = B where b r , = ^ ,. If r ^ j, then b r , = a r , since a' rk is 1 if At = r* and 0 otherwise. If 

it=i 

n 

r = j, then bj, = = ca,-» + aj s since a'j k is c if k = i, 1 if k = j, and 0 otherwise. Hence, 

k = l 

AijA is the matrix obtained from A by multiplying the i th row by c and adding it to the j th row. 



m 

56. Mi A = B where b r , = If r = i, then 6,j = ca,-, since m,-* is c if A? = i and 0 otherwise. 

t=i 

If r ^ i, then b rs = a rs since m r k is 1 if Ar = r and 0 otherwise. It follows that M, A is the matrix 
obtained from A by multiplying the i th row by c. 



H 2 



— 2 R\ 





1 2 
0 0 



57. 



1 2 

2 4 
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/ 1 — 3 3 \ / 1 — 3 3\ 

61. 0 -3 1 — [ 0 1 -1/3 

\1 02/ V° 0 0/ 

/ 1 0 0 \ /I 00\ / 1 0 0 \ / 1 — 3 3\ 

4= 010 0 -3 0 010 0 1 -1/3 

\1 0 1/ \0 01/ \ 0 3 1 / \0 0 0/ 



/ 100 \ / 1 0 0 \ 

62. 2 3 0 J -+ | 0 1 0 I 

\ — 1 4 0 / \0 0 0j 

/ 1 0 0 \ / 100 \ / 1 0 0 \ / 100 \ / 1 0 0 \ 
4=210)1 010030010010 

\0 0 1/ \ — 1 0 1/ \001/ \0 4 1/ \0 0 0/ 
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MATLAB 1.10 

1. (a) 



» A = 


round (10 


* (2*rand(4)- 


D) 


A = 


-6 


9 


-9 -10 




-9 


-2 


-9 -2 




4 


0 


1 -9 




4 


7 


3 -2 




w 


» F=eye(4) ; F(3,3) = 


4; '/, R3 — > 4 R3, i.e. mult ip! 




» F*A 
ans = 
-6 


9 -9 


-10 




-9 


-2 -9 


-2 




16 


0 4 


-36 




4 


7 3 


-2 


(ii) 




» F=eye(4) ; F(l,2) = 


-3; */. R1 — > R1 - 3 R2, i.e. : 




» F*A 
ans = 
21 


15 18 


-4 




-9 


-2 -9 


-2 




4 


0 1 


-9 




4 


7 3 


-2 


(iii) 




» F=eye(4) ; F([l 4], 


:) = F( [4, 1] . ; ; '/. interchange 




» F*A 
ans = 

4 


7 3 


-2 




-9 


-2 -9 


-2 




4 


0 1 


-9 




-6 


9 -9 


-10 


(b) 


» 


F=eye(4) ; 


; F(3,3) = 4; 


*/. Part (i): R3 --> 4 R3 


» 


inv(F) 




’/. This is R3 — > 1/4 R3, i.e. 


ans = 

1.0000 


0 


0 0 




0 


1.0000 


0 0 




0 


0 0 


.2500 0 




0 


0 


0 1.0000 
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(ii) 



2. (a) 



» F=eye(4) ; F(l,2) = -3; 
» inv(F) 
ans = 

13 0 0 

0 10 0 

0 0 10 

0 0 0 1 



*/. R1 --> Ri - 3 R2, 

*/. This is Rl — > Rl + 3 R2, i.e. add 3R2 to Rl . 



» F=eye(4) ; F(Cl 4],:) = 
» inv(F) 
ans = 

0 0 0 1 

0 10 0 

0 0 10 

10 0 0 



F([4,l],:); */. Part (iii): interchange 1 and 4. 

'/, This also interchanges 1 and 4. 



In each case, the inverse of F represents a row operation that is the reverse of the original 
operation. 



» A = [7 2 3; -1 0 4; 2 1 1] 
A = 

7 2 3 

-10 4 

2 11 
» B = A; 

» c = -B(2 , 1)/B(1 , 1) ; 

» FI = eye(3); Fl(2,l) = c; 

» B = F1*B ; 

» F = FI; 



'/, Store this matrix in B. 

'/, Compute the multiplier to eliminate B(2,l). 

'/, Generate the elementary matrix which 
'/. subtracts c*row 1 from row 2. 

'/, Apply the matrix FI to B. 

*/, F will be the product of the elemetary matrices. 



» c = -B(3 , 1)/B(1 , 1) ; 

» F2 = eye(3) ; F2(3,l) = c; 
» B = F2*B 
B = 

7.0000 2.0000 3.0000 

0 0.2857 4.4286 

0 0.4286 0.1429 

» F = F2*F ; 



'/, Compute the multiplier to eliminate B(3,l). 
'/, Form elementary matrix for R3-cRl 
'/, Finish column 1 elimination 



*/, F accumulates the product of Fi’s. 



» c = -B(3,2)/B(2,2); 

» F3 = eye(3) ; F3(3,2) = c; 
» B = F3*B 
B = 

7.0000 2.0000 3.0000 

0 0.2857 4.4286 

0 0 -6.5000 

» F = F3*F ; 



*/, Next, eliminate B(3,2) in column 2. 

’/. This will finish forward elimination 
*/, Apply the matrix to B. 



'/, Keep track of F. 



» c = 1/B(3 , 3) ; 

» F4 = eye(3); F4(3,3) = c; 
» B = F4*B; 

» F = F4*F ; 



*/, Start backelimination, divide row 3 by B(3,3). 

'/, Normalize R3 to start with 1. 

'/. Keep track of F. 
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» c = -B(2,3)/B(3,3) ; 

» F5 = eye(3) ; F5(2,3) = c; 
» B = F5*B; 

» F = F5*F ; 



'/, Now backeliminate B(2,3). 

'/ Apply the matrix to B. 

'/, Keep track ol F. 



» c = — B ( 1 ,3)/B(3,3) ; 

» F6 = eye(3) ; F6(l,3) = c; 
» B = F6*B 
B = 

7.0000 2.0000 0 

0 0.2857 0 

0 0 1.0000 



'/. Next, backeliminate B(l,3) to finish column 3. 
*/. Apply the matrix to B. 



» 


F = F6*F ; 




'/, Keep track of F. 


» 


c = 1/B(2,2) 


> 


*/, Next, divide row 2 by B(2,2). 


» 


F7 = eye(3) ; 


F7(2 , 2) = c; 




» 


B = F7*B 




'/. Normalize row 2 to start with 1 


B = 


7 2 


0 






0 1 


0 






0 0 


1 




» 


F = F7*F ; 




*/. Keep track of F. 



» c = -B(l,2)/B(2,2); 

» F8 = eye(3) ; F8(l,2) = c; 
» B = F8*B 
B = 

7 0 0 

0 10 

0 0 1 

» F = F8*F ; 



'/, Next, backeliminate B(l,2). 

'/, Apply the matrix to B, finish column 2. 



*/. Keep track of F. 



» c = 1/B(1 , 1) ; 

» F9 = eye(3) ; F9(l,l) = c; 
» B = F9*B 
B = 

10 0 
0 10 

0 0 1 

» F = F9*F 
F = 

0.3077 -0.0769 -0.6154 

-0.6923 -0.0769 2.3846 

0.0769 0.2308 -0.1538 



’/, Finish by dividing row 1 by B(l,l). 
'/, Normalizing Rl gives rref(A) 



'/. Form the final product of Fi’s. 



(b) 



F*A 






1.0000 


0.0000 


0.0000 


0.0000 


1.0000 


0.0000 


0 


0 


1.0000 


A*F 






1.0000 


0 


0.0000 


0 


1 . 0000 


0.0000 


0.0000 


0.0000 


1.0000 
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Both FA and AF are the identity, so F = A~' . FA = I 3 since F = (F9)(F8)(F7)...(F2)(Fl) 
and all the steps showing the row echelon form of A is I 3 can be achieved by one left multiplica- 
tion by F instead of the individual multiplications by FI, F2, .... 



» D = inv(Fl)*inv(F2)*inv(F3)*inv(F4)*inv(F5)*inv(F6)* . . . 
inv (F7 ) * inv (F8 ) * inv (F9 ) 

D = 



7 . 0000 


2.0000 


3 . 0000 


- 1.0000 


0 


4.0000 


2.0000 


1.0000 


1.0000 


D is the inverse of F 


= A~\ 


so it is the same as A. ( D = F~ l since the inverse of a product is 



the product of the inverses in the opposite order.) 



» A = [023; 114; 241]; 
» B = A; 



'/. Store this matrix in B. 

'/, Since B(l,l) is zero, we must first 
'/, interchange it with another row, so that we 
'/, do not get a divide by zero error. This is 
'/, the only difference between the steps here 
'/, and those in part (a) . 



» FI = eye(3) ; Fl([l,2],:) 

» B = F1*B 
B = 

114 
0 2 3 

2 4 1 

» F = FI; '/, F will be the product of the elemetary matrices. 

» c = -B(3, 1)/B(1 , 1) ; '/, Compute the multiplier. 

>> F2 = eye(3); F2(3,l) = c '/, R3-cRl will eliminate B(3,l). 



= FI ([2,1],:); */. interchange 1 and 2. 
'/. Apply the matrix FI to B. 



» B = F2*B 
B = 

114 
0 2 3 

0 2-7 



'/, Apply the matrix F2 to B. 



» F = F2*F ; 

» c = 1/B( 1 , 1) ; 

» F3 = eye(3) ; F3(l,l) = c; 
» B = F3*B 
B = 

114 
0 2 3 

0 2-7 

» F = F3*F; 

» c = -B(l,2)/B(2,2); 

» F4 = eye(3) ; F4(l,2) = c; 
» B = F4*B 



’/, F is now the product of F2*F1 . 
7, Next, divide row 1 by B(l,l). 

'/, Apply the matrix to B. 

’/, Keep track of F. 

'/, Next, eliminate B(l,2). 

'/, Apply the matrix to B. 
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B = 

1.0000 0 2.5000 

0 2.0000 3.0000 

0 2.0000 -7.0000 




» F = F4*F ; 


'/, Keep track of F. 


» c = -B(3,2)/B(2,2) ; 

» F5 = eye(3) ; F5(3,2) = c; 


i Next, eliminate B(3,2). 


» B = F5*B 
B = 

1.0000 0 2.5000 

0 2.0000 3.0000 

0 0 -10.0000 


'/, Apply the matrix to B . 


» F = F5*F ; 


'/, Keep track of F. 


» c = 1/B(2,2) ; 

» F6 = eye(3) ; F6(2,2) = c; 


'/, Next, divide row 2 by B(2,2). 


» B = F6*B 
B = 

1.0000 0 2.5000 

0 1.0000 1.5000 

0 0 -10.0000 


'/, Apply the matrix to B. 


» F = F6*F ; 


'/ Keep track of F. 


» c = -B(l,3)/B(3,3) ; 

» F7 = eye(3) ; F7(l,3) = c; 


'/, Next, eliminate B(l,3). 


» B = F7*B 
B = 

1.0000 0 0 

0 1.0000 1.5000 

0 0 -10.0000 


*/, Apply the matrix to B. 


» F = F7*F ; 


'/, Keep track of F. 


» c = -B(2,3)/B(3,3) ; 

» F8 = eye(3) ; F8(2,3) = c; 


'/, Next, eliminate B(2,3). 


» B = F8*B 
B = 

10 0 
0 10 

0 0 -10 


Apply the matrix to B. 


» F = F8*F; 


1, Keep track of F. 


» c = 1/B(3,3) ; 

» F9 = eye(3) ; F9(3,3) = c; 


'/, finally, divide row 3 by B(3,3) 


» B = F9*B 
B = 

10 0 
0 10 

0 0 1 


7. Apply the matrix to B. 


» F = F9*F 
F = 

-0.7500 0.5000 0.2500 

0.3500 -0.3000 0.1500 

0.1000 0.2000 -0.1000 


'/, Look at final form of F. 



As in (b) and (c), we find that F*A=A*F=I, so F is the inverse of A, and that D=inv(Fl)*inv(F2) 
* . . . * inv(F9) is the inverse of F, so it is the same as A. 




130 Chapter 1 Systems of Linear Equations and Matrices 



Instructor’s Manual 



3. (a) 

» A = C 1 2 3; 1 1 7; 2 4 5] ; 

» U = A; */. Store A in U; after reduction U will be "echelon" form 

» c = -U(2,l)/U(l,l) '/. Eliminate U(2,l). 

c = 

-1 

» FI = eye(3) ; Fl(2,l)=c 
FI = 

10 0 
-110 
0 0 1 

» U = F1*U '/. Apply FI to U, to eliminate. 

» c = -U(3, 1)/U(1,1) '/. Eliminate U(3,l). 

c = 

-2 

» F2 = eye (3) ; F2(3,l) = c 
F2 = 

10 0 

0 1 0 

-2 0 1 

» U = F2*U '/. Apply F2 to U. 

U = 

12 3 

0-14 
0 0-1 

Note: U is now upper triangular. 

» F = F2*F1 '/, F is the product of the elementary matrices. 

F = 

10 0 

-1 1 0 

-2 0 1 

(b) 

» L = inv(Fl)*inv(F2) 

L = 

10 0 
110 
2 0 1 

The matrix L is lower triangular. For each of the entries below the diagonal in L, we see that 
it is the same value as — c, where c was the multiplier used to eliminate the same entry in A. In 
fact the entries are just the entries in the inverses inv(Fl), inv(F2). Since we can recover FI, F2 
from L (move down column 1 for each successive Fi). 
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(c) 

» L*U */, This is the same as A. 
ans = 

12 3 

117 
2 4 5 



Since F is the product of the elementary matrices that reduce A to U , we know that FA = U. 
Since L = F ~ l , we have LU - LF A = F~*FA = A. 



» A = [6 2 7 3; 8 10 1 4; 10 7 6 
» U = A; */. 
» c = -U(2, 1)/U(1 , 1) •/. 



c = 



-1.3333 

» FI = eye(4) ; Fl(2,l) 
FI = 

1.0000 0 

-1.3333 1.0000 

0 0 

0 0 



= c 



0 

0 

1.0000 

0 



» U = F1*U; 

» c = -U(3, 1)/U(1 , 1) 
c = 

-1.6667 

» F2 = eye(4) ; F2(3,l) = c 



1.0000 


0 


0 


0 


1 . 0000 


0 


-1.6667 


0 


1.0000 


0 


0 


0 



» U = F2*U ; 

» c = -U(4, 1)/U(1 , 1) */. 

c = 

-0.6667 



» F3 = eye(4) ; F3(4,l) = c 
F3 = 



1.0000 

0 

0 

-0.6667 



0 0 

1.0000 0 

0 1.0000 
0 0 



» U = F3*U 
U = 

6.0000 2.0000 

0 7.3333 

0 3.6667 

0 6.6667 



7.0000 

-8.3333 

-5.6667 

4.3333 



8; 4895]; 

Store A in U, and work with U. 
Eliminate U(2,l). 



0 

0 

0 

1.0000 



Eliminate U(3,l). 



0 

0 

0 

1.0000 



Eliminate U(4,l). This finishes column 1. 



0 

0 

0 

1 . 0000 



3 . 0000 
0 

3.0000 

3.0000 
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» c = -U(3,2)/U(2 ,2) 


'/, Eliminate U(3,2), moving on to column 2. 


c = 






-0.5000 






» F4 = eye(4) 


; F4(3,2) = c 




F4 = 






1.0000 


0 0 


0 


0 


1 . 0000 0 


0 


0 


0.5000 1.0000 


0 


0 


0 0 


1.0000 


» U = F4*U; 






» c = -U(4,2)/U(2,2) 


'/, Eliminate U(4,2), finish column 2. 


c = 






-0.9091 






» F5 = eye(4) 


; F5(4,2) = c 




F5 = 






1.0000 


0 0 


0 


0 


1.0000 0 


0 


0 


0 1.0000 


0 


0 


0.9091 0 


1.0000 


» U = F5*U 






U = 






6.0000 


2.0000 7.0000 


3.0000 


0 


7.3333 -8.3333 


0 


0 


0 -1.5000 


3.0000 


0 


0 11.9091 


3 . 0000 


» c = -U(4,3)/U(3,3) 


*/, Eliminate U(4,3), moving on to column 3. 


c = 






7.9394 






» F6 = eye (4) 


; F6(4,3) = c 




F6 = 






1.0000 


0 0 


0 


0 


1.0000 0 


0 


0 


0 1.0000 


0 


0 


0 7.9394 


1.0000 


» U = F6*U 






U = 






6.0000 


2.0000 7.0000 


3 . 0000 


0 


7.3333 -8.3333 


0 


0 


0 -1.5000 


3.0000 


0 


0 0 


26.8182 


» F = F6*F5*F4*F3*F2*F1 


'/, The product of the F’s 


F = 






1.0000 


0 0 


0 


-1.3333 


1.0000 0 


0 


- 1.0000 


-0.5000 1.0000 


0 


-7.3939 


-4.8788 7.9394 


1.0000 


» L = inv(Fl)*inv(F2)*inv(F3)*inv(F4)*inv(F5)*inv(F6) '/, This is the inverse of F. 
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1.0000 0 

1.3333 1.0000 

1.6667 0.5000 

0.6667 0.9091 



0 0 

0 0 

1.0000 0 

-7.9394 1.0000 



If you print inv(Fi) you see its non- zero entry below diagonal is just the negative of correspond- 
ing entry in Fi and is equal to corresponding entry in L. 



» L*U 
ans = 

6 2 7 3 

8 10 1 4 

10 7 6 8 

4 8 9 5 



*/. As in (c), this is the same as A. 



As in (b), we see that both L and F are lower triangular, and that LU — A. As in (c), we each 
entry in L is the negative of the multiplier used to eliminate the corresponding entry in U , or the 
non-zero entries in inv(Fl),. . .,inv(F6). 
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Section 1.11 

In problems 1-8, the matrix L is constructed from the identity matrix in the following way: For each 
row operation on A of the form R\ —+ Ri—lcRj, put a A; in the i,j position of the identity matrix. 



1. A = 



1 2 \ R-2 ~ y ^2 — 3R\ 

3 4 



(i 



Thus L = 



1 0 

3 1 



2 A=lll)= u - ThusL = (j; 



3. A = 



-1 5 

6 3 



— *■ R-2 + f — 



- 1 b ) = 
0 33 ) 



0. Thus L = 



1 0 
-6 1 



4. A 



Ri 




R 2 — * R 2 ~ 2 K 1 




R3 — + R3 — 3 i?i 



0 0-3, 



r _ ion / 1 4 6\ /I 00' 

3 9 2 I 0 -9 -9 I = U. Thus L = ( 2 10 

3 fl 



1 4 6' 

0 -9 -9 
0 -10 -13. 



5. .4= ^ 435 ) 




R-2 — y R-2 ~ 2R\ 




= u. 



( 100 \ 

Thus L = 2 1 0 

V 1 0 1 / 



6. A = 




jF?2 — ► R 2 ~~ 2Hj 



H 3 



o _ 60 / 3 9-2 

3 0-2112 



0 




R3 —> R3 — ^Ri 



'3 9-2' 

0 -21 12 

.0 -6 f 



1. A - 



1 1 2 -1 4\ 

0-1 5 8 

2 3 14 

1 6 4/ 



Vi 




Ri — ► Ri — R\ 



1 2-1 4\ 

0-1 5 8 

0-1 3-4 

.0 -3 7 0 / 



/I 2-1 


4 \ 


Z 1 


2 -1 4\ 




1^3 + R 3 R 2 / 0 — 1 5 


8 


— ► /?4 — 3R 2 I 0 


-15 8 


R± — *■ R± — 4/?3 


, 1 n n _9 


-12 


f 0 


0 -2 -12 




VO -3 7 


0 / 


Vo 


0 -8 -24/ 





/i 


2 -1 


4 \ 




/I 0 0 0\ 


In 


-1 


5 


8 


= U. Thus L = 


r 0 1 0 0 | 


0 


0 -2 


-12 


i 2 1 1 0 


Vo 


0 


0 


24/ 




VI 3 4 1/ 
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8. A = 




R 2 — ■ + R-2 — 2R\ 



2 3-1 6 ' 

0 1 4-11 

-2 5-2 0 

0 -4 5 2/ 



R3 —* R3 + Ri 



(2 3-1 6 \ 

0 1 4-11 

0 8-3 6 

\o — 4 5 2/ 





/ 2 


3 


-1 


6 \ 




/ 2 3 - 1 6 \ 


/ 2 3 -1 


6 \ 




0 


1 


4 


-n 


— ► /?4 -f* 4i?2 


0 1 4-11 Rt ->Ri+ ffis 


0 1 4 


-H 


K 3 — ► i?3 — 8 H 2 | 


0 


0 • 


-35 


94 




n n -as 94 I 


0 0 -35 


94 




\0 


-4 


5 


2 / 




VoO 21 -42/ 


Vo 0 0 


72 / 
5 7 



U. Thus L = 



1 0 0 O' 

2 10 0 
-18 10 
0-4-f 1, 



9. From problem 1 we have A = LU where L — ( ^ and U = f q 2 ) • The system Ty = b, i.e. , 



3l) (y 1 ) = ( 4 ) y * e ^ s e 9 uat ions 2/1 = —2 and 3yi + 2/2 = 4. Solving we get y 2 = 10. Now, 
from Ux = y, i.e., ^ = ( lo) We °^ ta * n Xl + ^2 = — ^ anc * — ^x 2 = 10. Backsolving 

we get xi = 8, x 2 = —5. The solution is ^ . 



10. From problem 2 we have A = LU where L = ( g and U = f q 3 )• The system Ly = b, 



i.e., ^g (i/ 1 ) = ( 4 ) yields the equations j/i = — 1, 



y 2 — 4. Now, from Ux = y, i.e., 



1 2 \ (x x 
03j \x 2 

4 
3 



and x 2 — — . The solution is j 3 4 



^ j we obtain x\ + 2x 2 = —1 and 3 x 2 = 4. Backsolving we get x\ = — ^ 

-11 

4 
3 



11. From problem 3 we have A = L/7 where L = 
10 \(yi\ ( 0 



_g ^ j and U = ( q 33 ) • The system Ly = b, i.e., 



-6 1 / V 2/2 
from Ux = y, i.e. 



5 

-1 5 
0 33 



25 5 

Xl = 33 ’ 372 = 33 ' The solution is 



yields the equations y\ = 0, — 6j/x + 2/2 —3. Solving we get 2/2 = 5. Now, 

^ ^ = ^ 5 ^ we obtain — x\ + 5x 2 = 0, 33x 2 = 5. Backsolving we get 

" (!) 

' 33 / 



1 00\ /I 4 6 \ 

12. From problem 4 we have A = £17 where L = | 2 1 0 ),{/=( 0 — 9 — 9 I . The system Ly = b, i.e., 

VO 0-3/ 



1 , 




-1 



10 



2/2 1 = 1 7 yields the equations yi = -1, 2yi + t/ 2 = 7, 3j/i + — y 2 + 2/3 = 2. Solving 

3^1/ \ 2/3 / V 2/ 9 




